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The Diffraction of X-Rays by Liquid Na-K Alloy in a Magnetic Field 


C. W. Heaps, The Rice Institute 


(Received June 27, 1935) 


Because magnetoresistance is found in the liquid Na-K alloy there is reason to expect that a 
magnetic field might modify the x-ray diffraction curve of the liquid. Experiment indicates 
that the diffraction is not changed by as much as two percent in a field of 2700 gauss. A test 
made to detect magnetostriction in the alloy indicates no volume change per unit volume 
greater than 3X 1077 in a field of 7800 gauss. The cause of magnetoresistance in the liquid does 
not, therefore, appear to lie in changes of molecular grouping caused by a magnetic field. 


HE liquid alloy of sodium and potassium 
has recently been found to increase its 
resistance when placed in a magnetic field.' It 
also shows a fairly large Hall effect. The dis- 
covery of these galvanomagnetic effects in a 
liquid is a matter of some importance; it would 
now seem that a molecular lattice structure is 
not necessary for the production of the effects, 
although this structure has occasionally been 
suggested as fundamental in the phenomena. 
The empirical fact has been, heretofore, that 
strongly crystalline substances of low degree of 
symmetry (e.g., Bi, Sb, graphite) have the 
largest values of magnetoresistance. The Som- 
merfeld theory as developed by Houston’ as- 
sumes the resistance of a metal to be due to 
deviations of the lattice from symmetry, such 
deviations being caused by impurities or thermal 
motions. If a magnetic field in any way affects 
the regularity of lattice structure we might ex- 
pect magnetoresistance to appear. 

There is at present a considerable body of 
evidence in favor of the view that liquids possess 
characteristic molecular structures. G. W. Stew- 

1 Fakidow and Kikoin, Physik. Zeits. d. Sowjetunion 3, 


381 (1933); J. E. Armstrong, Phys. Rev. 47, 391 (1935). 
2W. V. Houston, Zeits. f. Physik 48, 449 (1928). 
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art’s idea*® is that large numbers of molecules 
may be cooperating at any one instant of time to 
form temporary lattice arrangements in the 
liquid. Although this viewpoint is not necessary 
in order to explain x-ray diffraction curves for 
liquids there are a number of reasons why it is a 
plausible hypothesis. If it is accepted we might 
expect a magnetic field to affect either the num- 
ber of molecules which at any one time are co- 
operating, the time of their cooperation, or the 
orientation (crystalline axis) of the group. 
Hence, following the general lines of Houston's 
theory, the field might be expected to change the 
resistance of the liquid. If the field produces 
any of the suggested effects on molecular group- 
ing, the x-ray diffraction curve of the liquid 
would be modified, provided the effects are large 
enough. It seemed worth while, therefore, to 
investigate the effect of a magnetic field on the 
x-ray diffraction curve of the liquid Na-K alloy. 
EXPERIMENT 

The x-ray spectrometer was equipped with 

Soller slits and an ionization chamber made of 


seamless steel tubing filled with argon at about 
20 atmospheres pressure. A celluloid window 


3G. W. Stewart, Rev. Mod. Phys. 2, 116 (1930). 
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Fic. 1. Effect of filter on the x-ray spectrum of Mo. Curve 
I, no filter; curve II, with ZrQsz filter. 


admitted the rays. The ionization current was 
measured by the use of a Western Electric 
D-96475 electrometer tube in an amplifying 
circuit as described by Harnwell and Van 
Voorhis.* The sensitivity was adjusted to be 
about 40,000 divisions per volt and an input 
resistance of 10'° ohms was used. The circuit 
was carefully shielded and readings were not 
taken until about two hours after the filament 
current had been turned on. Under these condi- 
tions the galvanometer drift was negligible and 
random fluctuations were seldom more than 
two mm. 

The x-rays were produced by a Coolidge tube 
with water-cooled molybdenum target. A po- 
tential of about 30 kv was used. The tube was 
enclosed in a lead box provided with a small hole 
through which the x-ray beam could pass. A 
filter of ZrO, placed over this hole produced 
approximate monochromatization of the beam. 
Fig. 1 shows the spectrum of the rays obtained 
with a calcite crystal, with and without the 
filter. 

The liquid alloy contained approximately 35 
percent sodium, this proportion having been 
found by Armstrong to give an alloy with fairly 
large magnetoresistance. The cell for holding the 
alloy was about 1 mm thick. It was made by 


4G. P. Harnwell and S. N. Van Voorhis, Rev. Sci. Inst. 
5, 244 (1934). 
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Fic. 2. Scattering by Na-K alloy. For the upper curve 
obtained without filter, the right-hand scale of ordinates 
is to be used. For the lower curve, obtained with filter, 
the left-hand scale is used. The experimental points all 
refer to the lower curve. 


bending a thin glass rod into suitable contour 
and cementing thin mica plates to both sides of 
this framework. The scattering of the empty 
cell was small; however it was carefully measured 
for the various angles and subtracted from the 
total scattering after the cell was filled. This 
difference was assumed to be the scattering of the 
alioy. To fill the cell it was placed mouth down- 
ward in an evacuated tube. The alloy, in a side 
bulb, was now decanted into the tube. When air 
was admitted the alloy was forced up into the 
cell. This cell was mounted on the spectrometer 
table between the poles of an electromagnet. 
Fig. 2 shows the scattering curves obtained 
with and without the filter. The points are 
averages of a number of observations taken with 
random distribution. A magnetic field of 2700 
gauss produced no measurable effect on the form 
of either curve. In testing the effect of the mag- 
netic field the following procedure was adopted. 
The galvanometer zero was first noted, then a 
lead shutter was opened to admit the x-ray beam 
through the Soller slits to the specimen. The 
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scattered beam passed through the second set of 
slits into the ionization chamber and produced 
a definite deflection of the galvanometer. The 
magnetic field was now switched on and careful 
observation made to see whether the deflection 
changed. Then the field was switched off, the 
shutter closed, and the galvanometer zero 
checked again. In no case could any effect of the 
magnetic field be detected. Tests were made at 
1° intervals for both curves of Fig. 2, a field of 
2700 gauss being used. An effect as large as two 
percent could have been detected in the case of 
the largest deflections. 


DISCUSSION OF RESULTS 


The diffraction of x-rays by Na-K alloys has 
been investigated previously by K. Banerjee,5 
who used the Mo Ka radiation and a photo- 
graphic method. Banerjee finds a maximum at 
11° 28’. In the present work the maximum is at 
about 9.5°. The scattering of the points at the 
maximum in Fig. 2 appears to be rather too 
large for high accuracy; however, repeated trials 
indicated a definitely smaller angle for the maxi- 
mum than was found by Banerjee. The reason 
for the discrepancy is not evident. It may be 
mentioned, however, that thin, flattened, glass 
tubes, such as Banerjee used to contain his 
alloy, could not be used in the present work 
because their scattering was comparable to that 
of the alloy. Using the Bragg formula the spacing 


> K. Banerjee, Ind. J. Phys. 3, 399 (1928). 
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of diffracting centers in the liquid turns out to be 
4.3 Angstrom units, according to the present 
data, instead of 3.5 as given by Banerjee. 

The absence of any effect of a magnetic field 
indicates that no appreciable modification of 
molecular groupings is produced. Hence the 
magnetoresistance of the alloy cannot be ex- 
plained on these grounds. The possibility re- 
mains that certain structural changes may occur 
which are too small to be detected by the x-ray 
method. Ferromagnetic crystals when magnet- 
ized give an unchanged x-ray spectrum although 
definite structural changes occur, as is evidenced 
by the presence of magnetostriction. 

An attempt was made to detect magneto- 
striction in the Na-K alloy. A quantity of the 
material under oil in a glass bulb was placed be- 
between the poles of an electromagnet and sub- 
jected to a field of 7800 gauss. Fhe oil meniscus 
in a fine capillary tube connected to the bulb was 
observed under a microscope. It was estimated 
that if any magnetostriction occurred the volume 
change per unit volume must have been less than 
3 10-", the limit of accuracy of the experiment. 
This is a smaller volume change than is found 
to exist in ferromagnetic materials. 

The conclusion appears to be that in this 
alloy a magnetic field does not produce a resist- 
ance change by modification of any molecular 
groupings or spacings in the liquid. The main 
cause of magnetoresistance is thus to be sought 
in the Lorentz forces which modify the electron 
free paths. 
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The Transmutation Functions for Some Cases of Deuteron-Induced Radioactivity’ 


ErNEst O. LAWRENCE, Epwin MCMILLAN AND R. L. THORNTON, Radiation Laboratory, Department of Physics, 
University of California, Berkeley 
(Received July 1, 1935) 


When nuclei are bombarded with deuterons, a commonly 
observed type of reaction is that resulting in proton emis- 
sion and the conversion of the nuclei into isotopes of mass 
number one unit greater. Many such reactions give beta- 
radioactive products. The variation of the transmutation 
cross section with deuteron energy has been studied for 
the radioactivity thus produced in Na, Al, Si and Cu, with 
deuteron energies ranging from 0 to 3.6 MV. In all these 
cases it is found that the excitation curve is too flat to be 


explained by the Gamow theory, while it is fitted very well 
by a theoretical curve derived by Oppenheimer and 
Phillips. This is based on the idea that the nucleus can cap- 
ture the neutron from the deuteron without requiring the 
penetration of the deuteron through the nuclear Coulomb 
barrier. The form of the theoretical curve depends on the 
binding energy of the deuteron; a value between 2.0 and 
2.4 MV is required by the experimental results. 
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INTRODUCTION 


N THE original investigations of the radio- 

activity induced in sodium? and aluminum? by 
deuteron bombardment, observations of the 
voltage excitation curves up to a deuteron energy 
of 1.9 MV were made. In this range the experi- 
mental curves can be fitted reasonably well by the 
Gamow‘ nuclear penetration function, which has 
been found to describe so well many other nuclear 
reactions produced by a-particle and proton bom- 
bardment. It was noticed, however, in the case 
of aluminum, for which the experimental data 
were more detailed, that the experimental points 
corresponding to the lowest deuteron energies 
lie above the theoretical curve when it is made 
to fit the high energy points. The true significance 
of this was not realized until the experiment was 
repeated with a deuteron beam of much greater 
energy, extending the excitation curve for 
aluminum up to 3.4 MY. 

It was then found that the differential excita- 
tion function for aluminum beyond 2 MY\, in- 
stead of rising with a rapidly increasing slope as 
predicted by the Gamow theory, goes up at a 
much more modest rate. The part of the curve 
between 2 and 3.3 MV is in fact a straight line 
when plotted against the range of the deuterons. 
Curves of similar shape were found for the sodium 
reaction, and some others that will be described 
later. (See Figs. 2—5.) 

These observations indicated that some new 
description of the mechanism of these reactions 
was necessary. If one should attempt to explain 
the experimental results while retaining the usual 
picture of a simple penetration of the bombarding 
particle through the nuclear potential barrier, it 
would be necessary to suppose either that a con- 
glomeration of unresolved resonance levels hap- 
pened to give this particular shape to the curve, 
or else that the deuterons were going over the top 
of the barrier at a remarkably low energy. The 

! This work was first reported at the Washington meeting 
of the National Academy of Sciences. Abstract in Science 
81, 421 (1935). 

* Ernest O. Lawrence, Phys. Rev. 47, 17 (1935). 

’ Edwin McMillan and Ernest O. Lawrence, Phys. Rev. 
47, 343 (1935). 

4*Gamow, Atomic Nuclei and Radioactivity (Oxford Uni- 
versity Press, 1931). This is based on the wave-mechanical 
idea of the ‘‘diffraction”’ of particle waves through potential 
barriers, which was first applied to a nuclear problem 
(a-radioactivity) by Gurney and Condon and independ- 
ently by Gamow. 
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first of these suppositions is made highly im- 
probable by the facts that the curves are smooth 
and monotonic, and are of similar shape for 
several different reactions. The second would 
require an abnormally large nuclear radius. 
Duncanson and Miller’ find 9 MV for the height 
of the Al potential barrier to a-particles; one 
would expect the height for deuterons to be half 
as great, and therefore beyond the energy range 
covered in these experiments. 

At this point Professor Oppenheimer became 
interested in the problem from a_ theoretical 
standpoint, and found that the experimental re- 
sults are not at all mysterious, but are just what 
must be expected for reactions of the type in- 
vestigated when the small stability of the deu- 
teron is taken into account. A full account of the 
theory is given by Professor Oppenheimer and 
M. Phillips in the paper immediately following 
this one. We shall return to a discussion of the 
mechanism of reaction involved and its applica- 
tion in this case after devoting some attention to 


the experiments. 


REACTIONS STUDIED 


The formation of radioactive bodies from 
sodium and aluminum by deuteron bombard- 
ment has been shown? * to occur almost certainly 


according to the reactions: 


Na**+ H*?-Na*t+H)!, (1) 
Al’ +H?—-AP8+H!. (2) 
Recently H. W. Newson,® working in this 


laboratory, has found that an active product of 
2.8 hours half-life, isotopic with silicon, is made 
from silicon under deuteron bombardment. The 
reaction concerned is presumably: 


Si*°+-H?—Si"" + H!. (3) 


Finally, we have observed that copper with 
deuterons gives an activity of about 12 hours 
half-life, which is probably the same substance as 
the active copper isotope produced by neutron 
bombardment, whose half-life is given as 10 
hours by Amaldi et al.? We also found a shorter 


®°W.E. Duncanson and H. Miller, Proc. Roy. Soc. A146, 
413 (1934). 

®°H. W. Newson, Bull. Am. Phys. Soc. 10, No. 4, p. 8 
(1935). 

7 E. Amaldi, O. D’Agostino, E. Fermi, B. Pontecorvo, F. 
Rasetti and E. Segré, Proc. Roy. Soc. A149, 522 (1935). 
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period, but cannot be sure that it arises from 
copper without further work, taking greater care 
to eliminate the ubiquitous carbon and nitrogen 
contamination effects. Only the long period activ- 
ity will be considered in this paper. It is weak 
compared to that of sodium; after equal exposure 
toa 3.6-MV deuteron beam a copper target shows 
about 1/200 the activity of a sodium target. 
Nevertheless it is very much larger than one 
would expect for an element of such high atomic 
number (29) without some modification of the 
ordinary theory of nuclear penetration. The reac- 
tion of formation is probably: 


Cu"+H?—-Cu"t!+HI!, (4) 


where m may be 63 or 65 depending on which of 
the two stable copper isotopes is the parent sub- 
stance. 

It will be noticed that all these reactions are 
of the same type, in which a deuteron is captured 
and a proton emitted, converting the bombarded 
nucleus into its next higher isotope. This sort of 
process can also be described by saying that the 
nucleus snatches a neutron away from the deu- 
teron, leaving the proton free to escape, and it is 
in terms of this description that the theory of 
Oppenheimer and Phillips is to be understood. 
Many such reactions with radioactive products 
are now known. The four listed above® were 
chosen for the present study on purely mechanical 
grounds, in that the substances concerned can be 
obtained in the form of very thin uniform foils, 
as required by the experimental technique used. 


EXPERIMENTAL METHOD 


In order to obtain a voltage-yield curve for a 
reaction giving a radioactive product, one meas- 
ures the intensity of activation of a number of 
samples exposed to bombardment by particles of 
different energies. The relative bombarding cur- 


5 One other case of particular interest has recently been 
examined by M. C. Henderson (Bull. Am. Phys. Soc. 10, 
No. 4, p. 6 (1935)). Magnesium with deuterons gives two 
active bodies, by the reactions: 


Mg**+H?—+Mg??+H!, (5) 
Mg?*+H?—>Na"‘+ Het. (6) 


Making use of the widely different decay periods, he ob- 
tained the separate excitation curves of the two reactions, 
and found that reaction (5) follows a curve of the type 
described in this paper, while reaction (6) follows a Gamow 
curve. Theoretically one expects just this sort of behavior, 
since in reaction (6) the deuteron as a whole must penetrate 
into the Mg** nucleus. 
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Fic. 1. Arrangement for bombarding targets. 


rent and exposure time must be accurately known 
for each sample, a requirement which is most 
easily and satisfactorily met by activating them 
simultaneously. This was done by using as a 
target a stack of thin foils of the substance being 
investigated, through which the activating deu- 
teron beam was sent. The energy of the deuterons 
traversing any particular foil is then determined 
by their initial energy and the stopping power of 
the preceding foils, and the measured activities 
give directly the differential excitation curve. 
An apparatus of the type developed by Lawr- 
ence and Livingston® was used to produce a beam 
of high speed deuterons.'? The beam was allowed 
to emerge from the vacuum through a thin 
aluminum window (8 mm air equivalent stopping 
power), mounted on a brass grid near the end of a 
tube built into the side of the apparatus, as 
shown in Fig. 1. The target was placed next to 
the window, and the space between them was 
evacuated during the bombardment. Deuteron 
currents between 1 and 2 microamperes and ex- 


posure times from 5 minutes to 1.5 hours were 


used. 
MEASUREMENTS 


The initial energy of the deuteron beam was de- 
duced from its range in air, determined visually 


* Ernest O. Lawrence and M. Stanley Livingston, Phys. 
Rev. 45, 608 (1934). 

‘0 Since we shall have many occasions in the future to 
refer to this apparatus, we feel that it should have a name. 
The term ‘‘magnetic resonance accelerator’’ is suggested. 
In this, the last two words imply the essential principle of 
operation, while the first word is added to distinguish it 
from the apparatus of Sloan and Lawrence (Phys. Rev. 38, 
2021 (1931)), which can be called a ‘‘linear resonance ac- 
celerator.”’ The word ‘“‘cyclotron,” of obvious derivation, 
has come to be used as a sort of laboratory slang for the 
magnetic device. 
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ic. 2. Differential excitation curve of the radioactivity 
induced in aluminum (Z=13) by deuteron bombardment. 
Circles, experimental points. Solid lines, Oppenheimer- 
Phillips functions with the values of J given by the num- 
bers, in MV. Dashed line, Gamow function. 


by observing the end of the streak of blue lumin- 
escence that it produces. Since the maximum air 
path available (see Fig. 1) was less than the 
range of the beam, aluminum foils were put next 
to the window, of such a thickness as to reduce 
the range to less than 4.5 cm. The stopping power 
of these foils was determined by weighing. The 
error in the range measurements made by this 
method was certainly not greater than 5 mm. 

The stopping powers of the foils which were to 
be activated were determined by means of a 
polonium source mounted on a micrometer screw 
and a Wynn-Williams counter, the displacement 
of the counts-against-distance curve on interpos- 
ing a foil being taken as its stopping power. An 
error of about 0.2 mm is possible in these 
measurements. 

The activities of the foils were observed by a 
quartz-fiber electroscope in an ionization chamber 
provided with a thin aluminum window to admit 
B-particles. All the activities were large compared 
to the natural background, and these determina- 
tions should be precise to within one or two per- 
cent. The time of making each reading was re- 
corded, making it possible to correct for the de- 
cay of the active products. In plotting the results, 
points with ordinates proportional to the cor- 
rected activities are put at abscissas correspond- 
ing to the centers of the foils. This method of 
plotting would be precise only if the curves were 
straight lines, but it is easy to see that the error 
introduced by the curvature is negligible with the 
spacing of points used. 

Details of the various cases investigated are 
given below: 
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Fic. 3. Differential excitation curve of the sodium radio- 


activity (Z=11). Notation same as in Fig. 2. The ordinate 
of the highest energy experimental point is probably too 
low. 


Aluminum (Fig. 2) 


Range of deuteron beam in air (at 20°C) =11 
cm. Stopping power of individual foils (pure 
aluminum) = 11.5 mm. Since the active product 
has a half-life of about 156 sec., it decayed by a 
factor of about 7 in the 7 min. needed to measure 
the activities of all the foils, but the half-life is 
well enough known that this can be corrected for 
with very little error. The part of the curve below 
4.5 cm deuteron range is taken from the results 
of reference 3, which were obtained with a 
shorter initial range, the ordinate scales being 
adjusted to make the two curves fit at 4.5 cm. 
The abscissas of the low-range points taken in 
this way are obviously more precise than those 
obtained by the passage of a long-range beam 
through a great thickness of foils. 


Sodium and silicon (Figs. 3 and 4) 


Range of beam =12 cm. Stopping power of 
foils =14.5 mm. The foils used were of mica, 
whose composition is given by the formula He 
(KK, Na) Als Sis Oi. This substance after exposure 
to deuterons contains several active bodies; there 
are some short periods that practically vanish 
within an hour, as well as the 2.8-hour silicon 
product, and the 15.5- and 13.5-hour sodium and 
potassium products. Of the latter two, the sodium 
activity is much the stronger, even though the 
amount of sodium is less than that of potassium 
in ordinary micas (Na/K~10-30 percent). A 
subsidiary experiment shows this. Targets of 
NaCl and KCl were given equal exposures 
(checked by the equality of the 40 min. chlorine 
activities induced), and observed 20 hours later. 
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Fic. 4. Differential excitation curve of the silicon radio- 
activity (Z=14). Notation same as in Fig. 2. The ordinate 
of the highest energy experimental point is probably too 
low. 


The NaCl was found to have 70 times the activity 
of the KCI. Therefore in mica the greater suscept- 
ibility of the sodium to activation far outweighs 
its smaller abundance, and we can assume with 
little error that the whole long period activity 
arises from sodium. 

The activities of the foils were measured 2 and 
20 hours after activation, in order to disentangle 
the silicon and sodium effects. The total activities 
(sum of all foils) observed for silicon and sodium, 
corrected for decay back to the end of the activa- 
tion (1.5 hours at 1 microampere) were 6 and 1.8 
divisions,‘sec., respectively, corresponding to 
about 3 and 0.9 microcuries of active substances. 
The saturation activities after a very long expo- 
sure would be 8 and 13 microcuries. 

The measured values of the activities of the 
highest energy foil are probably too low. This 
error arises from the forward recoil momentum 
imparted to the activated nuclei by the colliding 
deuterons, which, as shown in reference (3), 
drives the active material forward in the direction 
of the beam with an average range of the order of 
a millimeter air equivalent. As a result, an ap- 
preciable part of the activity of the first foil is 
transferred on,to the second one. For the second 
and succeeding foils the loss in activity on one 
side is nearly compensated by the gain on the 
other from the preceding one. Aluminum is not 
subject to this error, since even the first foil is 
preceded by the aluminum window. Copper has 
such a heavy nucleus that the recoil range is 
small, and the error produced is negligible. 
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Fic. 5. Differential excitation curve of the long-period cop- 
per radioactivity (Z=29). Notation same as in Fig. 2. 


Copper (Fig. 5) 


Range of beam=10 cm and 12 cm. Stopping 
power of foils (pure copper) = 9.8 mm. This curve 
was obtained in two sections. The part from 0 to 
10 cm was done at a time when the apparatus 
was not set to give a beam of the highest energy, 
and the extension up to 12 cm was made later, 
the two ordinate scales being adjusted to make 
the curves fit together at 9.5 cm. 

The activities were measured 3 hours after 
bombardment, when only the long period copper 
activity remained. The total activity (all foils), 
expressed as that after a very long exposure to 1 
microampere of 12 cm deuterons, was 14 divi- 
sions/sec. or about 7 microcuries. 

Since the stopping power of copper relative to 
that of air varies with the velocity of the particles, 
the air equivalent thickness of the foils measured 
with polonium a-particles is slightly in error for 
some parts of the deuteron range. The necessary 
correction, which increases the effective stopping 
power for the highest energy and decreases it for 
the lowest, was made according to the values 
given by Mano." The corresponding correction 
for aluminum and mica is negligible. 


DISCUSSION OF THE THEORY AND [ts FIT WITH 
EXPERIMENT 


The essence of the theoretical explanation of 
these results given by Oppenheimer and Phillips 
lies in the fact that reactions of the above type 
can occur without requiring the penetration of a 
charged particle through the nuclear Coulomb 
potential barrier which opposes the entry of the 


(4, Mano, J. de phys. et rad. 5, 628 (1934). 
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deuteron. The mechanism postulated can be 
described in the following way: Consider the 
deuteron as built of a proton and neutron, with 
the binding energy /. Because of the wave nature 
of the particles, this structure is not confined 
within a sharply limited volume, but there is a 
finite chance that the neutron or proton may be 
found some distance away from the center of 
mass. This chance decreases rapidly with the 
distance. When the deuteron is projected against 
a nucleus, the proton is held back by the Coulomb 
field. This distorts the probability distribution of 
the proton in such a way that for a given distance 
of the deuteron center of mass from the nucleus, 
the chance of finding it in the center of the 
nucleus is very much reduced. On the other hand, 
the neutron is immune to the action of the electro- 
static field, and its probability distribution can 
overlap the center of the nucleus; this distribu- 
tion is in fact distorted by the forces involved in 
the deceleration of the neutron, in such a way as 
to increase its density on the side toward the 
nucleus. Therefore the neutron has a relatively 
larger chance of being within the nucleus, and 
while there of undergoing a reaction in which it 
becomes bound, forming a new nucleus. The 
extension of the neutron distribution varies with 
I, becoming greater as J decreases. It is the rela- 
tively small value of J (~2 MV) that causes the 
slow energy variation of this type of reaction, 
by flattening out the neutron probability dis- 
tribution. 

The actual shape of the function giving the 
variation of the neutron-penetration with energy 
depends on the value of J. Since this most 
probably lies between 2.0 and 2.4 MV (see ap- 
pendix), the curves for these two cases have been 
compared with the experimental values. They 
are plotted as solid lines on Figs. 2—5. On Fig. 2 
(aluminum) the curve for J/=1.5 MV is plotted 
also. The ordinates of the theoretical curves are 
adjusted to make them fit at the highest experi- 
mental points, except in the cases of sodium and 
silicon where the measured values of these are 
almost certainly too low because of recoil effects 
(see the special discussion of these cases). 

The Gamow function: 

(1 ve sreZe2/ he (7) 


2 Eq. (1), reference 3, which purports to give this func- 
tion, contains a typographical error. 
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where Z =atomic member, v=deuteron velocity, 
is also plotted on the figures for comparison. In 
this connection it should be emphasized that the 
power of v which multiplies the exponential, al- 
though usually given as 1 2°, is not definitely 
determinable without a rather detailed knowl- 
edge of the forces acting during the reaction. 
Depending on these, it may vary between 1 and 
1 v°. We have used 1,v°, which gives the flattest 
curve, in order to show that even in this extreme 
case the function is too steep to fit the experi- 
mental results. 

On the other hand, the results of Oppenheimer 
and Phillips show a very satisfactory agreement 
with the observed values. Consider first the case 
of aluminum. This is the most trustworthy ex- 
perimentally; with pure aluminum foils the activ- 
ity is very large compared to that of any possible 
contamination, and there is no error in the 
highest point caused by recoil. It is seen that the 
Oppenheimer-Phillips function with J=2.0 fits 
extremely well, while with /=1.5 and J=2.4 it 
is off by more than the experimental error. 

The fit for sodium and silicon is also best for 
I= 2.0. Here the possibility of experimental error 
is somewhat greater, and the apparent deviation 
in the case of silicon between 1.5 and 2.5 MV 
may not be real. The results for these two ele- 
ments are of chief interest in showing that this 
type of excitation curve is not a characteristic of 
any particular reaction, but of a class of reac- 
tions. The comparative values for sodium and 
silicon show also that the steepness of the curves 
increases with Z, as it should. 

For copper the fit seems to be best with J = 2.4. 
Here, because the greatest variation of the func- 
tion occurs in a narrower voltage range, the 
different theoretical curves are not as widely 
separated as in the other cases and the apparent 
better agreement with J = 2.4 may not be real. We 
can give an explanation for it, which, however, 
we do not wish to emphasize unduly. Oppen- 
heimer and Phillips point out that the theory is 
subject to small modifications whose evaluation 
would require an intimate knowledge of the 
nuclear structure, and which would be in such a 
direction as to increase the value of J giving the 
best fit; the magnitude of these modifications 
would be least for elements of large Z. The result 
of this would be to make the form of the theory 


= 
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used here require a larger value of J for copper 
(Z=29) than for the other elements studied 
(Z=11, 13 and 14). 

Such refinements are hardly justified at pres- 
ent. We can, however, make some conclusions 
that we feel are definitely beyond the possible 
uncertainties of both the experimental results 
and the theory. They are first, that the theory 
proposed by Oppenheimer and Phillips is ade- 
quate to explain the observed excitation curves 
while the Gamow theory is not, and second that 
the value of the deuteron binding energy conso- 
nant with this explanation is 2 MV or slightly 
larger. 

The acceptance of this reaction mechanism as 
correct leads to another very important conclu- 


sion. Just as the variation of the transmutation 
yield with energy for a given atomic number is 
slower than expected, likewise the variation with 
atomic number for a given energy is slower. This 
has the consequence that it will be possible, with 
the available deuteron energies, to produce 
nuclear reactions much farther up the periodic 
table than one could have hoped before. 

The interest and cooperation of Professor 
Oppenheimer and M. Phillips in this work has 
been of great value to us, and it is with real 
pleasure that we acknowledge it. We also express 
our appreciation of the helpful financial support 
given by the Research Corporation, the Chemical 
Foundation, and the Josiah Macy Junior 
Foundation. 


APPENDIX ON THE BINDING ENERGY OF THE DEUTERON 


The binding energy of the deuteron can in principle be 
calculated precisely from the disintegration energies of 
various reactions. Actually the results are subject to some 
uncertainty, because of possible errors in the values given 
for reaction energies and nuclear masses. In particular, 
there are intrinsic difficulties in measuring the energies of 
reactions involving neutrons. As an example of this sort of 
calculation we have chosen a set of reactions which can be 
combined so as to make most of the masses cancel out, and 
in which the direction of the uncertainty in the measure- 
ment of the energies is known, so that definite upper and 
lower limits can be set. These are: 


Li? +H'—2He'+17.06+0.06 MV, (8)8 
Li®+H*+2He!+22.06+0.07 MV, (9)!3 
H?+H2+H!+H'+43.97+0.02 MV, (10)! 
Li? +H?+2He!+n+4 14.60+0.25 MV, (11)8 
H?+H2*+He®+n+2.8+0.2 MV, (12) 
Li®+H'—He+He'+3.540.3 MV, (13) 
Li®+n—H'+Het+4.6+0.3 MV. (14)"8 


The energies of the particles emitted in reactions (8) to 
(11) have been measured with the greatest precision, and the 
reaction energies deduced from these are known with cor- 
responding precision, except in the case of reaction (11). 
This isa three-body disintegration, and it was assumed that 
the most energetic a-particle comes off with the greatest 
energy allowed by the momentum laws. If this ideal limit is 
not actually attained, the reaction energy must be greater, 
and therefore the value given is a lower limit. 

The energy given for reaction (14) was calculated from 
the observed H®* range, using Mano’s" range-energy rela- 
tion. In a reaction of this type, involving the capture of a 


18 Oliphant, Kempton and Rutherford, Proc. Roy. Soc. A149, 406 
(1935). 

14 Dee and Gilbert, Proc. Roy. Soc. A149, 200 (1935) 

15 Oliphant, Kinsey and Rutherford, Proc. Roy. Soc. A141, 722 (1933). 

1% Chadwick and Goldhaber, Nature 135, 65 (1935); Taylor and 
Goldhaber, Nature 135, 341 (1935). 





neutron, it is known, particularly from Kurie’s'’ work, that 
some energy may be radiated during the capture; therefore 
the reaction energy given is a lower limit. 

Since three of these reactions involve a neutron, three 
independent values of the deuteron binding energy can be 
obtained from them. The combinations giving the greatest 
cancelation of masses are: 


(8) —(11), giving H'+n—H?=2.46+0.3, 
(13) —(9) —(12), giving H'+n—H? =2H?— He*— 21.36+0.4, 
(14) — (9) — (10), giving H'+n—H? =2H?— He*— 21.43+0.4. 


H'+n—H? is the binding energy sought. The only masses 
occurring on the other side of the equations are in the 
combination 2H?—He*, which is equivalent to 23.8+0.2 
MV using the new values,'* and to 23.4+0.2 MV using the 
old ones. Putting the first of tnese into the equations, the 
three values obtained for the binding energy are: 


(8)—(11), 7=2.46+0.3, 
(13) —(9) —(12), 7=2.44+0.45, 
(14) -—(9) —(10), 7=2.37+0.45. 


Recalling that the energies given for the reactions (11) and 
(14) are lower limits, one sees that the first value of J is an 
upper limit, and the third is a lower limit. 

These values can be compared with that obtained by 
Chadwick and Goldhaber'’ from the process: 


H?+hv(2.6 MV)—H!+n+0.5 MV, 


giving J=2.1 MV. The possible error in the measurement 
of the energy given to the ejected proton is uncertain, but 
may be some tenths of a million volts. The observed cross 
section for this reaction is smaller than that calculated from 
the theory of Bethe and Peierls,?° using 7 =2.1; the value 
of J required by the theory to fit the observed cross sec- 
ion is 2.4. 

17 Kurie, Phys. Rev. 47, 97 (1935). 

's Aston, Nature 135, 541 (1935). 


1® Chadwick and Goldhaber, Nature 134, 237 (1934). 
20 Bethe and Peierls, Proc. Roy. Soc. A148, 146 (1935). 
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Note on the Transmutation Function for Deuterons 


J. R. OprpENHEIMER AND M. PuILLips, University of California, Berkeley 
(Received July 1, 1935) 


We consider the effect of the finite size and ready polarizability of the deuteron on the 
probability of transmutations involving the capture of the neutron. These have as a conse- 
quence that the Coulomb repulsion of the nucleus is less effective than for alpha-particles or 
protons, and that the corresponding transmutation functions increase less rapidly with deuteron 
energy. We treat the collision by the adiabatic approximation and obtain quantitative results 
for this energy dependence which are in good agreement with experiment. 





ANY elements can be rendered radioactive 

by deuteron bombardment, the reaction 

involving the capture of a neutron and the libera- 
tion of a proton: 


oA” +-H2—,A"*! +H" (1) 


Four reactions of this type have been studied in 
detail by Lawrence, McMillan and Thornton.! 
The transmutation functions which are found in- 
crease smoothly with deuteron energy, but the 
increase is far less rapid than we should expect on 
the basis of the familiar considerations of 
Gamow? on the penetration of charged particles 
through the potential barrier of the nuclear 
Coulomb field. To account even roughly for the 
observations on this basis, we are forced to as- 
sume that the Coulomb potential breaks down at 
very large distances (~ 1.5 X 10~” cm for copper). 
The transmutation function is thus anomalous 
when compared to that for protons and alpha- 
particles, and it is natural to associate this 
anomaly with the structure of the deuteron, 
particularly its low stability. We wish to show 
in this paper that when this is taken into account, 
it does in fact provide a satisfactory explanation 
of the experiments. 

For neutron capture to be possible the neutron 
must have an appreciable probability of coming 
within the range of the nuclear forces. But this 
condition can be satisfied even when the center 
of mass of the deuteron lies beyond the range of 
these forces. It is this possibility which leads to 
an explanation of the fact that even with such a 


1 Lawrence, McMillan and Thornton, preceding article. 
We are greatly indebted to the authors for the opportunity 
of seeing their experimental results, and for many helpful 
discussions. 

2Gamow, Atomic Nuclei and Radioactivity (Oxford Uni- 
versity Press, 1931). 


highly charged element as copper nuclear trans- 
mutations can occur for deuteron energies of the 
order of 2 MV. 

The quantitative treatment of the correspond- 
ing collision problem is complicated, not only by 
our ignorance of the detailed forces involved, but 
by the complete inapplicability of the Born ap- 
proximation. For the velocity of the deuteron is 
not large compared to the internal velocity of 
proton and neutron; the effective time of collision 
of the deuteron is long compared to its period. 
We have thus to use the adiabatic approximation: 
the relative motion of proton and neutron is ap- 
proximately given by the solution of the wave 
equation when the distance of the center of mass 
of the deuteron from the nucleus is held fixed; 
the center of mass moves in an effective potential 
which is the energy &(X) of the relative motion; 
and the perturbation which induces the inelastic 
impacts of the transmutation is the term in the 
kinetic energy neglected in this zeroth approxi- 
mation. In fact the cross section for the transmu- 
tation is then given, with a proper normalization 
of the wave functions, by 


o=1/h?| Sfdpdny;(p, n) 
x [h?/4. Ax +W-I-8(X) Wilx, X) 2. (2) 


Here p and » are the coordinates of proton and 
neutron, x= p—vn their relative coordinates, Y 
the coordinates of the center of mass of deuteron 
measured from the nucleus as origin. Further y; 
is the approximate adiabatic wave function for 
the initial state (normalized to unit flux); ¥, the 
wave function for the final state, in which the 
neutron is captured and the proton is flying off 
with a considerable kinetic energy; and W, J, 2.\/ 
are kinetic energy, binding energy and mass of 
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the deuteron. We wish to evaluate the de- 
pendence of this cross section on the energy of 
the deuteron, insofar as this can be done without 
a detailed knowledge of the structure of the 
nucleus and the nature of the nuclear forces. 
The reactions (1) are highly exothermic, and 


we may therefore neglect the dependence of the 


final wave function for the proton on the deuteron 
energy; except in the immediate neighborhood of 
the nucleus we can take for the final wave func- 
tion of the proton a plane wave. Further, the 
final wave function for the neutron will vanish 
except in the immediate neighborhood of the 
nucleus. Since the neutron is far more stably 
bound in the nucleus than in the deuteron, it is 
reasonable to neglect the finite extension of the 
corresponding wave function. The effect of this 
neglect will be to make the value of the binding 
energy of the deuteron which best fits the ex- 
perimental curves somewhat lower than the true 
value, but the error should. be small, especially 
for an atomic number as high as that of copper. 
We shall take then 


yr~e' hem E)'y§(n), ” 


except for very small values of p. Here E is the 
energy imparted to the proton by the dis- 
integration. 

The essential energy dependence of ¢ is deter- 
mined by the form of y;. The adiabatic ap- 
proximation gives 


Yi=u(x, X)e(X), (4) 


where u(x, X) is the wave function for the rela- 
tive coordinates x, when the center of mass X is 
fixed, and y(X) is the wave function for the mo- 








f(a) =cos (a)'—(a(1—a))}, 


1+<2 
0(s) =—— cot“! (z)?—1, 
(z)! 
1+: 
=(2(1—2))?—14+ (tan ( 
(zs)! 


2s \} ; , 
a ) —tan “@)!), 
l-—s 


— 2g } 
—- ) —tanh '(=2)!). <0 
l-—z 


3 Bethe and Peierls, Proc. Rov. Soc. A148, 146 (1935). 


tion of the center of mass. If V» be the potential 
between neutron and proton, and Vy the nuclear 
potential, then these functions are determined by 
the wave equations 


(h?/MA,+8(X) — Vo(x) — Vu(x, X)}u=0, (5) 
(h?/4AMAx+W—I-8&(X)}o=0. (6) 


According to the familiar arguments of Wigner, 
Vo is given by an extremely narrow and deep 
trough, and is equivalent, as Bethe and Peierls® 
have observed, to the boundary condition for 


0 In u/dx) p90 = —(MI)5/h. (7) 


For Vy we shall take the Coulomb repulsion of 
the proton, and neglect the specifically nuclear 
forces in the immediate neighborhood of the 
nucleus; for these the adiabatic approximation 
can hardly be valid, and they could have a sen- 
sible effect on the transmutation function only if 
resonance occurred; the experimental curves 
afford no evidence for this. 

According to (3) we need consider u only for 
n=0, x= 2X. If we now determine u, & and ¢ by 
the Wentzel-Kramers-Brillouin method, we find, 
with sufficient approximation, 


2X(WI)! Ze ’ 
u(2X, X)=a(X) exp —_—___— “(= = ) | 


m \xr1 J) 
&(X) = —I+Ze?/X (8) 
ig ag, 
v | he Ze /} 


Here Z is the nuclear charge, v the deuteron 
velocity, W= Mv? the deuteron energy, and the 
functions f and @ are defined by 


s>1 


O<s<1 (9) 








502 (a 2 


Further a(Y), b(W, X) are slowly varying func- 
tions. Far from the turning points VY,=Ze® 2/, 
N.= Ze? /W, they are given by 

a(X)=(1-—Ze?/2NX) +; b(W, X) = (Ze /X —-W)-*. 
Since, as we shall see, the values of XY which 
contribute essentially to the matrix integrals for ¢ 
vary little with W, the contribution of a and 6 
to the transmutation function will be neglected. 

We can thus write 


Wi(2X,X)~1/v 


2Ze(M\) sW IX 
exp | —_-— ( ) r( _* ) , (10) 
a XJ Ize} | 
where 
F(q, v) = y@(1/v—1) +q-4f(qy). (11) 


Regarded as a function of Y, F has a rather 
sharp minimum. The position of this minimum 
varies little with deuteron energy in the range 
1.5 MV to 3 MY, and lies roughly at Ze?/2/. 
The breadth of the corresponding maximum of 
¥i(2X, X) decreases slowly with W, but since 
this breadth is of the same order as the wave- 
length #/(2.WE)* of the proton wave function, 
this variation will affect the magnitude of the 
matrix integral very little. The effect of the varia- 
tion of a and } and of the logarithmic derivatives 
of « and ¢ which appear in the integrals for o 
could in principle be taken into account by ob- 
taining accurate solutions of the wave equations 
(5), (6); but without a detailed knowledge of the 
energy of the proton £, and the form of the final 
neutron wave function, such refinements would 
be illusory. Throughout the range //2<W<2/, 
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the variation of o with W is given essentially by 


1 4Ze My? W 
exp ) — ( ) r( ) } (12) 
o? Ll aNg / 


where F(W/J) is the minimum value of F( WJ, 
IX /Ze*). A plot of #(W/7/) is given in Fig. 1. 

The form of the transmutation functions given 
by (12) still depends upon the value of 7. This is 
known to be roughly 2 MV, and this value gives 
satisfactory agreement with the experimental 
curves of Lawrence, McMillan and Thornton. In 
Fig. 2 of their paper Eq. (12) is plotted for three 
values of J (1.5 MV, 2 MV, 2.4 MV) for the 
aluminum reaction. It is seen that the agreement 
with the experimental values is best for /~2 
MV. The approximations in the theoretical treat- 
ment would tend to favor rather too low a value 
of the binding energy of the deuteron. 
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The Radioactivity of Solids Determined by Alpha-Ray Counting 


GLApys D. Finney, Physics Department, University of California at Berkeley 
AND 
RosbLEY D. Evans, Physics Department, Massachusetts Institute of Technology 


(Received June 11, 1935) 


The complete theory of alpha-ray counting for thin and thick radioactive sources is devel- 
oped. Tables of the numerical constants required for the interpretation of practical measure- 
ments are given. Apparatus and experimental verification of the theory are described. Several 
new methods for determining the thorium content of geological samples are described. Analyses 
of a suite of igneous rocks show Th contents between 3.9 and 16.5 10~°g Th per g, with an 
average Th/U ratio of 7. Therefore thorium is more important in age and geothermal problems 





than uranium. 





I. INTRODUCTION 


ov refinements in detection technique 
permit ionization measurements! on the 
alpha-rays emitted by all ordinary materials 
because of their natural contamination by 
traces of radioactive substances. Accurate alpha- 
ray measurements permit the determination of 
the thorium content of geological specimens, 
hence lead toward satisfactory solutions of such 
basic problems as the age of igneous rocks, the 
age of the stellar system as determined from 
meteorite analyses, and the excessive, radio- 
activity-derived internal heat of the earth. 

Feebly radioactive sources such as igneous 
rocks emit the order of one alpha-ray per cm? per 
hour. Any measurements are therefore subject 
to statistical fluctuations due to the random time 
of emission of these alpha-rays. Ionization 
measurements are further impaired by ioniza- 
tion fluctuations due to the variation in residual 
track length of the emitted alpha-rays, which 
must be compounded with the effects of random 
rate of emission in accord with the statistical 
theory of independently random processes.” 
lonization measurements must also deal with 
the unavoidable background due to cosmic-ray 
and local gamma-ray ionization in the detecting 
chamber. 

These difficulties are circumvented by count- 
ing the individual alpha-rays with a detector 
which is unresponsive to cosmic or local gamma- 
radiation. Precision is then limited only by the 
statistical fluctuations in the rate of emission of 


! Evans, Phys. Rev. 45, 29, 38 (1934). 
2 Evans and Neher, Phys. Rev. 45, 144 (1934). 


alpha-rays, and by the accuracy of the compu- 
tations dealing with the finite thickness of the 
source. The statistical problem concerning sev- 
eral alpha-emitters of widely different half- 
periods but in radioactive equilibrium has been 
treated by Adams.’ The equations describing 
thickness corrections for sources containing any 
number of radioactive elements are here de- 
veloped and confirmed by direct experiments. 


Il. THEORY OF PARTICLE COUNTING 


Consider the section of a radioactive solid 
bounded by the planes AA’ and BB’ in Fig. 1, 
and surmounted by absorbers of thickness a; and 
dz cm. The source emits V alpha-rays per cm* 
per sec., of mean range R air-cm (i.e., measured 
in air at 0°C, 760 mm Hg). The full alpha-ray 
range in the source would be u«R cm, and in the 
absorbers 4, R, and we, respectively. The source 
is wr cm (or 7 air-cm) thick and some of the 
alpha-rays which it emits are observable in the 
detecting chamber located above the plane DD’. 














D t D’ 
a 
e P2 i J e 
= a she B’ 
Mor vl, 
dy 
dA 
ral A’ 





Fic. 1. Geometry for Eqs. (1) et seg. Errata: Subscripts 
for a; and a, should be interchanged. 


3 N. 1. Adams, Jr., Phys. Rev. 44, 651 (1933). 
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A. Thin sources 

There are N dy dA alpha-rays per sec. emitted 
in all directions from the volume element dydA 
located y cm below BB’. Alpha-rays which have 
equivalent air ranges of 7 or less below the inter- 
face DD’, hence residual ranges equal to or 
greater than (R—r) in the detecting chamber 
above DD’, will be emitted within a cone of half- 
angle 6, or fractional solid angle (1—cos 6)/2, 
where r cos 0=y/u+ai/uitde/u2. Hence the 
number of such rays from the elementary volume 
is N dy dA(1—cos 6)/2, and n,’ the total number 
of such rays per cm? of surface of a source of 
large area (where edge effects are negligible) is 
obtained by integrating over y between the 
limits: y=0, and y=yur 


N ¢* 
n,=— f (1—cos @)dy 
2 9 


N *" y/et+a1/uitde/ pe 
=— [ [1 ——_—_— a (1) 
2 Yo 


r 


Nurf2(r—a)—r 
= =| —_— | where 
4 r 


@=41/Mitde/pe-+:aircm. (2) 


Now if the detecting instrument above DD’ 
responds only to alpha-rays having a residual 
range equal to or greater than p air-cm, then the 
number of counted alpha-rays per cm? per sec. 
is given by Eq. (2) when r= R—p, that is: 


n' = (Nur/4)-[2(R-—p—a)—7]/(R-p), (3) 


which is the counting rate for thin sources 
surmounted by absorbers.* 

If the thin source has an area A then the total 
alpha-ray counting rate is m,/A. When d is the 
density of the source, its total mass m=yrAd, 
and n,/"urAd=n,'A, where n,” is the total 
alpha-ray counting rate per g of source. There- 
fore: 





n'=n,'/urd 
= (N’’/4)[2(R—p—a) 
—(m/uAd)]/(R—-p), (A) 


where N/d=N” alpha-rays emitted per sec. per 
g of source. When an infinitely thin source is 


4 Evans, Phys. Rev. 45, 29 (1934). 
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used without an absorber r=m=a=0, and Eq. 
(4) reduces to n’’=N’/2. Only in this limit- 
ing case is the counting rate independent of R. 
Under these conditions the alpha-ray counting 
rate is a direct measure of the rate of helium 
production, hence could be directly combined 
with a helium analysis on a bulk rock specimen 
to determine its age. 

N and N” may be expressed in terms of the 
fractional concentration of radioactive material 
Q g per g, its decay constant \ sec.~', its atomic 
weight W, and the density d of the source; thus: 


N= 6.06 X 10°QAd/W, (5) 
N" = N/d= 6.06 X 10"QX/W. (6) 
When the source thickness increases to 


t= R—p-—a, any further increase in thickness 
will not increase the counting rate because even 
the alpha-rays directed normal to the plane DD’ 
will emerge with residual ranges less than p 
and hence cannot be counted. The thickness 
7=R-—p-—a is the transition between thin 
sources, whose theory has been given above, and 
thick sources, which are treated later. The varia- 
tion in counting rate with source thickness may 
be obtained by employing Eq. (3) to derive the 
ratio of n,’ the counting rate above a thin source 
to mn, the rate above a thick source, both having 
the same absorber a superposed. 

The counting rate m, for thick sources (i.e., 
72> R—p-—a) follows from Eq. (3) for a source of 
critical thickness r= R~—p—a. Substituting, we 
find: 

na= (Nu/4)(R—p—a)?/(R—p). (7) 


Dividing Eq. (3) by Eq. (7) we obtain the varia- 
tion in counting rate with source thickness. 


na’ /Ma=1[2(R—p—a) —7]/(R—p—a)* 
=1-—(1-—17/(R—p—a)} (8) 
=1-(1-FY, (9) 


where F, is the fractional thickness r/(R—p—a) 
of the source, and 0< F,<1, as O& (n,'/ma) <1. 
In Fig. 2, the ordinates at the right are n,'/ma, 
when the abscissae are F,. 


B. Thick sources 


Sources are called thick when 72 R—p—a, in 
which case the counting rate is independent of 
source thickness and is given by Eq. (7). 
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Fic. 2. Left-hand ordinates show fractional counting 


rate from thick sources surmounted by absorbers of frac- 
tional thickness F,z=a,(R—p); see Eq. (11). Right-hand 
ordinates show fractional counting rate from thin sources 
having a fractional thickness F; =7/(R—p—a); see Eq. (9). 


The variation with absorber thickness of the 
counting rate due to a thick source is best 
computed in terms of the counting rate n for 
zero absorption. Thus from Eq. (7): 


n,/n=(R—p—a)?/(R—p)? 
=[1—a/(R—p) F (10) 
ares (11) 


where F,, is the fractional thickness a/(R—p) of 
the absorber, and 0S F, £1, as OS (m,/n) <1. In 
Fig. 2, the ordinates at the left are n./n when 
the abscissae represent F,. 

Combining Eq. (8) for a=0, with Eq. (10), 
we find: 


n'/n+n,/n=1—[1—7/(R—p) ? 
+[1—a/(R—p)}. (12) 


Hence if the thin source has a thickness 7, and no 
absorber, its counting rate n’ may be immediately 
obtained from curves for the counting rate n, 
due to a thick source with an absorber of thick- 
ness a equal to +r. Or, analytically, Eq. (12) 
reduces to: 


n'/n+n./n=1 when r=a. (13) 


Eq. (13) holds for each particular alpha-ray 
range R, hence also for the integrated effects of 
several alpha-ray emitters in the same source. 
It is the basic correspondence theorem correlat- 
ing the theories of thin and thick sources, and its 
use halves the numerical calculations required 
for applications to sources emitting alpha-rays 


from many or all members of the U, Th and Ac 
series. It can, of course, be verified by inspection 
of the geometry involved. 

It has been shown previously‘ that the de- 
creased counting rate due to alpha-particle 
straggling in the solid is almost completely offset 
by the gain in counting rate occasioned by the 
symmetrical character of the error function 
describing straggling. The fractional loss is 
An/n = }(p3)?, where pr is Brigg’s range straggling 
coefficient. The fractional loss is thus always less 
than 0.3 percent. 


III. TotraL CounTING RATE FOR SOURCES 
CONTAINING MANY RADIOACTIVE 
ELEMENTS 


Measurements on chemically prepared ma- 
terials containing only one alpha-ray emitting 
element are easily interpreted by Eqs. (3) and 
(7). A more complicated case is met in measure- 
ments of the total alpha-ray emission from 
materials containing all members of the U, Th 
and Ac radioactive series, such as rocks. We 
shall consider the case of thick sources; values 
for thin sources may be readily obtained by 
use of the correspondence theorem of Eq. (13). 

We shall regard the specimen as sufficiently 
old to have reached radioactive equilibrium. 
Then the number of alpha-rays produced per 
cm’ per sec. is essentially the same for each 
member of a radioactive decay series. Substi- 
tuting® the decay constant and atomic weight of 
U and of Th in Eq. (5) we find: 


6.06 X 10° X 4.77 X10-"S 


[hig tiicmnctatscsnententns t 
238.14 
= 1.2110! U-d. (14) 
6.06 X 107° & 1.7 & 107'8 
Weg teen a 
232.12 
= 0.4410! Th-d. (15) 
Nac=0.04Ny = 0.05 X 104 U-d. (16) 


®* Kovarik and Adams, Phys. Rev. 40, 718 (1932) give 
\=4.865 X 1075 sec.~! for Ul. When corrected for AcU, 
this value becomes 4.77 X 107'* sec.~'. Geiger and Ruther- 
ford, Phil. Mag. 20, 691 (1910), and H. Fesefeldt, Zeits. f. 
Physik 86, 605 (1933), find \=1.7X107'*§ sec. for Th. 
A. v. Grosse, Phys. Rev. 42, 565 (1932) reviews the work 
on Ac and recommends 0.04 for the Ac/U activity ratio, 
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Since each alpha-ray body in the U series 
emits 1.2110'U-d alphas per sec., the total 
counting rate for a thick source containing only 
the U 


(1.21 108 U-u-d/4)=(R—p—a)*/(R—p) 


series would be: 


where the summation extends over all the range 
values R which are members of the uranium 
series. 

In a natural source, such as a granitic rock, 
the Th and Ac series will be present also, hence 
similar expressions must be written for them. 
The total observed counting rate 7, will then be 
due to terms from the U plus Ac series and from 
the Th series, or 


n,=v,U +w,Th (17) 


where U and Th are the uranium and thorium 
concentrations in g per g of source. Writing 


I1=S(R—p—a)?/(R—p) (18) 
we have: 
Ve= (1.21 XK 10%ud 4)//¢ 
+ (0.05 X 104ud 4)17 4. 
=[3.03//¢ +0.12/7,. J10%ud (19) 
Wa=[1.1177 7, 10% ud (20) 


where the summations are made over the mem- 
bers of the U, Th and Ac series as indicated, and 
IT for the CC’C” bodies considers their respective 
branching ratios. 

The computation of /7y, /7,. and J/7,, may be 
considerably simplified by the following devices. 
Let Aj be a numerical parameter defined as the 
ratio of the number of alpha-rays per second 
from any element, such as Ra, to the number 
from its parent, in this case U I. Then, without 
practical error, Aj=1 for all alpha-emitting 
elements, except those in the CC’C”’ branches, 
when radioactive equilibrium is present. For 
Ra C, Aj=0 (accurately, 0.0004); for Ra C’, 
Aj=1 (0.9996); for Th C, Aj=0.35; for Th C’, 
Aj=0.65; for Ac C, Aj=1 (0.997); and for Ac C’, 
Aj=0 (0.003). To account analytically for 
branching we therefore introduce the param- 
eter Aj into Eq. (18), obtaining: 


II=SAj(R—p—a)?/(R—p) (21) 
= YAj(R—p)[1—a/(R—p)}?? 
= SR-Aj+aAj/(R—p)—(p+2a)=Aj. (22) 
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TABLE I. Mean range R tn air at 0°C and 760 mm Ig; 
and Aj, j, Pj, Qj (see Eg. (23)) for all alpha-ray emitters 
in the U, Ac and Th series. 


URANIUM SERIES 


R Aj Jj P; QO;(p =0.5) 
Ra C’ 6.51 0.99982 1 651 0.1665 
Ra A 4.38 1 2 10.89 0.425 
Ra C (av) 3.80 0.00022 2 10.89 0.425 
Rn 3.80! 1 3 14.69 0.728 
Ra F 3.61! 1 4 18.30 1.050 
Ra 3.165 1 5 21.46 1.426 
ULI 3.058 1 6 24.51 1.818 
lo 2.978 1 7 27.48 2.223 
Ul 2.515 1 8 29.99 2.721 
ACTINIUM SERIES 
Ae C’ 6.17! 0.0035 0 0 0 
Ac A 6.08! 1 1 6.08 0.179 
An (av) §.291 1 2 11.37 0.388 
Ac C (av) 5.04) 0.9973 3 16.41 0.6008 
Rd Ac 4.3063 1 4 20.77 0.867 
Ac X 4.07 1 5 24.84 1.147 
Pa 3.423 1 6 28.26 1.489 
Ac U 3.04 1 7 31.26 1.889 
THORIUM SERIES 

Thc’ 8.991 0.658 0.65 5.26 0.0857 
Tha 5.31! 1 1.65 10.57 0.294 
Tn 4.71! 1 2.05 15.28 0.532 
Thc 4.45! 0.354 3 16.84 0.620 
Thx 4.005 1 + 20.90 0.901 
Rd Th 3.743 1 5 24.04 1.210 
Th 2.065 1 6 27.30 1.673 


! Lewis and Wynn-Williams, Proc. Roy. Soc. A136, 349 (1932). 
? Rutherford, Wynn-Williams, Lewis and Bowden, Proc. Roy. Soc. 
A139, 617 (1933). 

> Rutherford, Chadwick and Ellis, Radiations from Radioactive Sub- 
stances (1930), p. 26. The extrapolated ionization ranges R; there 
given are converted to mean ranges using the relation® Rj —R =0.042 
(Ry —- 1.51)? cm. 

4 Evans, Phys. Rev. 45, 29 (1934), estimated from the half-period, 
A. v. Grosse, Phys. Rev. 42, 565 (1932), and the Geiger-Nuttall relation 
using Geiger's, Zeits. f. Physik 8, 45 (1921), data for the Ac series. 

5 Kurie and Knopf, Phys. Rev. 43, 311 (1933), give the extrapolated 
ionization range R; in air at 0°, 760 mm as 2.71+0.03 cm. The mean 
range is estimated by deducting 0.042 (Rj —1.51)4 cm [Briggs, Proc. 
Roy. Soc. All14, 313 (1927); Kurie, Phys. Rev. 41, 701 (1932); Evans, 
Phys. Rev. 45, 29 (1934)]. Using Lewis and Wynn-Williams (refer- 
ence 1, Table I) relationship, Re—R =0.8p —0.06, where p is the strag- 
gling coefficient (p =psx used by Briggs) one obtains R =2.67 for Th. 
G. H. Henderson and Nickerson, Phys. Rev. 36, 1344 (1930) obtained 
a lower value Rj =2.59+0.05 cm at 15°C, 760 mm, while Geiger, Zeits. 
f. Physik 8, 45 (1922) found Rj =2.90 cm at 15°C, 760 mm. 


Now we define P;=R- Aj, Q;=Aj/(R—p), and 
j=XAj; then Eq. (22) becomes: 


H= P;+a°Q;—j(p+ 2a). (23) 


In evaluating P;, Q; and 7 the summation ex- 
tends only over those alpha-rays having R2a+p, 
since only these contribute countable rays. 
Table I gives the values of R, Aj, j, P;, and Q; for 
all alpha-rays in the U, Ac and Th series. 

As an example, consider a thick radioactive 
source containing only the Th series. Let the 
minimum detectable track p=0.5 air-cm and the 
absorber a=4 air-cm. Then only alpha-rays hav- 
ing R2 4.5 air-cm may be counted. From Table I 
this is seen to include Th C’, Th A and Tn, but 
neither Th C nor shorter members. Reading 
across the row corresponding to the shortest 
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TABLE II. Dependence on absorber thickness, a air-cm, of 
the relative counting rate vq/v for the U plus Ac series in 
equilibrium, wa/w for the Th series, and vq/wWa the ratio of 
U’ plus Ac series to Th series counts, All values correspond 
to p=0.5 air-cm. Values of vq and wae (see Eqs. (19) and 
(20)) are also given. 








a ta ‘10° ud wa 108 ud ta/ta ta/? Wa /w 
0 82.1 2 7 0 3.04 1.000 1.000 
1 40.4 15.52 2.60 0.493 0.575 
2 15.62 7.75 2.02 0.1905 0.287 
3 5.63 3.36 1.68 0.0687 0.124 
4 2.12 1.40 1.51 0.0258 0.052 
5 0.529 0.633 0.84 0.00645 0.0236 
6 0 0.237 0 0 0.0089 





countable emitter, Tn, we find j7=2.65, P; 
=15.28, Q;=0.532. Hence by Eq. (23), Hn 
= 15.28 + (4)? X 0.532 — 2.65(0.5 + 2 K 4) = 1.26. 
Finally, from Eqs. (17) and (20), the number of 
alpha-rays counted per sec. per cm* of source 
surface is 2,=w,Th=1.11 X 10°//7,ud-Th = 1.40 
K10%ud-Th, when a=4, p=0.5 air-cm. This 
value will be found in Table II. 

Writing subscripts 1 and 2 to characterize two 
values of absorber thickness, a; and a2, we may 
write the equations for absorption measure- 
ments which might permit separation of the U 
and Th effects, and hence provide a direct 
measure of either or both of these elements. 

my=vuU+wiTh, ne=veU+weTh. (24) 


Combining these equations we obtain: 


Th = (11 /ws) (2/01 — 02/01) / 
(We/w1—V2/t1), (25) 


—n»/n,)/ 
(w '9/W,—Ve/V1), (26) 


Th/U = (2: /w1) (m2/1 — 02/01) / 


(w 2/W;—N2/ny), (27) 


n,/Ne= (w1/We)(21/@1+Th/U)/ 
(v2/we+Th/U). (28) 


U=(n, i 2) (We/ u 


Here ; and m2 are values which can be observed 
experimentally, while 7, v2, @:, @2 are computed 
theoretically, or obtained experimentally by 
internal standardization, as described in Sec- 
tion IV. 

Numerical examination of Eqs. (25), (26), 
(27), (28) quickly shows that an impossibly high 
experimental precision is required if U and Th 
concentrations are to be obtained from experi- 
ments in which the two absorbers differ in thick- 
ness by only about an air-cm. If, however, a; 


approaches zero, while a, is made the order of 5 
or 6 air-cm, then separation of U and Th be- 
comes possible even for quite feebly radioactive 
sources. Due to the low background counting 
rate attainable in counting experiments these 
measurements become possible, whereas they 
are forbidden by the relatively high gamma- 
and cosmic-ray background if ionization methods 
are employed. 

It will be noticed by comparison with Eqs. 
(19) and (20) that Eqs. (25) and (26) express the 
Th and U content in terms of (1/ud) since the 
product ud occurs once in each term involving 
v or w. Experimentally, the effective stopping 
power (1/u), is not readily determinable for 
rocks, particularly if they are coarse grained. 
The product zd can be approximated, however, 
from the Bragg-Kleeman rule which gives 


ud,=d.W,'/W.’, (29) 


where d,= density of air at 0°C, 760 mm; 
d,=density of the mineral grain; W,!=sum of 
the atomic percent times the square root of the 
atomic weight of each of the chemical consti- 
tuents of the mineral grain; W,)=the similarly 
computed mean square root of the atomic 
weight for air. Direct measurements of ud for 
mica have been made by Briggs® and by G. H. 
Henderson*® who obtained 1.547 X 10-*g/cm? and 
1.659 X10-%g/cm?, respectively. Computations 
of ud for mica have been given by G. H. Hender- 
son® and by Kerr-Lawson,’ who found 1.629 
X 10-%g/em? and 1.639 x 10-*g/cm?, respectively. 
Since W,! is about 4.79 for mica (Kerr-Lawson) 
we may say that within an uncertainty of 
about +3 percent 


ud, =3.4X10-W,}. (30) 


This value is also obtained when the known 
values of d, and W,' are substituted in Eq. (29). 

For practical ge ae of Eqs. (24) to 
(28) we need tables of v2/2:, we/w,, and the 

values of v and w when a =0. These are obtained 
from the data of Table I combined with Eqs. 
(19) and (20), and are given in Table II for the 
case in which p=0.5 air-cm, which is an easily 


6G. H. Briggs, Proc. Roy. Soc. Al14, 341 (1927). G. H. 


Henderson, Proc. Roy. Soc. A145, 563 (1934). D. E. Kerr- 
Lawson, Univ. of Toronto Studies, Geological Series, No. 
24, p. 54 (1927). Each of the values of ud quoted above 
has been corrected to 0°C, 760 mm Hg pressure. 




















508 G. D. FINNEY 
1,0 T T T my T T T 
08 a 
0.6 d 

r TH 
0.4 “”~ d 
Us#ac 
O2- PA] a 
° 1 1 1 ; _ 
0 ' 2 ‘.* 5 6 7 a 


lic. 3. Dependence of the relative counting rate vq v for 
the U plus Ac series in equilibrium, and w./w for the Th 
series, on absorber thickness a air-cm, when p=0.5 air-cm. 


obtained experimental condition. Table II cor- 
responds to Eqs. (24) to (28) when a;=0, a2=a; 
0, = =U, Wi = Wo =W; and V2=Va, We=Wa. Fig. 3 
shows graphically the dependence of v./v and 
uw’, /w on a. 


IV. DETERMINATION OF THE THORIUM CONTENT 
OF GEOLOGICAL SPECIMENS 


Completely satisfactory experimental methods 
are available for determining the radium (hence 
uranium) content of rocks.’ Geological age and 
thermal problems also require accurate knowl- 
edge of the thorium content of rocks. The results 
of Section III lead to several new methods for 
determining the thorium content of rocks® by 
counting the alpha-rays emitted from the rock 
surface. 

Method A rests on absorption measurements 
and the use of Eq. (25). The value of ud for the 
source may be estimated from Eq. (30), or 
obtained experimentally by internal standardiza- 
tion, that is, by adding a known amount of a 
finely powdered Th mineral such as monazite 
sand to a powdered rock specimen and determin- 
ing the change in the observed counting rate. 

Method B also involves absorption measure- 
ments, but employs Eq. (27) to evaluate the 
Th/U ratio. This method is independent of the 
value of ud for the source, but must be combined 
with a separate determination of U by emanation 
measurements. These can be made with high 
precision. The method B is really equivalent to a 


7 Evans, Rev. Sci. Inst. 6, 99 (1935). 
* Evans, Finney, Kip and Mugele, Phys. Rev. 47, 791 
(1935), 
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Fic. 4. Left-hand ordinates: ratio mg," of counting rates 
with and without absorbers; curve m; » for a=5 air-cm, 
curve m for a=6 air-cm. Right-hand ordinates: ratio 
nen; of counting rates with a=6 toa=5 air-cm absorption. 
Abscissae: Th, U ratio in g per g. 


combination of method A with a determination 
of ud by means of Eq. (26) and the independently 
measured uranium content. Fig. 4 shows the 
dependence of the relative counting rate n,/n 
with 5 and 6 air-cm absorbers to the counting 
rate with no absorber, and also n¢/n; the rela- 
tive counting rate with 6 air-cm and with 5 
air-cm absorbers, each as a function of the Th/U 
ratio. It will be noted that the method is best 
suited to Th/U ratio less than about 6; for 
higher Th/U ratios, method C is preferable. 
Method C is especially applicable to very 
weak sources in which absorption measurements 
would reduce the counting rate to an imprac- 
tically low value. It is, therefore, generally ap- 
plicable even when methods A and B become 
objectionally inaccurate due to statistical fluctua- 
tions in the counting rates. It consists of simply 
measuring the total counting rate from a source 
with no absorber, hence by Eq. (24) and Table II: 


n= (82.1U+27.0Th) 10%ud. (31) 


As in method A, ud is evaluated by internal 
standardization using either a thorium or a 
uranium mineral, or ud is estimated from Eq. 
(30). In the latter case, it must be remembered 
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that if a coarsely powdered rock is used as a 
source, the density of the source d in Eq. (31) is 
less than the density d, of the individual grains 
treated in Eq. (30). Since the ratio d/d, is not 
easily measured, the method of internal stand- 
ardization gives a more reliable value of ud. 
A separate measurement of the uranium content 
is then used with Eq. (31) for the determination 
of the thorium content of the specimen. Since 
knowledge of the U content is always required in 
age or thermal problems, this involves no addi- 
tional work and is therefore not objectionable. 

Where ATh is the known change in Th con- 
tent, introduced as an internal standard, and 
An is the corresponding change in the counting 
rate n, we have 


Th=n(ATh/An) — (v,4/w.a)U (32) 
=n(ATh/An) —3.04U(when a=0). (33) 


Method D. For sufficiently strong sources a 
fourth method is available. This consists of a 
single measurement of the counting rate with an 
absorber sufficiently thick (a+p=6.5 air-cm at 
0°C 760 mm) to stop all but the 8.08 air-cm 
alpha-ray from Th C’. This method is really a 
special limiting case of method A, for when 
ve=0, i.e., complete U series absorption, Eq. (25) 
leads to: 


Th=n2/ws=4.22 X10. /ud 
(for a=6, p=0.5), (34) 


where 7, is the counting rate per cm’ per sec., and 
ud is obtained either by internal standardization 
or from Eq. (30). Taking ud =1.63 x 10-*, as for 
mica, we then have Th=2.6n,. If the source is a 
granitic rock having 10-°g Th per g the total 
counting rate through the 6 air-cm absorber will 
then be 10-> 150 X 3600/2.6 = 2.8 alphas per hr. 
per 150 cm’, which is detectable on the apparatus 
described in Section V. For sources containing 
as much as 10 g Th per g the method becomes 
quite accurate, as statistical fluctuations are 
then much less pronounced. 


V. APPARATUS FOR COUNTING ALPHA-RAYS 
FROM LARGE SURFACE AREAS 


The counting chamber (Fig. 5) is a parallel 
plate condenser having a lower collecting elec- 
trode 10 cm in diameter. This electrode is sur- 
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Fic. 5. Construction of the counting chamber, source 
holder, and electrometer-tube housing. Accurately to scale. 


rounded at a distance of 1 mm by a circular 
guard ring 2.5 cm wide. The edges of both are 
beveled to reduce the electrostatic capacity of the 
collector. The upper electrode is 2 cm above 
the lower collecting plate, and is mounted on a 
brass ring supported by Bakelite blocks. In the 
early experiments a grid of parallel, gauge 40, 
German silver wires spaced 5 mm apart was used. 
The alpha-ray source was mounted immediately 
over it with the active face downward and the 
active area defined by a circular opening in a 
brass mask which rested on the upper electrode 
ring and supported the source. 

The source is either a flat surface polished 
onto a hand specimen of a rock; a chemically 
prepared precipitate painted onto an aluminum 
disk by using ethyl alcohol as a diluent; or a 
finely ground rock powder painted on in the same 
way. The source is covered with a sheet of thin 
aluminum or gold leaf (stopping power 0.03 air- 
cm) fastened to the periphery of the source disk 
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with vaseline. This covering provides a uniform 
electrostatic field normal to the surface of the 
source, and prevents minute free grains of the 
source from dropping into the counting chamber. 
Gold leaf will not prevent the escape of Rn and 
Tn from polished rock surfaces, but direct tests 
by several methods show that the powdered 
rock sources do not permit the escape of a 
significant amount of emanation. 

Satisfactory electrical operation can be realized 
without the use of this covering foil, even when 
the wire grid is removed from below it, but we 
have always used the foil in order to provide 
exactly reproducible. collecting fields for all 
specimens. In this case it is best to remove the 
wire grid, making the foil which covers the 
source act as the upper electrode of the collecting 
chamber. This procedure also eliminates the 
absorption of alpha-rays by the wires of the grid 
and by the gas between the source and the grid, 
and avoids the uncertainties of penetration of 
the collecting field beyond the wires of the grid. 
The equations of Section III permit accurate 
correction for the small absorption in the foil. 
For absorption measurements, aluminum foils 
are placed directly between the source and the 
chamber, and nearly in contact with the source. 

In addition, powdered sources similarly pre- 
pared on Al foils may be placed directly on the 
lower collecting electrode with their active sur- 
faces facing upward into the counting chamber. 
It is thus possible to utilize over 150 cm? of 
active source material. 

The lead from the lower collecting electrode 
passes through an amber insulator into an 
evacuated iron tube-housing patterned after 
Bearden’s® apparatus. The amplifying circuit 
follows the methods of DuBridge and Brown," 
and the output galvanometer deflection is con- 
tinuously photographed with a telechron driven 
drum camera nearly identical with that de- 
veloped for string electrometer measurements 
by Evans.? Eastman P.M.C. No. 2 contrast 
bromide emulsion on smooth single-weight 
paper, 5 inches wide, is used for recording. 
The drum is 12 inches in diameter and rotates 
once an hour, hence the travel is 1.6 cm per min. 


9 Bearden, Rev. Sci. Inst. 4, 271 (1933). 
 DuBridge and Brown, Rev. Sci. Inst. 4, 532 (1933). 
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This high film speed permits recording the 
full form of the alpha-ray deflection and gives a 
distinctive hump with a steep initial slope, cor- 
responding to the ballistic period of the gal- 
vanometer; a maximum height corresponding 
to the range of the alpha-ray in the chamber; and 
an exponential return to the undeflected position, 
corresponding to the electrostatic time constant 
of the collecting system. From measurements on 
the shape of the return curve, and knowledge 
of the value of the grid resistance (4 10!'! ohms), 
the electrostatic capacity of the collecting system 
and its leads is found to be 36 cm. 

The sensitivity" is such that a 1 cm alpha-ray 
produces a 1 cm deflection of the galvanometer. 
The maximum noise level fluctuations are about 
2 to 3 mm hence, while shorter alpha-rays can be 
clearly identified, we have taken p as 0.5 air-cm 
(0°C 760 mm of air) for then there can be no 
ambiguity. 

Radioactively dead nitrogen, from tanks 
stored longer than one month, is flowed con- 
tinuously through the collecting chamber and its 
cylindrical brass housing at a rate of about 20 
cm*® per minute, which causes no appreciable 
pressure rise in the chamber. This excludes room 
air from the apparatus and eliminates resultant 
contamination. Great care must be taken in 
preparing sources, and in all operations with the 
apparatus, to avoid traces of contamination. 
Numerous control and background runs ac- 
company every measurement of any source, 
in order to avoid contamination effects. The 
parts of the counting chamber are readily re- 
movable for cleaning, and average background 
counting rates of 5 alphas per hour are realizable, 
while 10 per hour is an easily reproduced routine 
standard. The maximum resolvable counting 
rate is the order of 350 alphas per hour, and 
various classes of rock specimens normally run 
between 20 and 200 alphas per hour per 50 cm? 
of surface area. 

The entire apparatus except the camera is 
housed in a double-walled masonite case for 
protection from temperature variations, and all 
electrical parts are electrostatically shielded. 

't The apparatus has been moved from the University of 
California to the Massachusetts Institute of Technology, 
where, with the assistance of Mr. Howard Tatel, a Western 


Electric 7567-A electrometer tube has been substituted for 
the F P-54 formerly emploved. 
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The galvanometer spot is automatically moved 
across the camera slit so that successive hours’ 
records lie side by side on a single photograph. 
Two methods, each telechron controlled, have 
been successfully used. By means of a 22 contact, 
rotary, magnetic telephone relay, which is fed 
hourly impulses from contacts on the minute and 
second hands of a telechron clock, the plate 
resistance and hence the galvanometer zero may 
be shifted over suddenly once each _ hour. 
Alternatively, the galvanometer lamp, mounted 
on a traveling carriage, may be moved uni- 
formly by means of a screw-thread driven at any 
speed by a gear on the minute hand shaft of a 
telechron motor. The former method was used 
for over 2200 hours of photographic records 
without any failures, but the latter method has 
proved preferable at very high sensitivities, 
because of the low and critical value of the plate 
circuit resistance. 


VI. EXPERIMENTAL RESULTS 


Analyses of suites of rock specimens, protac- 
tinium sources, kolm, etc., will be separately 
reported in other journals. In addition we have 
measured the alpha-ray emission from artificial 
rock sources in order to verify the theoretical 
equations of Section III. The specimens were 
prepared by mixing small amounts of a weak 
uranium ore (carnotite) and a weak thorium ore 
(monazite) with plaster of Paris, and have been 
previously described in connection with ioniza- 
tion measurements on alpha-rays emitted by 
solids."* The uranium series sources contain 
2.5310 g U per g and 2.53X10-* g U per g; 
the thorium series sources 4.25 X10 g Th per g 
and 4.25 10~* ¢ Th per g. 

Taking the value of ud, from Eq. (30), and 
measuring d for the plaster of Paris (=1.32 g 
cm“), it is found that the observed counting 
rate, within the statistical fluctuation of +3 to 


2? Evans, Phys. Rev. 45, 38 (1934). 


5 percent occasioned by the random emission of 
alpha-rays, is equal to that predicted by the 
equations of Section ITT. 

Further calibration tests were made by adding 
known small amounts of the uranium and 
thorium series to powders of igneous rocks, in the 
course of repeated internal standardization tests. 
Here again the experimental observations are in 
close agreement with Eqs. (30), (32), and the 
analog of Eq. (32) for the uranium series. The 
agreement in the case of the thorium standards 
tends to support the Geiger-Rutherford value of 
the Th decay constant which was chosen® in 
Section III as more probable than Kirsch’s 
lower value." 

In general, the thorium analyses of igneous 
rocks have shown values between 3.9 and 
16.5 <10-* g Th per g, and Th/U ratios ranging 
between 3.2 and 10.6, and averaging about 7. 
These results have been confirmed by inde- 
pendent measurements of the U and Th content 
of these specimens by the emanation method" 
and the y-ray method."® They focus attention on 
the great importance of a reliable routine method 
for the determination of the Th content of 
geological materials, because for Th/U ratios 
greater than 4 the Th series produces more 
helium and more thermal energy than the U 
series, and hence is of greater importance 
geophysically. 

The experimental work and a preliminary 
theory for thin absorbers which is equivalent to 
the present general theory was completed at the 
University of California during the National 
Research Council Fellowship tenure (1932-1934) 
of one (E) of the authors. The generous support 
and encouragement given this program by 
Professor Leonard B. Loeb is particularly ap- 
preciated. 


18 (4. Kirsch, Physik. Zeits. 31, 1017 (1930). 

4 Evans, Finney, Kip and Mugele, Phys. Rev. 47, 791 
(1935); Physics (1935), in preparation. 

Evans and Mugele, Phys. Rev. 47, 427 (1935). 
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Double Beta- Disintegration 
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From the Fermi theory of 8-disintegration the probability of simultaneous emission of two 
electrons (and two neutrinos) has been calculated. The result is that this process occurs sutti- 
ciently rarely to allow a half-life of over 10'* years for a nucleus, even if its isobar of atomic 
number different by 2 were more stable by 20 times the electron mass. 





1. INTRODUCTION 


N a table showing the existing atomic nuclei 

it is observed that many groups of isobars 
occur, the term isobar referring to nuclei of the 
same atomic weight but different atomic number. 
It is unreasonable to assume that all isobars have 
exactly the same energy; one of them therefore 
will have the lowest energy, the others are un- 
stable. The question arises why the unstable 
nuclei are in reality metastable, that is, why, in 
geologic time, they have not all been transformed 
into the most stable isobar by consecutive 8-dis- 
integrations. 

The explanation has been given by Heisen- 
berg! and lies in the fact that the energies of 
nuclei of fixed atomic weight, plotted against 
atomic number, do not lie on one smooth curve, 
but, because of the peculiar stability of the 
a-particle are distributed alternately on two 
smooth curves, displaced by an approximately 
constant amount against each other (the mini- 
mum of each curve is therefore at, roughly, the 
same atomic number). For even atomic weight 
the nuclei of even atomic number lie on the lower 
curve, those with odd atomic number on the 
higher one. One §8-disintegration then brings a 
nucleus from a point on the lower curve into 
one of the upper curve, or vice versa. The nuclei 
on the upper curve are all of them unstable. 
But it may happen that a nucleus on the lower 
curve, in the neighborhood of the minimum, 
even though it is not the most stable one, cannot 
emit a single 8-particle, since the resultant isobar, 
whose energy lies on the upper curve, has higher 
energy. This nucleus would then be metastable, 
since it cannot go over into a more stable one by 
consecutive emission of two electrons. This 
explanation is borne out by the fact that almost 


1W. Heisenberg, Zeits. f. Physik 78, 156 (1932). 


only isobars of even difference in atomic number 
occur. 

A metastable isobar can, however, change into 
a more stable one by simultaneous emission of 
two electrons. It is generally assumed that the 
frequency of such a process is very small. In 
this paper the propability of a disintegration 
of that kind has been calculated. 

The only method to attack processes involving 
the emission of electrons from nuclei is that of 
Fermi? which associates with the emission of an 
electron that of a neutrino, a chargeless particle 
of negligible mass. Thereby it is possible to ex- 
plain the continuous $-spectrum and yet to 
have the energy conserved in each individual 
process by adjusting the momentum of the 
neutrino. In this theory the treatment of a 
8-disintegration is very similar to that of the 
emission of light by an excited atom. 

A disintegration with the simultaneous emis- 
sion of two electrons and two neutrinos will then 
be in strong analogy to the Raman effect, or, 
even more closely, to the simultaneous emission 
of two light quanta,* and can be calculated in 
essentially the same manner, namely, from the 
second-order terms in the perturbation theory. 
The process will appear as the simultaneous oc- 
currence of two transitions, each of which does 
not fulfill the law of conservation of energy 
separately. 

The following investigation is a calculation of 
the second-order perturbation, due to the inter- 
action potential introduced by Fermi between 
neutrons, protons, electrons and neutrinos. As 
far as possible the notation used is that of 
Fermi. For a more detailed discussion and justi- 
fication of this mathematical form and the as- 
sumptions involved reference must be made to 
Fermi’s paper. 


2 E. Fermi, Zeits. f. Physik 88, 161 (1934). 
5’ M. Goeppert- Mayer, Ann. d. Physik (V) 9, 273 (1931). 
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2. Tut MATHEMATICAL: APPARATUS 


The nucleus is assumed not to contain any 
electrons and neutrinos but to be built up out of 
neutrons and protons only. Neutron and proton 
are regarded as not essentially different from one 
another, but to represent two different quantum 
states of the heavy particle. The two kinds of 
light particles outside of the nucleus, the elec- 
trons and neutrinos are treated according to the 
method of superquantization. The stationary 
states of the electrons are taken to be those of 
positive energy J/, in the Coulomb field of the 
nucleus, described by four Dirac functions 
¥.=(¥.', ¥.2, v8, v4). Since the neutrinos are not 
affected by the field of the nucleus their eigen- 
functions are represented by plane Dirac waves 
vo=(¢o', ¢o?, ¢o*, ¢e*) with energies K,. The 
Pauli principle is assumed to hold for both 
electrons and neutrinos, so that the number .V, of 
electrons in a stationary state s and the number 
JJ, of neutrinos in a stationary state ¢ can be 
0 or 1 only. 

A quantum state .V of the total system is then 
determined by these numbers .V,, ./, and by the 
quantum state #, with energy W,, of the nucleus; 
n means in this case not only the quantum states 
in the ordinary sense but includes the quantum 
states of all possible isobars. The energy of a 
quantum state .V is given by 


Ey=>NU.4+M.Kot+ Wa. (1) 
8 ¢ 
The interaction energy between the heavy par- 


ticles within the nucleus and the light particles 
without is constructed in such a way that the 





transition from neutron to proton is necessarily 
accompanied by the emission of an electron and 
a neutrino and vice versa. A matrix element of // 
corresponding to the transition of a neutron with 
eigenfunction u, to a proton with eigenfunction 
Um is different from 0 only if at the same time two 
numbers N,, \/, change from 0 to 1, and is then 
given by 
N neeeDgreeO gees 
FL yy = TD ...1,.001¢00 
=(-1 Nites tN ort Miter + Mo ly Bett um, (2) 
with the abreviation 
yig*= — Yl pt? fF? Fl 4 VV ¥3 D4 — 4 D8 (3)* 
Hrum=2S im*u,dr. (4) 


The functions ¥ and ¢ are taken at the place of 
the nucleus, the assumption being made that 
they do not vary considerably over its range. 
For an ‘‘allowed”’ transition fv,,*u,dr, taken 
over the volume of the nucleus, has the order of 
magnitude 1. The value of the proportionality 
constant g can be roughly obtained by equating 
calculated and observed intensity of 6-disinte- 
gration. Fermi determines it to be 


g=4xX10~-°° cm’ ergs. (5) 


3. CALCULATIONS 


The probability of simultaneous emission of 
two electrons and two neutrinos is obtained 
from the second order terms of the perturbation 
theory. In this approximation the amplitude of 
the transition probability is given by 


HyX Hy je@ri (BeBe — | @@eilh(Bu- BR) 1) 





ay= y atetaednals 


® Ex—Ey | 


We want to consider transitions of such a type 
that in the beginning state, V, no electrons or 
neutrinos are present; in the end state, -V/, 
the charge of the nucleus has been increased by 
2, that is two neutrons have been transformed 
into protons, and two electrons have been emitted 
into states s and ¢, two neutrinos into states o 
and +r. The intermediate state K then must be 
such that one neutron has changed into a proton, 
one electron and one neutrino have been emitted. 
This means that the intermediate state k of the 


ieiesieniinienias en (6) 


Bante Ey—Ex 


nucleus is a quantum state of the isobar with 
atomic number larger by 1 than the original 
nucleus. The assumption that ordinary §-dis- 
integration is energetically impossible means 
that W,.,,= W.— W,, > —mce?; the first denomina- 
tor Ex—Ey in (6) will therefore always be 
positive. The process of double 8-disintegration 
is, however, energetically permissible if W,,,, 
=W,,—W,<—2mc*. In this case the second 


* A star denotes the conjugate complex of a quantity. 
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denominator, y—Fy, the difference of energy The summation over the intermediate states 
of the light particles can be made explicitly; 
there are only four intermediate states possible, 
namely, electron s may have been emitted before 


between beginning and end state, will, for those 
processes for which the energy is conserved, 
approach zero, and only then has ay* appreciable 
values. As in the theory of radiation the second — or after electron ¢, and neutrino o before or after 
term in the bracket in (6) can then be neglected neutrino 7+. The power of —1 in the operators 
a,*, b,* introduces a difference in sign in the 


with respect to the first, since /24,— /x« has always 
different cases. The result is: 


a finite value. 


| 1 . : 
y= EM uthall + serbe.* 
LWeetletKe WentI+K, 
1 1 | exp [2 h)( Won tH tM + Ko +K,)t]—1 
6 . ; 


Wait. t+HitKo+K, 


= + (7) 
F mtlHt+K, Wi.tlitKe 

To obtain the transition probability ay* * has to be summed over all possible states s, o, ¢ and r. 

The first step to this is to average the absolute value squared of (7) over all possible directions of the 

plane waves of the neutrino. This is done most easily by assuming that the neutrinos are confined to 


a large space 2. In that case, neglecting the neutrino mass, ¢.’¢o*=(1 42)6,,. We obtain terms of 


the type F . a 
Pidge* 2=(1 42)(a.), Prdee* W*5e.= (1/42) (vd), 


where the abbreviation (¥.¥,) = >-,¥.’*¥:" has been used. 
Furthermore, from the properties of the y’s, the relativistic eigenfunctions of the Coulomb field, 
it follows that, if the spin of the state s is parallel to that of ¢, 


(va) 2=(ba) (bad; 


if the spins of states s and ¢ are antiparallel 


Vv, 2 =0.T 
We therefore obtain 
a. 1 2, Asin? £09k) (Want tM AK +Ky)t] 
ay z= — A(vws) (Wad -— = . 
162° (Wrartls+ih+Ko+K,)? 
where A is the following abbreviation: 
For parallel spin of states s and ¢, 
| 1 1 1 1 2 
A=A,= EM illo} ~~ + — —--- . (8a) 
k Win»tll,+Kao Winrtlht+kK, Wintll +k, Wintlhit+kK,5| 


For antiparallel spin 


A=A2= 


i ' - 
? Wen tH +Ke Wen tl+K,J| 





| 1 1 2 
+| Eadie} noes + ||. (8b) 
| 3 WinthHl+K, Wintli+K-.J| 


Here the sum over & runs over all possible quantum states of the intermediate isobar. 


+ This is true only for ¥ functions with the angular which we are concerned. But all other ~’s are negligible at 
quantum number j=1, 2 in the neighborhood of r=0, with — these distances from the origin. 











9 





4 
— 
41 


DOUBLE BETA-DISINTEGRATION 


The summation over all possible states of the neutrinos is now replaced by an integration over 
the energies K,, A,. The number of neutrino states whose energy lies between K, and K,+dK, is 
given by 


(8rQ/hec*)K dK. 


In the first integration, i.e., over K,, only those energies contribute which lie in the neighborhood of 
the place Ey, — Ey =0, which corresponds to conservation of energy. The integral becomes propor- 
tional to the time ¢ in the usual manner. By neglecting the dependence of A on the energies, the 
probability of emission of one electron into state s, one into the state ¢ becomes 


Pst=(1/t) > |ay |?= (824 hic) A (Webs) Wd S Ke(Wnn tl +i t+Ks)dKe 
= (4r4/15)(A/hic*)(v.v.) Vad (Wamn-H.—I)*. (9) 


Finally the summation over s and ¢ is replaced by a double integration over the energy, that is over 
the continuous spectrum of the electron in the Coulomb field of the nucleus. It is somewhat more 
convenient to express the energies in units of mc*?, namely, to introduce the dimensionless quantities 
mee= Wan >0, meh, =I1,, meh,= . 

According to Fermi the value of y at the outer edge of the nucleus, that is at a distance p=9 X 107" 
cm from the center has to be used; the sum of (¥,¥,) over all states with a fixed direction of spin and 
energy in the range dh, is given by: 


167 mc3 44amcp\ 75 =H 
> v.(p)¥.(p) = ( - ) h(he—1)**8 


x n(3+2S)\2 hk? XA 
h, h, |? 
exp| ~ | r(ts+in — ) dh,, (10) 
(h2—1)!J| (h.?—1)'7 | 


where Z is atomic number, y=Z, 13, S=(1—-y?)!—1. S is for any occurring value of nuclear charge 
negative and small compared to 1; y=0.3, S= —0.05 correspond to Z =41. (10) can therefore be 
approximated by a much simpler form by neglecting S against 1/2 in the exponent of (4,2—1) and S 
against 1 in the ['-function, since 


h, | h, |? h, 
exp (= - ) r(1+in- - ) =2ry—— s 
(h,2—1)*7 | (h?—1)*7 | (h.?—1)? 


The transition probability of the nucleus obtained by summation of (9) over all states and the two 
cases of parallel and antiparallel spin of the electrons, then turns out to be 


4678-2 m''cl9 74armcp\ 48 * 
P= : : (- *) (A,+Az2) | 
15,;1(34+25S)\/4 his h 7} 


1 at—hs 
h | h?(e—h,—h,)*dhidh, (11) 


46754? mc! 74armcp\ ** 
- ( —) (A,+A2)F(e—2), 


6-7-15/1(3+25)|4 ae \ hk 


4 1 1 
Pa)=s(1 +—x+-1?+ x*+ “). 
2 7] 9.10 2-9-10-11 


where 
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e is the difference of energy between the begin- 
ning and end states of the nucleus in units mc’. 
A, and Az are given by (8); they contain the 
dependence of the process on the energy levels 
of the intermediate nucleus. It turns out, how- 
ever, that the values of those energies do not 
greatly influence the probability. Roughly one 
can put 


9 


LT ndTim g* 


W.-W, 





’ 
m*c4 


~ 


A,+A2= 








with g given by (5). Numerical evaluation for 
Z=31 leads to 
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BARKER 


P=1.15 X10-*°F(€—2) sec. 


= 3.6K 10° F(e—2) year“. (12) 


As seen from the general formula (11) P is 
almost independent of Z. 

The value of F for some arguments is given in 
the following table: 


e= 4 6 8 10 
F(e+2)= 0.37 X10? 9.2 «108 3.4 10° 4.2 10' 
«= 12 20 
P(e—2)= 3.3 X 108 1x10" 


The author wishes to express her gratitude to 
Professor E. Wigner for suggesting this problem, 
and for the interest taken in it. 
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The Infrared Absorption Spectra of the Linear Molecules Carbonyl Sulphide and 
Deuterium Cyanide 


Pact F. BARTUNEK AND E. F. BARKER, University of Michigan 
(Received June 28, 1935) 


Carbonyl sulphide. The discovery of new vibration- 
rotation bands in the infrared absorption spectrum of car- 
bony] sulphide has made it possible to determine the vibra- 
tional energy level scheme of the molecule. The agreement 
between theory and experiment is quite satisfactory. 
Deuterium cyanide. The fundamental bands of deuterium 
cyanide w. at 570 cm™ and w; at 2630 cm™ have been 
measured with grating spectrometers. The former has a 
strong zero branch at 570.16 cm™, and the fine structure 


HE infrared absorption spectrum of car- 
bonyl sulphide (COS) has been investi- 
gated by Cassie and Bailey' who found ten 
unresolved bands in the region 1—20u using a 
prism spectrometer. Vegard? has deduced from 
X-ray measurements that the molecule is linear 
with interatomic separations C—O of 1.10 
Angstrom units and C-S of 1.96 Angstrom units. 
The calculated moment of inertia is 178 10-*° 
gram cm’, 
The perturbed expression for 
tional energy is 


the vibra- 


E (vibration) =», Vi +v.Ve+v3V3+XiuV2 
+ X09 V2? + X33V3? +X 2ViV2t+XisViVs 
+ Xo Ve VstX, J?+ constant 


! Bailey and Cassie, Proc. Roy. Soc. A135, 375 (1932). 
2 Vegard, Zeits. f. Krist. 77, 411 (1931). 


lines on each side are well separated. From the line spacing 
the moment of inertia of the molecule is found to be 
22.92 10-*° gram cm*. A comparison of this value with 
the corresponding one for HCN, i.e., 18.72107* gram 
cm*, permits the calculation of the internuclear distances. 
From the positions of these two bands, together with those 
of the observed fundamentals of HCN, the zeroth order 
quadratic potential energy expression is computed. 


where the V’s are the vibrational quantum num- 
bers, / is the azimuthal quantum number, and 
v, to X,, are constants. If ten bands which in- 
volve these constants in an independent way are 
located experimentally the ten constants may 
be calculated and the energy level system for the 
molecule determined. This is of fundamental 
importance because it provides the correlating 
network and serves to predict new absorption 
bands. 

This problem is similar to that of HCN which 
has been developed through a number of experi- 
mental researches to a fairly complete solu- 
tion.®: 4:5. 6. 7.8 Recently Herzberg and Spinks® 


3 Burmeister, Verh. d. D. Phys. Ges. 15, 589 (1913). 

4 Barker, Phys. Rev. 23, 200 (1924). 

5 Badger and Binder, Phys. Rev. 37, 800 (1931). 

® Brackett and Liddel, Smith. Inst. 85, No. 5 (1931). 

7 Choi and Barker, Phys. Rev. 42, 777 (1932). 

8 Adel and Barker, Phys. Rev. 45, 277 (1934). 

® Herzberg and Spinks, Proc. Roy. Soc. Al47, 434 (1934). 
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have photographed two HCN bands at A10,385 
and \11,645 under high resolution. The moment 
of inertia obtained from the fine structure ob- 
served is 18.7010-* gram cm?*. By assuming 
the C-H distance to be the same as in acetylene, 
i.e., 1.06 Angstrom units, they found the atomic 
distance C—N to be 1.15 Angstrom units. These 
dimensions may be determined in a more reliable 
manner by comparing the spectra of HCN and 
DCN. We have observed two of the fundamental 
bands of the latter molecule, and have completely 
resolved one of them, thus obtaining a precise 
measure of the moment of inertia. 

Adel'® gives an analysis of vibrational isotope 
effects in three particle systems showing how 
the infrared spectra of two or more isotopic 
molecules can be applied to the analysis of their 
quadratic potential energy function making the 
assumption that the molecular force field is in- 
variant under an isotopic change of mass. We 
are now able to compute these potential con- 
stants for HCN and DCN. 

Three infrared spectrometers were available. 
The first, which is equipped with a rocksalt 
prism and gratings having 2400 and 4800 lines 
1 Adel, Phys. Rev. 46, 222 (1934). 


TABLE I. Observed bands of carbonyl sulphide. 
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per inch has been described by Barker and 
Meyer." It was used for the region 1—10u. The 
second is a recording instrument with a KBr 
prism set up by Randall and Strong” and is 
particularly adapted for rapid exploration out to 
20u. The third spectrometer has a grating with 
1200 lines per inch and a KBr foreprism, the 
optical system having been described by Hardy." 
It was used for analyzing the bands at 17.54 and 
19x. 

Two absorption cells were used. The first is 
twenty cm in length with KBr windows, and the 
second is twenty-five cm in length with windows 
of rocksalt. 

Carbonyl sulphide was prepared by treating 
50 cm’ of concentrated water solution of KSCN 
with a cold mixture of 520 grams of concentrated 
H.SO, and 400 cm’ of water. The gas COS is 
given off at 25°C (contaminated with a few 
impurities in slight amount namely CO:, CS» 
and H,S) according to the reaction;'* KSCN 
+ H.O +2H.SO,— COS + KHSO,+ NH«HSO,. 

!! Barker and Meyer, Trans. Faraday Soc. 25, 912 (1929). 

2 Randall and Strong, Rev. Sci. Inst. 2, 585 (1931). 


'S Hardy, Phys. Rev. 38, 2162 (1931). 
4 Mellor, Comprehensive Treatise in Chemistry. 
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JBSERVED BY | RAMAN LINES OBSERVED BY 
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6 2050.5 14 $21.5 6 2079 w3 14 527 w | 6 678 2w1 —2w2 
7 1892 15 $14.3 || 7 1898 2w2 +w1 15 _ 1120 . 2 524 w? 
s 1710 ] 8 1718 2w1 
if | 





TABLE II. Carbon 





OBSERVED CALCULATED ASSIGNMENT 





yl sulphide bands. 


| OBSERVED CALCULATED 








4101* 4101 2ws $14.3 516 11-20 
3739 3779 2w1i +ws RAMAN LINES 
3096 3094 2w2+ws 2233 2227 2wi +w2 
2919* 2919 wi +w3 2055 2050.5 w3 
2575* 2575 w2+ws3 1383 1380.7 wi +wr 
2050.5* 2050.5 w3 1041 1037.5 2we 
1892* 1892 2w2 +wi 859 859.2 wl 
1710* 1710 2w1 678 672.5 2w1 —2w2 
1558.6* 1558.6 3w2 524 §21.5 w2 
1530.4 1529 w3 —w? CONSTANTS 
1047.4 1037.5 2we vi = 863.4 Xu = —4.2 Xie = —2.35 Xu =2.65 
859.2* 859.2 Wi ve= 518.95 X22 = —.10 X13 =9.3 
526.6* 526.6 11-22 v3 = 2050.5 Xa; =0 Nas =3. 
§21.5* §21.5 w2 
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Fic. 1. Carbonyl sulphide bands and the deuterium cyanide band @"*. 


The products were passed from the generating 
flask through a trap immersed in a mixture of ice 
and salt which condenses the greater part of the 
CSe, then through a 33 percent solution of KOH 
to precipitate CO, and H.S, and finally through 
traps containing CaCl, and P,O;. The dry gas 
was then condensed with liquid air and stored in 
a five liter reservoir. 

The deuterium cyanide was supplied by the 
California Isotope Company in a small glass 
capsule containing a mixture of about 40 percent 
DCN and 60 percent HCN. 

The readings were taken in the usual in and 
out method for all the bands with the exception 
of the 17.54 band of DCN, the 14u band of HCN 
and the 19u region of COS. In these cases the 
cell was in the beam permanently during the time 
deflections were taken. It was then removed and 
the ‘“background”” mapped on the same day. 
Care was always taken of course to insure a 
pure spectrum. The Nernst glower was used as a 
source of continuous radiation. The calibration 
of the spectrometers is in terms of the water 
vapor lines in the region of 6y. 

The results for carbonyl sulphide are shown in 
Fig. 1 and Tables ' and II. Table I gives the 
observed bands in terms of wave numbers with 


the bands observed by Cassie and Bailey! for 
comparison. At the end of the table the Raman 
lines observed by Dadieu and Kohlrausch” are 
listed. 

The 2.674 band is not reproduced. It is very 
irregular due to superposition of the CO, and 
H.S bands in the same region since it is impos- 
sible to remove every trace of these substances 
from the gas. It is extremely difficult to judge 
the position of the center but it seems to be 
at 3739 cm~. Cassie and Bailey! found a band 
in this region which they list as 3742 cm™. This 
band was not used in calculating the constants. 
The two bands 1558.5 cm and 1530 cm™ are 
also not reproduced. Both have sharp zero 
branches similar to the 2575 cm~! band. Two 
of the three zero branches in the 19y region are 
upper stage transitions from the first excited 
state to the second excited states. These are 
clearly shown on the photographic record, Fig. 1. 
The proportion of molecules in the first ex- 
cited state at room temperature is 17 percent 
which is sufficient to account for the upper stage 
absorption. 

Table II gives the observed and calculated 
positions of the infrared bands and Raman lines, 


1 Dadieu and Kohlrausch, Physik. Zeits. 33, 165 (1932). 
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Fic. 2. Rotation lines of the 17.54 region of deuterium cyanide. 


the calculated values of the ten constants being 
obtained from the bands marked with asterisks. 
These constants are given at the end of the 
table. 

The results for deuterium cyanide are shown 
in Fig. 2 indicating the fine structure of the 
17.54 region; in Fig. 1 where the fundamental 
band w; appears; and in Tables III and VY. 
Table III gives the observed and calculated 
positions of the rotation lines of DCN while 
Table IV gives similar results for HCN. The 


calculated values in these tables were obtained 
from the formula v=vp+Am+Bm?+Cm' the 
constants vo, A, B and C being obtained empiri- 
cally. The ordinal number ‘‘m’’ is counted from 
the band center. The constants are for DCN 
vo =569.10; A =2.425; B=0.0013; C= —0.00003 
and for HCN») =712.07; A =2.952; B=0; C=0. 

The part of the band immediately to the low 
frequency side of the zero branch is somewhat 
irregular due to the superposition of the upper 
stage bands of DCN. The rather sharp intense 


TABLE III. Rotation lines of the 17.54 band of DCN zero branch 570.16 em™, 








R Brancn AJ=-+1 P Braxncn AJ=-—-1 R Brancu AJ=-+1 P Brancn AJ 1 
Os. CaALc. Res. m Ors. CALc. REs. ORs. CALC. Res. m Ons. CALC Res 
571.53 1 600.76 600.78 —.02 13 537.86 §37.87 —.01 
573.74 573.96 —0.22 2 564.52 564.26 0.26 603.17 603.22 05 14 535.53 535.48 0s 
576.28 576.39 —.11 3 561.81 561.83 —.02 605.66 605.67 —.01 15 533.17 §33.11 06 
578.84 578.82 02 4 559.41 559.42 —.01 608.05 608.11 — .06 16 530.83 530.75 OS 
581.06 581.26 —.20 5 556.98 557.00 —.02 610.45 610.54 — .09 17 528.36 528.40 — O04 
583.67 583.69 — .02 6 554.62 554.61 O01 612.94 613.00 — .06 18 526.02 526.04 — 02 
586.25 586.13 -12 7 $52.29 §52.19 .10 615.60 615.44 .16 19 523.86 §23.71 .15 
588.59 588.56 03 8 549.98 549.80 .18 617.99 617.88 11 20 §21.51 §21.36 15 
591.15 591.02 .13 9 547.69 547.41 .18 620.56 620.32 .24 21 
593.60 593.45 mi 10 545.10 545.01 .09 622.71 622.76 —.05 22 
595.94 595.90 O04 11 $42.59 §42.63 —.04 624.88 625.21 —.33 23 
598.31 598.34 — .03 12 540.37 $40.24 .13 
TABLE IV. Rotation lines of the 14u band of HCN zero branch 712.28 cm™. 
R BrancHh AJ=+1 P BRancnh AJ=—1 R Brancu AJ =-+1 P Brancn AJ 1 

Oss. CALc. Res. m OBs. CaLc. REs. Oss. Cale. Res. m Orns. CaLe RES 

- 715.02 — 1 - 753.35 753.40 — .05 14 670.69 670.74 05 
717.86 717.97 —0.11 2 706.33 706.17 0.16 756.29 756.35 — .06 15 667.66 667.79 13 
720.81 720.93 —.12 3 703.23 703.21 02 759.25 759.30 —.05 16 664.606 664.84 —.18 
723.78 723.88 —.10 4 699.60 700.26 —.66 762.23 762.25 —.02 17 661.69 661.89 20 
726.77 726.83 —.06 5 697.25 697.31 —.06 765.24 765.21 03 18 658.75 658.93 18 
729.80 729.78 02 6 694.24 694.36 —.12 768.28 768.16 -12 19 656.08 655.98 10 
732.85 732.73 12 7 691.38 691.41 — .03 770.99 771.11 —.12 20 653.19 653.03 16 
735.93 735.69 24 & O88.42 O88.45 — .03 is. 74.07 774.06 01 21 650.33 650.08 25 
738.72 738.64 OS 9 685.48 685.50 —.02 777.18 77701 17 22 647.02 647.13 11 
741.85 741.59 6 10 682.57 682.55 02 779.97 779.97 00 23 644.17 
744.69 744.54 15 11 679.69 679.00 09 782.78 782.92 14 24 641.22 
747.55 747.49 06 12 676.57 676.05 — .O8 785.84 785.87 03 25 638.27 
750.44 750.45 01 13 673.49 673.609 —.20 788.81 788.82 Ol 26 
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TABLE V. Average separations of the rotation lines. 


HCN 14 Banp 
R(m) —P(m) 


DCN 17.54 BAND 














Rim) —P(m) 
2m m 2m 

2.305 cm=! 2 2.883 cm7! 
2.412 3 2.930 
2.429 4 3.023 
2.408 5 2.952 
2.421 6 2.973 
2.426 7 2.962 
2.413 & 2.969 
2.414 9 2.958 
2.425 10 2.964 
2.425 11 2.955 
2.414 12 2.958 
2.419 13 2.959 





DCN 17.54 Bann 


HCN 144 BAnp 








R(m) —P(m) R(m) —P(m) 
2m m 2m 
2.416 14 2.952 
2.416 15 2.954 
2.413 16 2.956 
2.414 17 2.957 
2.414 18 2.958 
2.414 19 2.953 
2.412 20 2.954 
-- 21 2.946 
— 22 2.958 


Average 2.955 cm~! 


Average 2.414 cm7! 
J =18.72 X10~® g cm? 


=22.92 X10~-” g cm? 








line at 561.07 cm™ is probably the zero branch 
of the transition 1,—2») where the subscript indi- 
cates the / value. The upper stage transition 
1,—22 probably lies somewhere within the zero 
branch of that from the ground state 0 —1). 
At room temperature about 14 percent of the 
molecules are in the first excited state. This is 
sufficient to account for the presence of upper 
stage zero branches. 

Table V gives the average separations of the 
lines from which the moment of inertia about 
the center of gravity is calculated. Values of 
(R(m)—P(m))/2m are given for both HCN and 
DCN, with the calculated moments of inertia. 

A search was made for the other two funda- 
mentals w; and ws; as well as for harmonics. 
Because of the rather small amount of gas avail-’ 
able the only other band observed was w3 at 
2630 cm™. 

In the case of HCN because of the strong 
absorption of carbon dioxide in the background, 
readings could not be taken for a short space as 
the deflections dropped to zero. Because of the 
regularity of the lines the missing line was inter- 
polated and found to agree very well with the 
calculated value. 

The expression for the moment of inertia of 
this linear triatomic molecule about its center of 
gravity is 


T=m,x?+my* 
1 





(m?x?+ me2?y? — 2m ymexy) 
mi+me+ms; 


where x is the H—-C distance and y is the C-N 
distance. m,, m2 and m3; are the masses of the 
‘hydrogen, nitrogen and carbon atoms, respec- 
tively. Introducing the isotopic masses m, and 


m,’ and solving the two equations simultane- 
ously we obtain the interatomic separations, 
x=1.0610-§ cm and y=1.15 10-5 cm. 

The theory of small oscillations may be ap- 
plied, considering the molecule as a group of 
uncoupled simple harmonic oscillators of which 
two are linear and one is plane isotropic in a 
plane perpendicular to the axis of the molecule. 
The normal modes of vibration are well known. 
We shall use coordinates somewhat different 
from those defined by Adel'® in order to pre- 
serve the same meaning of these coordinates for 
both isotopes. Let 7 represent the distance from 
the central atom to the line joining the end atoms 
measured parallel to the displacement of the 
particles in the perpendicular vibration. Let q 
measure the relative separation of the hydrogen 
and carbon atoms, and z measure the relative 
displacement of the carbon and nitrogen atoms. 
In the perpendicular vibration let y:, ye and y3 
represent the displacements of the hydrogen, 
nitrogen and carbon atoms. Then the kinetic 
energy is 

T= mY r+ 2myet+ mys. 
The expression for 7 in terms of y;, ye and y; and 
the internuclear distances is 


r= i¥/ro+yex/ro—¥s Where rm=x+y. 


Since the momentum is conserved as well as the 
moments of momenta about any point one has 
two additional relations 

Mii +mye+myYs=0 and miyix—meyoy=0. 


In terms of 7 the kinetic energy becomes T = }u,i* 
where 


1yx\? 1yy\? 1 
ASA) 
Me2X\7o m,\ To Ms 
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For the longitudinal vibrations x, x2 and x3 repre- 
sent the displacements from equilibrium of the 
hydrogen, nitrogen, and carbon atoms, respec- 
tively. Then 


Q=X3— Xi; S=%2.— Xs; 
and 

MX, + Meke+m3%3 = 0. 
In this case the kinetic energy is T=}m,é? 
+ 3me%2? +3m3%;" or in terms of g and z, T=} u,¢* 
+3u.2°+}-2u,.g2 where 


m,(m2+ms3) mo(m,+ms3) m\Ms 
= Ss oe es Kqz:> 


S S S 








and S=m,+m.+ms3, the total mass of the mole- 
cule. The potential energy is subject to the 
geometric symmetry of the molecule. It may be 
written V= (Rig? t+keor?+ 2kisg2+k3327). The 
secular determinant of the normal coordinate 
transformation factors into a linear and a quad- 
ratic expression, and yields the following set of 
equations, 


MoRsatmekii— 2g 2k13 
Ai +A3= : ’ 
Mgkz—Maz 


Risks3— kis? 





AiA3= 
MoMe—Ma? 


Ao = koo/ pr. 
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In normal coordinates the kinetic and potential 
energies assume the forms 


T=VO+R+4Z2), V=}OiQ+eR?+2s2?). 


The normal vibrations expressed in wave num- 
bers are v;=(1/2mc)(A;)'. The constant ke. may 
be computed from we of HCN and DCN and is 
nearly the same for both as would be expected. 
Using the fundamentals w; and w; of HCN and 
the fundamental w;’ of DCN the unobserved 
fundamental w,;’ of DCN can be calculated as 
well as the other three constants kj), ki3 and 33. 
There are two sets of constants because the 
equations are quadratic. These are 
ke. = 0.8071 105 ko= 0.8071 105 
ku =5.6987 * kui= 5.6981 “ 


O 
kis=6.443 "his = —0.30491 “ 
ky=25.34 k= 18.065 “ 


dynes/cm 


ki3 represents the interaction between the ex- 
treme atoms, and should be small. Hence the 
second set is the more plausible. 

The observed and calculated values of the 
wave numbers for the fundamental bands of 
HCN and DCN are as follows: 


Observed wave numbers Calculated wave numbers 


w, =2089.0 cm w, = 2089.6 cm™ 
w, = 712.1 w, = 712.8 
w;3 = 3312.9 ws =3312.3 
w= — wit = 1896.7 
we'= 569.1 w'= 568.6 


ws! = 2630.0 w;! = 2629.5. 
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Absorption Spectrum of Tin Sulphide 


E. N. SHAWHAN, Mendenhall Laboratory of Physics, Ohio State University 
(Received July 8, 1935) 


The absorption spectrum of tin sulphide between 3200A 
and 4600A has been photographed at moderate dispersion 
and between 4000A and 4600A at high dispersion. The 
system in the near ultraviolet may be represented by the 
formula: 


v = 28,336.2+332.77(v’ +3) —1.42(v’ +4)? 
—488.25(v" +3) +1.47(0"+4)°. 


This formula is not in agreement with that given by 


Butkow and Tschassowenny who seem to have incorrectly 
assigned the system origin. A new system has been located 
in the visible which is limited by predissociation to two 
v’’ progressions. These bands are given by the formula: 


v = 23,287.67 +367.5(v'+ }) —488.25(v"’ +3) +1.47(0''+})?. 


’ The observed isotope structure of the strong bands in this 


system confirms their vibrational assignment. 
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INTRODUCTION 


UR knowledge of the band spectra of the 

metallic sulphides seems quite meager when 
compared with the information available on the 
spectra of the oxides of the metals. This is 
perhaps in part due to the difficulties involved in 
obtaining spectra of sulphides in emission. Some 
progress has been made in several cases where it 
is possible to obtain absorption spectra. Wilhelm! 
has reported the spectrum of MgS and given its 
vibrational analysis. The absorption spectra of 
ZnS, CdS and HgS have been studied by Sen- 
Gupta.’ He observed only regions of continuous 
absorption, which indicates that the potential 
curves of the excited states of these molecules 
either have shallow minima or none at all. 
Shapiro, Gibbs and Laubengayer*® made a vibra- 
tional analysis of the absorption spectrum of 
GeS and identified heads due to four isotopes of 
germanium. Rochester and Howell‘ have given a 
rather complete vibrational analysis of the ex- 
tensive absorption spectrum of PbS and indicate 
that the rotational analysis is in progress. 

A short report on the vibrational analysis of 
the band spectrum of SnS in absorption has been 
published by Butkow and Tschassowenny.°® They 
report only two band systems, one lying between 
2580 and 3164A and the other between 3277 and 
3799A. They give the band-head formula in each 
case without recording data on the bands ob- 
served. One would expect the spectrum of SnS to 
resemble that of PbS. The fact that the spectrum 
of PbS extends into the near infrared and includes 
at least five systems suggested that a further 
critical examination of the spectrum of SnS 
might reveal other systems lying farther toward 
the visible. 

EXPERIMENTAL 


Chemically pure tin sulphide (SnS) contained 
in a graphite tube about 50 cm long and 2 cm in 
internal diameter, open to the air at the ends, 
was heated by means of a Hoskins carbon 
resistance furnace operated at 25 volts a.c. and 


1H. A. Wilhelm, Iowa State Coll. J. Sci. 6, 475 (1932). 

2 P. K. Sen-Gupta, Proc. Roy. Soc. A143, 438 (1934). 

’C. V. Shapiro, R. C. Gibbs and A. W. Laubengayer, 
Phys. Rev. 40, 354 (1932). 

4(4. D. Rochester and H. G. Howell, Proc. Roy. Soc. 
A148, 157 (1935). 

°K. Butkow and W. Tschassowenny, Zeits. f. Physik 90, 
53 (1934). 
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about 300 amperes. The temperature was varied 
between 800 and 1100°C. A five-hundred watt 
projection lamp provided a source of continuous 
radiation which was passed through the tin- 
sulphide vapor and then focused on the slit of the 
spectrograph. 

The spectrum between 3200A and 4000A was 
photographed on an E-1 quartz spectrograph 
giving a dispersion of 6 to 12A/mm over this 
range. The times of exposure with this spectro- 
graph ranged up to a half-hour. A_ Littrow 
mounted concave grating with a radius of curva- 
ture of ten feet, giving a dispersion of 5.6A/mm 
in the first order was used to photograph the 
region between 3400A and 4650A. In order 
properly to develop the bands near the red end 
of the spectrum it was necessary to increase the 
vapor pressure until the light on the slit appeared 
distinctly yellow. At a somewhat lower tempera- 
ture the violet end of the spectrum could be 
photographed in about an hour. 

For still higher dispersion between 4000A and 
4650A this portion of the spectrum was photo- 
graphed in the first order of a 21-ft. concave 
grating in a Paschen-Runge mount. This grating 
is ruled with 30,000 lines to the inch and gives a 
dispersion of 1.25A/mm in the first order. The 
SnS was introduced into an iron tube approxi- 
mately 80 cm long and 1 cm in internal diameter. 
This tube was heated to about 1000°C in two 
nichrome-wound Hoskins furnaces placed end to 
end. As in the previous case a five-hundred watt 
projection lamp provided the continuous radi- 
ation which was passed through the iron tube 
containing the SnS vapor. 


RESULTS 


Two band systems were observed, one in the 
ultraviolet extending from about 3250A to 
4300A, and the other in the visible between 
approximately 4200A and 4600A. The system in 
the near ultraviolet probably includes the bands 
observed in this region by Butkow and Tschasso- 
wenny.°® Their near ultraviolet system extended 
as far toward shorter wave-lengths but not as far 
to the red. The system in the visible does not 
seem to have been previously reported. It is 
confined to a few fairly strong bands degraded to 
the red as are all other bands in the spectrum. 
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The rotational structure of the stronger bands is 
resolved on the high dispersion plates. 

A number of weak heads, scattered through the 
portion of the spectrum explored, were measured 
but could not be classified. Those falling in the 
region photographed on the 21-ft. grating were 
examined on the high dispersion plates. Most of 
them were found to be due to the accidental 
grouping of rotation lines occurring in the 
complex structure resulting from the overlapping 
of bands corresponding to the numerous isotopes 
of tin. The few unclassified heads that cannot be 
accounted for in this manner may belong to 
fragmentary systems. A number of unclassified 
heads® have also been observed in the spectrum 
of SnO. 


Near ultraviolet system 


Table I gives the wave numbers of the heads 
in the near ultraviolet system, together with 
their estimated intensities, the values of v’ and v” 
and the departures of the observed wave numbers 
from those calculated from the formula: 


v= 28,336.2 +332.77(v' +3) —1.42(v'+})* 


— 488.25(0'° +3) +1.47(0'"+3)?. 

At a temperature of about 1000°C this system 
is characterized by strong v’ and v” progressions 
which allow the values of v’ and v” to be assigned 
with confidence. At lower temperatures the bands 


TABLE I. Bands observed in near ultraviolet system. 








Yobs p> Yobs = 

I v 9” "obs. Veale. I = Yobs ¥ealc 
0 1 10 23865.5 —2.1 2 4 3 27637.3 —0.3 
2 1 9 24329.0 +2.5 10 0 1 27773.1 —0.1 
0 3 10 24520.6 +1.6 4 3 3 27791.8 —0.9 
0 2 9 24650.1 —2.5 1 5 4 27957.8 +1.6 
0 1 8 24788.3 0.0 10 1 1 28102.9 —0.2 
1 0 7 24923.6 +0.5 4 0 0 28259.1 +0.6 
0 3 9 24975.2 —2.6 4 3 2 28271.4 —0.9 
1 2 & 25115.2 —0.2 8 2 1 28430.2 0.0 
3 1 7 25253.8 +0.8 0 10 6 28559.7 —2.5 
1 4 9 25296.3 —2.9 8 1 0 28589.2 +0.8 
5 0 6 25391.5 +0.8 6 4 2 28594.3 +0.8 
6 1 6 25720.5 —0.2 1 3 1 28753.4 —1.1 
7 0 5 25862.8 +0.5 7 2 0 28915.8 +0.3 
4 1 5 26190.7 —0.6 1 4 1 29075.2 —0.7 
8 0 4 26337.5 +0.6 6 3 0 29239.5 —0.3 
3 3 6 26369.4 —2.6 2 5 1 29393.4 —1.0 
1 2 5 26517.7 —0.7 1 10 4 29504.5 —1.8 
0 5 7 26548.5 —2.8 5 4 0 29560.5 —0.7 
10 0 3 26811.0 —0.1 2 9 3 29677.7 —0.8 
3 3 5 26845.6 —1.5 5 6 1 29710.3 +0.1 
4 2 4 26989.2 —1.5 1 & 2 29852.0 +1.3 
0 1 3 27139.8 —1.5 5 5 0 29880.1 +0.3 
1 4 5 27167.3 +2.7 3 7 1 30022.1 —1.0 
10 0 2 27291.7 +0.9 1 9 2 30159.8 +1.8 
2 3 4 27315.7 —0.5 3 6 0 30195.5 0.0 
5 2 3 27468.5 +0.1 2 8 1 30332.8 —0.3 
5 1 2 27621.0 +0.2 1 7 0 30508.2 —0.2 





: 
| 
| 


® F.C. Connelly, Proc. Phys. Soc. 45, 780 (1933). 





TABLE II. Jntensity distribution in near ultraviolet system. 








- 0 1 2 3 4 5 6 7 8 9 10 

0 4 10 10 10 8 7 5 1 

1 10 10 5 0 4 6 3 0 2 0 
2 10 8 5 4 1 1 0 

3 10 1 4 4 2 3 3 0 0 
4 8 1 6 2 1 1 

5 5 2 1 0 

6 3 5 

7 1 3 

8 2 1 

9 1 2 
10 1 0 








corresponding to the higher v’’ values are notice- 
ably weaker. Higher temperatures failed to reveal 
additional bands. 

Watson and Loomis’ found difficulty in meas- 
uring the heads of the bands in the far ultraviolet 
system of SnO because of the overlapping of 
bands corresponding to the five abundant and 
six less abundant isotopes of tin. The values of p 
for the various tin isotopes are larger for SnS than 
for SnO, but most of the bands in this system are 
so near the system origin that the heads corre- 
sponding to the several isotopes are nearly 
coincident and the observed heads are fairly 
sharp. Near the limits of the system the heads 
appear somewhat diffuse but the spreading is not 
sufficient to allow the individual isotope heads to 
be distinguished. 

The intensity distribution of the bands in this 
system is shown in Table II. There is a rather 
open Condon parabola to be expected from the 
difference in the values of w,’ and w,’’. Inside this 
parabola, and separated from it by weaker 
bands, is another group of fairly strong bands. 
This effect* is quite marked in the A system of 
PbS. 

Comparison with earlier work. Butkow and 
Tschassowenny® gave the following formula for 
the bands between 3277A and 3799A: 


v= 30,280.5+315.6(v' +3) —0.59(v' +3)? 
—484.5(0"’ +3)+0.50(0'’ + })?. 


This formula accounts for only a few of the 
strong bands observed in this region of the 
spectrum during this investigation. It may be 
that the lower temperatures used by these 
investigators, 600 to 840°C, were not sufficient to 
develop properly the bands near the red end of 


7F. W. Loomis and T. F. Watson, Phys. Rev. 45, 805 
(1934). 
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the system, leading to an incorrect assignment of 
the system origin, which necessarily gives an 
incorrect value for w,’. The (0,0) band of the 
formula as given by Butkow and Tschassowenny 
with origin at 30,280.5 becomes the (6,0) band of 
the system proposed here with vy, as 28,336.2. 
Their formula fits only bands which may be 
assigned values of v’ equal to or greater than 6 in 
our formula. 

Butkow and Tschassowenny represented the 
bands between 2580A and 3164A by the formula: 


y = 34,868.3 + 288.3(0’ +3) —1.255(0’ +3)? 
— 484.5(0' +3) +0.50(0" + 3)°. 


At the temperatures used in this investigation 
SnS vapor is almost opaque to light in the far 
ultraviolet. No attempt was made to repeat their 
observations in this region. 


Visible system 


In the vicinity of 4300A there are three strong 
bands and six less intense, all degraded toward 
the red. These bands may be represented by the 
formula: 


v = 23,287.67 + 367.5(v' +3) 


—488.25(0" +3)+1.47(0"+ 3). 
The vibrational assignments of the bands, esti- 
mated intensities, and differences between the 
observed and calculated wave numbers are given 
in Table III. No bands have been observed which 


TABLE III. Bands in visible system. 








Y’ 0 1 2 3 4 5 
(10) (10) (6) (1) 
23531.89 23046.51 22563.78 22084.34 
(+0.09) (+0.02) (-—0.34) (—0.35) 
(10) (2) (4) (1) (0) 
1 23899.34 22931.78 22452.38 21976.65 21501.87 
(+0.04) (+0.16) (+0.19) (—0.05) (—0.28) 





seem to correspond to values of v’ higher than 1. 
In view of the fact that the (0,2) band is fairly 
strong, this can be explained only by supposing 
that the potential curve of the upper state of this 
system is crossed by a repulsive state just beyond 
v’ = 1. The conclusion that there is predissociation 
is supported by the fact that the rotational 
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structure of the (1,0) band seems to terminate 
abruptly about 200 cm"! from the head. 

The vibrational assignments of the three 
strong bands are qualitatively confirmed by the 
observed isotope structure. The lines of the 
(0,0) band are sharp and single near the head, 
but become broader and appear to split into 
several partially resolved components at a con- 
siderable distance from the head. On the other 
hand, the lines of the (1,0) band are rather diffuse 
near the head but become sharper toward the 
tail of the band. The structure of the (0,1) band 
is quite complex. The displacements of the lines 
of any isotopic molecule from those of the most 
abundant molecule may be expressed by the 
approximation, 


v'—pvy=(p—1)v,+2(p—1)p,. 


The value of x,’w,’ is not known but it is certainly 
small and the system origin must be in the 
vicinity of 23,580 cm~'. Hence v, and v, have 
opposite signs for the (1,0) band and the same 
signs for the (0,0) and (0,1) bands. The lines of 
all isotopes should coincide because of the 
cancelation of the rotational and vibrational 
displacements at about 160 cm™ from the head 
of the (1,0) band. This cancellation does not 
occur in the (0,0) and (0,1) bands. The lines of 
the (0,0) band are fairly sharp because of its 
proximity to the system origin but become 
broader at increasing distances from the head due 
to the combined rotational and_ vibrational 
displacements. 

The rotational analysis of the bands in this 
system is in progress and will be the subject of a 
later report. 

Grateful acknowledgment is made to Professor 
Alpheus W. Smith and Dr. R. V. Zumstein for 
their many helpful suggestions. 

Note added in proof: While this paper was in 
the press Rochester (Proc. Roy. Soc. A150, 668 
(1935)) also reported on the band spectrum of 
stannous sulphide. He assigns the origin of the 
visible system at 23,211.8 cm! whereas our 
analysis indicated that it is at about 23,580 
cm~!, Our data on these band heads taken from 
the high dispersion spectrograms are incon- 
sistent with his assignment. 
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A mass spectrograph, designed for the identification of the 
products of ionization and dissociation by electron impacts 
in molecular gases is described. 

Ionization potential of nitrogen. .\ direct comparison in a 
nitrogen-argon mixture showed that the ionization poten- 
tial of nitrogen is at least 0.04 volt less than that of argon. 
The ionization potential is therefore not greater than 15.65 
volts, though it may be lower by as much as 0.1 volt. 

Analysis of ionization and dissociation products in C,H», 
NO, and C:N,. The ions formed by electron impacts in these 
gases together with their appearance potentials in volts 
are as follows: 


Ion C:H2* C:H*t C2* CHt c Ht 
App. pot. 11.2 17.8 23.8 22.2 24.5 21.7, 25.6 
Est. error +0.1 +0.2 +0.3 +0.5 +1.0 +10 +1.0 
Ion NOt Nt NO*+ 

App. pot. 9.5 22.0 44 

Est. error +0.1 +0.5 +1.0 

Ion C2Ne* C:N*t CN* C2* 

App. pot. 14.1 19.8 21.3 18.6 

Est. error +0.1 +0.5 +0.3 +0.5 


Negative O~ ions appear in NO as reported by Tate and 
Smith. An ion of mass 27 found in C.Hz is ascribed to 
C®8C4H,. The abundance ratio C'® : C® was calculated to 
be 1 : 100. 





I. Tue Mass SPECTROGRAPH 


N THE studies here to be described of the 

effects of electron impacts in various gases a 
mass spectrograph was used to identify the ion- 
ized products of the impacts. 

The mass spectrograph which is shown dia- 
grammatically in Figs. 1 and 2 is similar to that 
described by Bleakney.' Since many of the gases 
to be studied are dissociated at the hot filament 
it was necessary to shut off the filament, except 
for the first small diaphragm in the cylinder A, 
in a separate compartment provided with a 
separate pump. The scheme is shown in Fig. 1. 
Since the pressure is kept lower on the filament 
side of A than in the main tube the products of 
dissociation were kept from diffusing into the 
region where impact ionization and dissociation 
were being studied. 


i IW, Bleakney, Phys. Rev. 40, 496 (1932). 


! | to pump 


A magnetic field of about 900 gauss parallel to 
the axis of the tube was provided by a large 
solenoid surrounding the tube and coaxial with it. 
This field served the double purpose of maintain- 
ing the electron beam from the filament in line 
and of providing the magnetic field for the 
analyzer. 

Electrons from the tungsten filament F are 
accelerated by a constant drawing out potential 
between Fand A. S,;is at the same potential as A. 
A variable potential V, is applied between the 
filament F and the slits Ss, S; and S; which are 
all at ground potential. The electrons with 
energy of approximately V,, electron-volts then 
enter the region between the plates P; and P» 
which form part of the analyzer shown in end 
section in Fig. 2. Here a few of them suffer im- 
pacts with the gas molecules but the majority go 
straight through to the electron trap which is, in 


DOS S,54 : R, 
—_ - ] iV roa (UO 4-2 _ 
st — * ; 
BR? CF | 
ade 3 e Llectron trap 


Analyser 


Fic. 1. Diagram of the apparatus as viewed from the side. 
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ic. 2. Diagram as viewed from the end with electron 
gun removed. S;, So, Ss; and S;—1X6 mm; S;—0.2 x50 
mm; Ss and Sy—0.155 X50 mm. 


turn, connected with a galvanometer. The pres- 
sure of the gases to be studied is so low (10~4 to 
10-* mm Hg) that a very small fraction of the 
electrons experience single impacts, and multiple 
impacts are negligible. 

The ions produced in the region between P; and 
P, are drawn downward by a small potential Vo 
(about 2 volts) applied between P; and P2. This 
potential is so applied that the space through 
which the electron beam passes is maintained at 
ground potential. There is a long slit in P»2 (Fig. 2) 
through which a narrow sheet of ions may pass. 
These are further accelerated by a variable poten- 
tial V; between P2 and P;, are bent around in the 
magnetic field and pass through the slit S, to the 
collector P;. The current to P; is measured with 
a Compton electrometer or vacuum tube ampli- 
fier. The resolving power of the analyzer (ratio of 
A(m/e) to (m/e) is about 1/40). 

All of the metal parts of the spectrograph are 
of tantalum or tungsten with Pyrex insulation. 
The tube is mounted in an electric furnace which 
in turn is mounted in the center of a water- 
cooled solenoid. 

The tube was evacuated and baked at a tem- 
perature of about 400°C, for many hours and the 
baking process was repeated at frequent intervals. 
The spectrograph itself forms a good ionization 
gauge. The residual pressure, before admitting 
the gas to be studied, was less than 10-7 mm Hg. 
The gases were admitted through a long capillary 
and allowed to flow continuously out through the 
pumps. 

In making an m/e analysis of the ions their 
velocity was varied by varying V; in steps of 0.1 
or 0.05 volt and the current to P; plotted as a 
function of V; (see Fig. 3, for example). When the 
appearance potential of a particular type of ion 
is to be determined J; is set for the collection of 
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that ion and the electron velocity varied by 
varying V’,. The ion current is then plotted as a 
function of V,. This type of curve is shown in 
Fig. 4. 


II. IoNIZATION POTENTIAL OF NITROGEN 


The precise value of the ionization potential of 
the nitrogen molecule has for several years been 
a matter of some concern both to those interested 
in electron collision phenomena and to band 
spectroscopists. Determinations of the threshold 
energies electrons must have in order that they 
be capable of producing ions in a gas are subject 
always to uncertainties introduced by the shape 
of the ionization curves near the threshold. These 
uncertainties have two causes: the first a lack of 
sharpness in the velocity distribution of the 
electrons themselves and the second a lack of 
knowledge about the actual probability-of- 
ionization as a function of velocity in the neigh- 
borhood of the threshold. 

The first cause of uncertainty can be removed 
by an independent determination of the velocity 
distribution but the second is more refractory and 
is particularly troublesome in the case of mole- 
cules for which the sharpness of onset of ioniza- 
tion depends on the shapes of the potential-energy 
curves for the molecule and molecular ion and the 
relative positions of their minima with respect to 
internuclear separations. The result is that the 
estimate of ionization potential will depend on 
the sensitiveness of the measuring instruments 
and any value so obtained usually represents an 
upper limit. 

The present determination of the ionization po- 
tential of nitrogen was made by using a mixture 
of nitrogen and argon in the mass spectrograph 
and by measuring simultaneously the number of 
N.*+ and of A+ ions collected for a number of 
electron energies near the threshold energy. The 
advantage of using argon as a calibrating gas is 
that its ionization potential differs from that of 
nitrogen by only about 0.1 volt yet the mass 
spectrograph permits a complete separation of 
the At and N¢* ions and a determination of their 
relative abundance. 

This method of estimating the ionization po- 
tential of nitrogen eliminates all uncertainty rela- 
tive to contact potentials and velocity distribu- 
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Fic. 3. Mass spectrum of the products of ionization in a mixture of argon and nitrogen. Electron 
energy V,=100 EV. H=900 gauss. 


tion of electrons. There remains only the uncer- 
tainty about the relative behavior of the ioniza- 
tion-probability curves for argon and nitrogen. 
It is probable that the onset of ionization in argon 
is somewhat sharper than in nitrogen. 

The mass spectrum of an argon-nitrogen mix- 
ture is shown in Fig. 3. The heights of the A* and 
of the N.* peaks are plotted against electron 
velocity in Fig. 4. It will be seen that the ioniza- 
tion potential of nitrogen is a little lower than 
that of argon. Twelve curves of the sort shown 
in Fig. 4 give for the difference in ionization 
potentials an average value of 0.04 volt. The 
mean deviation from this average was 0.04 volt. 
If we assume that the ionization-potential of 
argon is 15.69 volts, the apparent value for nitro- 
gen is 15.65 volts. We believe that the true ioniza- 
tion potential is certainly not greater than this 
by more than 0.02 volt but may be less by per- 
haps 0.10 volt. 

In Table I are given the values of previous 
measurements of the ionization potential of 
nitrogen. The fact that electron impact measure- 
ments lead to upper limits is indicated by the 
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Fic. 5. Mass spectrum of the products of ionization in acetylene. Electron energy 1,=100 EV. 


way in which, as the experimental technique is 
refined the values given become lower and lower. 


TABLE I. Observed ionization potential for molecular nitrogen. 


(eoite) OBSERVER REFERENCE 

18 Davis and Goucher Phys. Rev. 13, 1 (1919) 

18 Smyth Phys. Rev. 14, 409 (1919) 
16.9 Mohler and Foote Phys. Rev. 15, 555 (1920) 
17.2 Stead and Gosling Phil. Mag. 40, 413 (1920) 
15.8 Found Phys. Rev. 16, 41 (1920) 
16.95 Brandt Zeits. f. Physik 8, 32 (1921) 
15.8 Boucher Phys. Rev. 19, 189 (1922) 
16.2 Duffendack Phys. Rev. 20, 665 (1922) 
16.3 Mackay Phys. Rev. 24, 319 (1924) 
15.8 Turner and Samson’ Phys. Rev. 34, 747 (1929) 
15.8 Vaughan Phys. Rev. 38, 1687 (1931) 
15.7 Tate and Smith Phys. Rev. 39, 270 (1932) 
15.65 Tate, Smith and This paper. 


Vaughan 


Recently Mulliken? concludes from Hopfield’s* 
study of a Rydberg series of absorption bands 
that 15.51 volts is a minimum value for the 
ionization potential of nitrogen. We can therefore 
assert with some confidence that the ionization 
potential lies somewhere between 15.51 and 
15.65 volts. 

As a matter of general interest we looked for 
negative ions of nitrogen but found none at any 
velocity of impact. 


*R.S. Mulliken, Phys. Rev. 46, 144 (1934). 
3 J. J. Hopfield, Phys. Rev. 36, 789 (1930). 


III]. PRopucts OF IONIZATION AND DISSOCIATION 
PRODUCED BY ELECTRON IMPACT IN ACETY- 
LENE, NITRIC OXIDE, CYANOGEN AND 
CARBON MONOXIDE 


A. Acetylene 

The gas was prepared by repeated fractional 
distillation of the product. The 
measurements were made with the acetylene at a 
pressure of about 10~> mm Hg. The mass spec- 
trum of the products of ionization and dissocia- 
tion are shown in Fig. 5. Eisenhut and Conrad‘ 
found the same primary ions. The appearance 
potentials and estimated accuracy are obtained 
from curves of the type shown in Fig. 6 and are as 
follows: 


commercial 


Ion C;H.+ C:-H+ C.+ CH*+ Ct H* 
App. pot. 

(volts) 11.2 17.8 23.8 22.2 245 21.7 25.6 
Est. error +0.1 +0.2 +03 +05 +1.0 +1.0 +1.0 


It will be noted that the peaks for the C+, CH* 
and H* ions trail out to the left. This is due to the 
fact that in the dissociation process the ions pick 
up kinetic energy. The energy lost by the electron 
at impact is greater by several electron volts than 
the sum of the ionization potential and the heat 
of dissociation and the difference appears as 
kinetic energy of the dissociation products. This 


4Q. Eisenhut and R. Conrad, Zeits. f. Elektrochemie 36, 
654 (1930). 





[ON 


nal 
The 
ita 
eC - 
‘ja- 
ad‘ 
nce 
1ed 


as 


om 


H* 
the 
ick 
“on 
an 
pat 
as 


his 





PRODUCTS OF TONIZATION IN GASES 





























/] 
+— 
HZp At 
ra 
> | 
=} | —+ oe 
| | 
c| 
>| 
=| | | 
1 
| | [ 
| 4 / 
r) p 
7 9 12 4 | 15 19 23 ef 


Electron enerey (volts) uncorrected 
Sy 


F1G. 6. Typical ionization potential curves in acetylene. 
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Fic. 7. Heights of ion peaks as a function of the electron energy in acetylene. 
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ion in nitric oxide. Electron energy V,=100 EV. 
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Fic. 9. Heights of the ion peaks as a function of the electron energy in nitric oxide. 


effect will be the more marked the smaller the 
mass of the ion relative to the molecule. The C* 
and CH?* ions have as much as 1.5 EV kinetic 
energy while the H* ion has up to 6 EV. Fig. 7 
shows the heights of the various peaks (probabil- 
ity of formation) as functions of the electron 
energy in volts. No negative ions were observed. 

There are two peaks labeled C,H2*. The one 
on the left (m/e= 27) is presumably (C®C!H,)*. 
This ion appears at the same electron velocity as 
the ion of mass 26 and the relative heights of the 


two peaks remain the same as the electron veloc- 
ity is varied. The relative abundance of C'* to C® 
is approximately 1 part in 100. Jenkins and Orn- 
stein® have reported a value of 1 part in 106. A 
more recent and refined determination® gives 1 
part in 91.6. Aston’ has found a ratio of 1 part 
in 140. 


5F. A. Jenkins and S. S. Ornstein, Proc. Acad. Sci. 
Amsterdam 35, 1212 (1932). 

6A. L. Vaughan, J. H. Williams and J. T. Tate, Phys. 
Rev. 46, 327 (1934). 

7F. W. Aston, Proc. Roy. Soc. A149, 396 (1935). 
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B. Nitric oxide 


The sample was prepared by the action of 
HNO; on Cu. The NO was passed through water 
to remove any NOs, dried over solid NaOH, col- 
lected over liquid air, and purified by repeated 
fractional distillation. After all these precautions 
there still remained approximately one percent of 
N.O as determined by the mass spectrograph. 
Fig. 8 shows the mass spectrum. Fig. 9 is the 
height of the peaks plotted as a function of the 
electron energy in volts. The negative ions 
found by Tate and Smith® showed up strongly 
and were identified as O~ ions. The results are 
given in the following table. 


lon NO* * NO** 
App. pot. (volts) 9.5 22.0 44 
Error (volts) +0.1 +0.5 +0.1 


A very faint trace of O* ions appeared at 21 volts. 
The ionization potential of NO is taken as 9.5 
volts from the work of Tate and Smith® and this 
value was used to calibrate the scale. 

In Fig. 8 there is an ion having an m/e value 
of 31. This peak seems to bear a constant relation 
to the NO* peak of about 1/2 to 1/3 percent 
independent of electron velocity and is probably 
due to the N'* isotope. A more careful study® of 
this has given a ratio of N'4, N! of 265. 


C. Cyanogen 


The sample of C:Ne was obtained from the 
Department of Chemistry. The results may be 
compared to those of Dorsch and Kallman.’ 





8]. T. Tate and P. T. Smith, Phys. Rev. 39, 270 (1932). 
® kK. E. Dorsch and H. Kallman, Zeits. f. Physik 60, 376 
(1930). 
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F1G. 10. Heights of the ion peaks asa function of the electron 
energy in cyanogen. 


lon : C.N.* CNt C,* Ct C.N* 
App. pot. ) Authors 14.1 21.3 18.6 19.8 
(volts) » Dorsch and 
Kallman 13.5 18 17 22.5 


Fig. 10 shows the heights of the peaks as a func- 
tion of the electron velocity in volts. C+ and N* 
peaks were found, but their origin was uncertain 
since the sample of cyanogen contained small 
amounts of CO and N, as impurities. 


D. Carbon monoxide 

The work on carbon monoxide agrees with that 
published by Vaughan.'® The O* ion which was 
masked by H2O* in that work appears here but is 
very small. It appears at 27+1 volts and is only 
about 10 percent as plentiful as the C* ion at 100 
volts electron velocity. 

The authors are very grateful to C. H. Dane, 
E. F. Greinke and William B. Haliday for their 
help and cooperation. The senior author wishes 
to express his appreciation to the Graduate 
School of the University for a grant from the 
Fluid Research Funds. 


1” A. L. Vaughan, Phys. Rev. 38, 1687 (1931). 
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Effect of a Magnetic Field Upon the Polarization of Impact Radiation 


McCorke,* Physical Laboratory, University of Iowa 
(Received July 8, 1935) 


W. H. 


Sodium PD line radiation excited by electrons moving parallel to a magnetic field of 300 gauss 
or more shows about 16 percent polarization, the major electric vibration being parallel to the 
direction of the exciting electron beam. In fields below 300 gauss the polarization is less, de- 
creasing to zero in the absence of a field. That the polarization changes in a magnetic field 
must be due to the fact that the 3°P3/2 level isa hyperfine multiplet which undergoes a Paschen- 





Back effect. 


INTRODUCTION 


N most cases the application of the principle 

of conservation of angular momentum to the 
excitation of radiation by electron impact leads 
to results which are at least in qualitative agree- 
ment with experiment. The failure of sodium 
D line impact radiation to show polarization! * 
might be due to exchange, operating in such a 
way that the upper magnetic levels are always 
equally populated. However, it may be doubted 
that the probability of exchange should vary 
with velocity in such a way that this should be 
true over a considerable velocity range. Another 
possible explanation of the low 
that it results from the hyperfine structure of the 
levels involved. It is known that this is responsi- 
ble for greatly lowering* the polarization of 
resonance radiation. Any effect of the hyperfine 
structure should however, disappear for the im- 
pact radiation just as for the resonance radiation‘ 
in a magnetic field strong enough to produce a 
complete Paschen-Back effect of the upper 
hyperfine multiplet. 

With this in mind the writer has looked for 
polarization in the D line impact radiation 
excited in a magnetic field parallel to the direc- 
tion of motion of the exciting electrons. 


polarization is 


APPARATUS AND PROCEDURE 


The impact radiation was excited in a tube 2 
cm in diameter which had on one side a window 
blown and flattened and on the opposite side a 


* Now at Agricultural and Mechanical College, College 
Station, Texas. 

' Kossel and Gerthsen, Ann. d. Physik 77, 273 (1925). 

* Ellett, Foote and Mohler, Phys. Rev. 27, 31 (1926). 

* A. Ellett, Phys. Rev. 35, 588 (1930). 

'L. Larrick, Phys. Rev. 46, 581 (1934). 


light trap in the form of a horn. A cylinder of 
sheet molybdenum with two rectangular aper- 
tures was inserted in the tube so that the aper- 
tures came opposite the window and light trap, 
respectively. The molybdenum cylinder con- 
tained the electron gun, which consisted of an 
oxide-coated platinum filament placed about 2 
mm from a pair of 3-mm diaphragms arranged to 
project the electron beam into the space between 
the window and the light trap. The molybdenum 
cylinder and diaphragms were placed at the 
same electrical potential with the idea of provid- 
ing a field-free space into which the electron 
beam should be projected. However at low elec- 
tron velocities the potential difference due to 
space charge became troublesome, and in order 
to compensate this to some extent a third inde- 
pendent electrode was introduced, see Fig. 1. 
The entire electron gun except the filament was 
constructed of sheet molybdenum and arranged 
as shown in Fig. 1. 

The oxide coating of strontium, barium and 
calcium oxides was applied by the method 










































































(a) | (6) 


Fic. 1. Tube used for the excitation of impact radiation. 
(a) Sectional side view of tube. (b) Front view of tube. 
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POLARIZATION OF 


described by Arnold.’ The filament used would 
yield a maximum current of 1 to 1.5 milliamperes 
in the electron beam with the accelerating 
potentials employed. 

The heating chamber which was an asbestos 
cylinder painted on the inside with lampblack 
and shellac having a rectangular opening in one 
side to permit observation of the radiation ex- 
cited by the electron impacts, extended to a gas 
burner 3 feet below. A water-cooled Helmholtz 
coil of copper tubing placed with its axis vertical 
and coincident with the direction given by a 
line drawn through the center of the filament and 
the defining diaphragms surrounded the heating 
chamber and provided the variable magnetic 
field to which the impact radiation was subjected. 
The entire assembly was mounted within a large 
Helmholtz coil adjusted to neutralize the earth's 
magnetic field. 

The sodium which had been distilled in a 
vacuum and collected in small glass tubes was 
placed in one of the side tubes of the observation 
chamber and when the system had been suffi- 
ciently evacuated was redistilled into the obser- 
vation chamber until a light coat of pure metallic 
sodium became deposited on the top and sides 
of the main tube. 

A “sticking’’ vacuum, as indicated by the 
McLeod gauge, was maintained by means of a 
Cenco Hyvac and a mercury diffusion pump. 

With the heating chamber adjusted for any 
desired temperature, observations were made of 
the impact radiation. The observations were 
made perpendicular to the electron beam and 
magnetic field, first, visually by means of a 
quartz wedge® and Ahrens prism, and later, 
photographically by the same method, and also 
photographically by the Cornu method employ- 
ing a Wollaston double image prism and an 
Ahrens prism. The position of the Ahrens prism 
which gave equal intensity of the two images was 
determined by microphotometer measurements 
of the two images obtained for various settings 
of the prism. The angular position of the Ahrens 
prism was read from a vernier scale arranged to 
read to tenths of a degree of arc. 

“TH. DD. Arnold, Phys. Rev. 16, 70 (1920), 


"R. W. Wood, Physical Optics, third edition, pp. 339 
340. 
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RESULTS 


From the work of Larrick‘ it is known that a 
magnetic field of 300 gauss is sufficient to com- 
plete the Paschen-Back effect of the hyperfine 
structure for the sodium 3°P3,2 hyperfine multi- 
plet and produce the maximum polarization of 
resonance radiation. Hence observations of the 
impact radiation of sodium were first made in a 
magnetic field of 300 gauss. Preliminary observa- 
tions showed that below sodium vapor pressures 
of 210-7 mm of Hg, absorption and secondary 
resonance radiation did not influence the polar- 
ization of the impact radiation. 

When the electrons in the incident beam had 
very nearly the critical velocity corresponding 
to the excitation potential the polarization of the 
impact radiation was found to increase with the 
applied magnetic field to a value slightly above 
16 percent, with the maximum polarization ap- 
proximately realized at 200 gauss and remaining 
nearly constant for higher fields as Fig. 2 shows. 
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ic. 2. Polarization vs. magnetic field for electrons having 
energy within one volt of the resonance potential. 
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ic. 3. Polarization es. magnetic field for electrons with 
energy 3-5 volts above the ionization potential. 
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As one might expect the speed of the incident 
electron was an important factor in determining 
the magnitude of the polarization. The speed of 
the incident electron was controlled by the 
applied accelerating potential and the space 
charge effect in the electron beam. For accelerat- 
ing potentials in excess of the ionization potential 
the space charge disappeared and there was a 
decrease of the polarization for all magnetic 
fields as indicated by Fig. 3. If the accelerating 
potential and electron current were fixed to give 
the space charge effect for higher fields but not 
adjusted for the lower fields there was a decrease 
of the value of the lower field polarization caused 
by the change of space charge with magnetic 
field. This is shown by comparison of Fig. 4 and 
ig. 2. 

Ata temperature of 130°C the impact radiation 
was apparently limited to a well-defined beam, 
while at higher temperatures a nebulous glow 
of secondary radiation could be seen to surround 
this beam. Observations made by means of the 
quartz wedge and Ahrens prism above 130°C 
showed polarization of the secondary radiation 
as may be seen in Fig. 5. Observations of the 
main beam by this method failed to show definite 
evidence of polarization either in the field of 300 
gauss or in the absence of a magnetic field. 
However, when the Cornu method was employed 
at a somewhat lower temperature, polarization 
of the impact radiation was clearly indicated. 
The degree of polarization was still very small 
even in the field of 300 gauss. It was observed, 
however, that when large electron currents were 
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hig. 4. Polarization vy. magnetic field showing decrease 
of low tield polarization due to modification of the space 
charge by the magnetic field. 
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a b c 


Fic. 5. Photograph of impact radiation using a quartz 
wedge and Ahrens prism. (a) Fringe system with slightly 
convergent light incident on the quartz wedge, with the 
sodium vapor at 160°C in a magnetic field of 400 gauss. 
(b) Fringe system in secondary radiation when the image 
of the main beam coincides with the plane of the wedge 
and the sodium vapor is at 150°C in a magnetic field of 
400 gauss. (c) Same as (b) above with the sodium vapor 
at 130°C, 





a b 


Fic. 6. Beam of radiation excited by electron impact. 
(a) Effect of space charge in the electron beam. (b) Space 
charge removed by higher accelerating potential. 


used there occurred a break in the beam of im- 
pact radiation (see Fig. 6) and at the same time 
the polarization increased. 

For magnetic fields of various intensities with 
the sodium vapor at 110°C it was found that the 
polarization in any given field was a maximum 
when the space charge in the electron beam was 
sufficient to render a portion of the beam non- 
luminous. The radiation excited by the electrons 
on either side of the nonluminous portion of the 
beam is attributed to electrons with energy close 
to the resonance potential, since in the non- 
luminous part of the beam their energy must be 
insufficient to produce excitation. To verify this 
the polarization of the impact radiation was 
determined in magnetic fields of 400 and 500 
gauss for accelerating potentials less than one 
volt above the resonance potential. The same 
strong field polarization was observed under these 
conditions. The exposure time required for a 
picture used in the microphotometer measure- 
ments was however increased from a few minutes 
to 8 or 9 hours with this adjustment of accelerat- 
ing potential and necessarily smaller electron 
current, consequently the observations with the 
low accelerating potential were not applied to 
other points of the curves. 
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lic. 7. Polarization vs. vapor pressure for sodium vapor in 
a magnetic field of 300 gauss. 


From consideration of the principle of con- 
servation of angular momentum the polarization 
should be expected to be greatest when the 
energy of the exciting electron is equal to or just 
exceeds that necessary for excitation. This is 
confirmed by the above observations of increased 
polarization in the region where the electron 
beam was retarded by space charge. The effect of 
the space charge on the electron beam was de- 
creased by increasing the vapor pressure of the 
sodium since more ionizing collisions occurred. 
In addition to the decrease of space charge with 
increase of vapor pressure there was an increase 
of the depolarizing effect of the secondary radia- 
tion. The result of these combined effects 
caused the polarization to decrease from about 
16 percent at 105—110°C to nearly zero at 130°C 
as shown by Fig. 7. 


DISCUSSION 


It is evident that measurements of the polar- 
ization of the radiation excited by electron im- 
pact on sodium vapor in a magnetic field must be 
made below 130°C to escape the effect of the 
polarized secondary radiation. There is a de- 
polarizing effect due to the secondary radiation 
appreciable at temperatures as low as 115°C but 
the effect is exaggerated by Fig. 7 which must be 
to a considerable extent characteristic of the 
apparatus, since space charge is affected by the 
vapor pressure of the sodium due to increase in 
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IMPACT RADIATION 


the number of positive ions with increase of 
vapor pressure. 

There is qualitative agreement between the 
polarization of impact radiation from sodium 
observed perpendicular to the applied magnetic 
field and electron stream, and the polarization 
of the resonance radiation excited by the sodium 
D lines polarized with E //. The separation of 
the hyperfine levels appears to be completed at 
about the same value of the magnetic field for the 
impact radiation as for the resonance radiation 
with E /I but the value of the polarization 
reaches a maximum of only 16 percent for the 
impact radiation compared with about 50 per- 
cent for the resonance radiation. 

If the electron beam consisted only of electrons 
having energy equal to the resonance potential 
and moving accurately parallel to the magnetic 
field, it would seem upon the basis of a simple 
theory that no angular momentum parallel to 
the magnetic field could be imparted to the ex- 
cited sodium atom. The excited atom should 
then be raised to the same state as when ex- 
cited by light polarized with its electric vector 
parallel to the magnetic field. That is, in strong 
fields only the 17 ;= +3 levels of the 3°P3,. state 
should be populated, and the light emitted by 
the transitions from the 3°P3,2 should be 60 per- 
cent polarized. The polarization observed should 
be less than this because of unpolarized emission 
from 3*P,, but probably not less than 33 percent, 
corresponding to equal excitation of the /;= +} 
levels of ?P 3 and *P,. 

The low (16 percent) strong field polarization 
observed shows that the 1/;=+3/2 levels of 
*P3 are excited, and this requires the com- 
munication of angular momentum parallel to the 
magnetic field. The writer was not able to work 
with accelerating potentials less than one-half 
volt above the excitation potential, for as the 
excitation potential is approached more closely 
the intensity drops to so low a value that polariza- 
tion measurements are not feasible. 

In conclusion the author wishes to take this 
opportunity to thank Professor A. Ellett for his 
assistance and suggestions during the course of 
this research. 
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Atomic Energy States for Excited Helium 


Wa. S. Witson, Sterling Chemistry Laboratory, Yale University 
(Received July 8, 1935) 


By means of Slater’s method the total atomic energies for the excited states, 1s?, 1s2s, 1s2p, 
2s*, 2s2p und 2p? of helium have been calculated with the use of Hartree functions previously 


published. 





INTRODUCTION 
pecs method! of calculating total atomic 


energies has been applied to several atoms.’ 
The success of these calculations suggested the 
desirability of determining theoretically the 
more important term values of the unobserved 
states of doubly excited helium, by using the 
results of self-consistent field calculations pre- 
viously published.* The results are compared 
with experimental data as far as possible. 


METHOD OF CALCULATIONS 


The method used in the present calculations 
was to consider the two-electron problem as a 
special case of the n-electron problem described 
by Slater, the procedure being similar to that 
used by Hartree and Black. Since no essentially 
new features were developed for this work, any 
description of the method will be omitted because 
the method is amply described in the two refer- 
ences given. However, there is one point to be 
noted; that is, if we omit the Hartree energy 
parameter from the definition of J(@) as given by 
Hartree and Black, we can at once write the 
equality: 


I(a) = 1(8) = Fo(aB) = Fo(Ba) 


for those states of two electron atoms or ions 
for which the Hartree functions are already 
orthogonal. In the present work this equality 
was used to check the numerical work by calcu- 
lating a given quantity in two _ essentially 
different ways. Should the Hartree functions be 
changed by orthogonalization, it is necessary to 
use the method of treatment given by Hartree 
and Black. This was here necessary only for the 
(1s) (2s) 4S level. 

In the actual calculations the full five place 
figures used in the Hartree calculations as well 

1]. C. Slater, Phys. Rev. 34, 1293 (1929). 

2See, for example, D. R. Hartree and M. M. Black, 
Proc. Roy. Soc. A139, 311 (1933); J. Macdougall, ibid. 
A138, 550 (1932). 


*W.S. Wilson and R. B. Lindsay, Phys. Rev. 47, 681 
(1935). 


as the values at the intervals omitted in the 
tabulation of references were used. These extra 
intervals enabled all integrals to be evaluated 
conveniently by direct quadratures rather than 
with the more elaborate method described by 
Macdougall for large k. The energy values were 
obtained from the integrals by using Eq. (8) of 
Hartree and Black in the manner outlined by 
Slater in the latter part of his article. 


RESULTS AND DIscUSSION 

The results of these calculations are presented 
in Tables I and II. Table I gives the mean values 
of the various Slater integrals, Table II the 
energy values obtained from these, together with 
experimental values and the results of other com- 
putations where such are available. The fourth 
column of this table gives the values obtained if 
the smaller Hartree energy parameter is assumed 
to be equal to the first ionization potential of that 
state. The agreement of the calculated energy 
values with the experimental ones where such 
are available is satisfactory, there being a maxi- 
mum deviation of 1.5 percent for (1s)? '\S and a 
much smaller deviation for other levels. The 
agreement with calculations is very good, save 
for 'S states. An exception to this last statement 
is the discrepancy of the (2p)? °P level when com- 


TABLE I. The Slater functions for the several excited states 
(in units Rhe where R= 109,737). 











CONFIGU-  ELEC- I(n,}) =1(n'l’) 
RATION TRON € =Fo(nl, n’l’) G(nl. nl’) 
(1s)? Is 1.836 2.051 
(2s)? 2s 0.4606 0.5189 
(1s) (2s)’S 1s 3.469 0.5317 Go =0.0740 
2s 0.3068 
(1s)(2s)3S Is 3.469 T(1s) =0.5317 Go= .0123 
I(2s) =0.4911 
2s 0.3068 Fo(1s2s) =0.4764 
(1s) (2p) 1s 3.496 0.5070 Gi= .0261 
2p 0.2522 
(2s)(2p) 2s 0.4195 0.5403 Gi= .2942 
2p 0.4443 
(2p)? 2p 0.3861 0.5761 Same as F;(2p:, 2p) 
F2(2p, 2p) =0.2675 
Bett(1s)? 1s 11.344 4.558 
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TABLE IL. The atomic energies (units Rhe). 





TABLE III. Observed lines from the helium spectrum. 








Conricgu- = MutL- EXPERI- OTHER 
RATION TIPLET ENERGY HIARTREE MENTAL! METHODS 
- (Us)? "Ss —5.723 —5.836 —5.807 —5.8092 
(2s)? Ss — 1.4400 — 1.4605 — 1.494 
(1s)(2s) ‘Ss —4.232 —4.3062 —4.2919 —4.2902 
—4.2893 
as —4.334 —4.3499 — 4.3345 
(1s) (2p) - —4.246* —4.2517 —4.2473 — 4.2455 
3P —4.264 — 4.2600 —4.2625 
(2s)(2p) -_ — 1.3060 — 1.4194 — 1.3014 
— 1.3086 
3P — 1.5022 — 1.4984 
— 1.5046 
(2p)? 3p — 1.4018 — 1.3860 — 1.6626 
’D — 1.3376 — 1.3206 

"Ss — 1.2323 

Be**{ is}? ‘S — 27.246 —27.344 —27.312 — 27.3087 





1 Bacher and Goudsmit, Alomic Energy States. 

2 E. A. Hylleraas, Zeits. f. Physik 65, 209 (1930). 

3 J. P. Vinti, Phys. Rev. 37, 448 (1931). 

4F. G. Fender and J. P. Vinti, Phys. Rev. 46, 78 (1934). 
§C. Eckart, Phys. Rev. 36, 878 (1930). 

6 Ta-You Wu, Phys. Rev. 46, 239 (1934). 

7E. A. Hylleraas, Zeits. f. Physik 83, 739 (1933). 


pared with Wu's value. The value listed by him 
gives too large a multiplet separation.‘ Further 
as we will see later the present value is in line 
with some experimental predictions. The simple 
Hartree values in most cases gave very closely 
the center of gravity of the multiplet system as 
calculated by the Slater method for states having 
but two levels. The values for Be*+*(1s)? 1S° are 
included to show the trend for the more closely 
bound beryllium charge. The fractional error is 
smaller than for normal helium. Whether this is 
generally true for other levels of Be**+ cannot be 
stated at present. 

The assignment of the newly calculated levels 
to experimental results is restricted because of 
the lack of data. If we calculate transitions from 
singly to doubly excited states using the known 
experimental values for the singly excited states, 
we may note the following results. Wilson and 
Lindsay’s assignment of the transition (1s)? 1S 
—(2s)?'S to the energy lost observed by 
Whiddington and Priestley seems to be further 
verified. It might be that the line is due to 
transitions from 1s?'S to both (2s)?'S and 


4 If we apply the rule that the *P —'D separation is three 
halves the 'D—'S separation we find, using Wu’s value for 
this separation, the 4S level to lie at —0.777 Rhe, which 
indeed is unlikely since this is greater than the ionization 
potential of the ion He*(2p). 

5 Based on the wave function recently published by 
D. R. Hartree, Proc. Roy. Soc. A149, 210 (1935). 





LINE PREVIOUS ASSIGNMENT PRESENT PAPER 
357.507! (1s)(28)’S—(2s)? “1 “none 
322.517! (1s)(2p)"P —(2 p)?’p none 
321.186! (1s)(2p)’P —(2p)?2’S! none 
320.392! (1s)(2 s)8S —(2p)?8P1 (1s)(2s)8S —(2s)(2p)3P 


(1s)(2 s)8S —(2s)(2p)8P3 
309.042 (1s)(2p)’P —(2s)(2p)'’P8 


1P. G. Kruger, Phys. Rev. 36, 855 (1930). 
2K. T. Compton and J. C. Boyce, J. Frank. Inst. 205, 497 (1928). 
3 F. G. Fender and J. P. Vinti, Phys. Rev. 46, 78 (1934). 


(1s) (2s)8S —(2p)? 3P 











(2p)? *P which lie very close together. This would 
also explain the indications of structure observed. 

The results of a study of the observed lines in 
the helium spectrum and their previous assign- 
ments are given in Table III. In this respect the 
results of these calculations verify those given 
by Fender and Vinti with the one exception 
noted. Calculation of transitions between the 
doubly excited levels fail to fit any of the corona 
lines discussed by Rosenthal. The values thus 
calculated might, however, be in error by an 
amount which would mask any agreement and 
too definite conclusions cannot be drawn from 
these transitions between doubly excited states. 

It is believed that the error in the doubly 
excited states is much less than that in some of 
the singly excited states. This belief is founded 
on the excellent agreement of calculations for 
higher atoms by other workers with observed 
values which involve 2s and 2 electrons, as 
well as the agreement with other calculations 
noted in Table II. Further, should it prove that 
the absolute errors in these levels remain essen- 
tially the same for all levels of higher two elec- 
tron ions as it does for Be** (1s)? 4S, this would 
introduce a great simplification in the calculation 
of the relative positions of many closely lying 
transitions for the helium atom. For in that case 
one need only compute the larger and more 
easily obtainable similar separations for higher 
ions and contract them suitably, thus rendering 
the error small. The Hartree-Slater method, 
though having the great drawback of involving 
considerable work seems on the whole to give 
quite tenable results for two electron atoms and 
ions. 

It is a pleasure to thank Professor Henry 
Margenau of the Sloane Physics Laboratory for 
his many helpful discussions and suggestions. 


6 A. H. Rosenthal, Zeits. f. Astrophysik 1, 115 (1930). 
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Dielectric Constant of Water Vapor 


J. D. SrRANATHAN, Department of Physics, University of Kansas 


(Received June 27, 1935) 


The frequently discussed deviations from linearity in the 
curves of (A—1)/(A+2) versus pressure for water vapor 
have been investigated. It has been shown: (1) That with 
proper location of the insulators of the vapor condenser 
in weak fields, the dielectric constant curves for the vapor 
are jinear until saturation is approached closely; they are 
linear over the entire range of pressures at which breaks 
have previously been reported; (2) that neither pronounced 
breaks similar to those reported by several workers, nor 
the small residual deviations always present near satura- 


On THE DEVIATIONS FROM LINEARITY 


ARLY dielectric data for water vapor ob- 

tained by Jona! did not conform to the 
theoretical Debye line at the lower temperatures; 
the deviations were attributed to association of 
vapor molecules. Zahn* obtained data which 
showed the measured polarization of water vapor 
not at all proportional to the pressure at tempera- 
tures of 40°C or less. Pronounced breaks occurred 
at a pressure of the order of 1 cm of Hg; the 
lower the temperature the more pronounced was 
the break. Zahn attributed these breaks to the 
probable existence of a film of adsorbed water on 
the surfaces of the condenser plates. He pointed 
out that if observations were confined to pres- 
sures below those at which the breaks occurred, 
the data, while not accurate, were not incon- 
sistent with the Debye theory. Wolf* confirmed 
the general nature of the breaks reported by 
Zahn, but observed that they were entirely 
independent of plate spacing. He reported also 
that the magnitude of the break depended upon 
the insulating material used in the vapor con- 
denser, being least for quartz; it appeared that 
even the slope of the line below the break was 
somewhat dependent upon the kind of insulator 
used. He further observed that the d.c. con- 
ductivity of the vapor condenser varied with the 
type of insulator used and with the pressure of 
the vapor in the condenser. The fact that the 
conductivity increased rapidly at approximately 
those pressures for which the breaks in the 

! Jona, Physik. Zeits. 20, 14 (1919). 


2 Zahn, Phys. Rev. 27, 329 (1926). 
* Wolf, Ann. d. Physik 83, 884 (1927). 


tion are due to the conductivity effect as commonly sup- 
posed; (3) that these breaks are due wholly to the added 
polarization contributed by the film of water molecules 
adsorbed on the insulator surfaces. Accurate dielectric data 
have been obtained for water vapor at 14 temperatures 
ranging from 21.3°C to 197.9°C; they can be represented 
by the equation [(A—1)/(A +2) ]RT7/p’= (4.03 +0.39) 
+ (20,710+140)/7T. These data fall accurately along a 
Debve line, vielding an electric moment of (1.831 +0.006) 
x 1078 e.s.u. for the water molecule. 


dielectric curves occurred for the several insu- 
lators, led him to attribute the existence of the 
breaks to the conductivity effect. Zahn‘ ap- 
parently confirmed these observations and was 
led to the same conclusion. Similar breaks, 
though occurring much nearer saturation, have 
since been observed for many materials. Maske? 
reported them for benzophenone vapor; Knowles® 
reported them for ethyl alcohol; and the author 
has found though not reported them for all 
vapors studied near saturation, including numer- 
ous alcohols, toluene and the nonpolar vapor 
benzene. Knowles studied variations in d.c. 
conductivity in the region of the breaks for ethy] 
alcohol vapor, and concluded that both experi- 
mentally and theoretically the breaks could not 
be caused by a conductivity effect; the con- 
ductivity effect of the proper 
magnitude nor, for his experimental arrangement, 


was neither 
in the proper direction. 

Since there appears some evidence conflicting 
with each of the three suggested interpretations 
of the breaks (association, the added polarization 
due to adsorption on the condenser plates, and 
conductivity across the insulator surface), the 
author has undertaken to get at the real cause. 
He has chosen to study water vapor, that 
material yielding the most pronounced breaks. 
By use of an unusually stable heterodyne beat 
apparatus’ operating entirely on alternating 
current, along with a carefully calibrated variable 

* Zahn, Phys. Rev. 35, 1047 (1930). 

5 Maske, Physik. Zeits. 28, 533 (1927). 


® Knowles, J. Phys. Chem. 36, 2554 (1932). 
7 Stranathan, Rev. Sci. Inst. 5, 334 (1934). 


538 





p- 
led 
les 
ata 
res 
red 
39) 


16) 


ng 
ns 
Qn 
id 
ne 


at 
S. 
at 
ig 
le 





DIELECTRIC CONSTANT OF WATER VAPOR 539 





Water Vapor — 23.3" c 










os en ° |e Quarfe “aes 


. + Pyrex Insulators 






| 
} 
| 





eZahns Data 
120| $———-.——_ 
525000 Ohins ” areMe 


me —vadde d Pyrex Sutface| 







= 
} 
| 































Fic. 


capacity standard,’ both previously described by 
the author, capacity changes as small as a few 
thousandths of a wuf due to the vapor could be 
observed accurately. It was found’ immediately 
that by using a vapor condenser similar to that 
described by Knowles, in which care is taken to 
locate the insulators in weak fields, the dielectric 
polarization of water vapor varies linearly with 
pressure up to pressures far higher than those at 
which pronounced breaks have previously been 
reported. Large deviations could be observed 
only if the vapor were allowed to enter the 
condenser too rapidly; these could not be re- 
peated and were probably due to local conden- 
sation. The addition of a tightly packed glass 
wool plug through which the vapor must pass in 
going from the vapor system (which was that 
used by Knowles") to the condenser, completely 
eliminated any such possibility of error. Curve 1, 
Fig. 1 shows the results at 23.3°, using quartz 
insulators. Plotted in the same figure, curve 4, is 
the pronounced break reported by Zahn? for the 
same temperature, and confirmed approximately 
by Wolf.’ It is apparent that the present curve 
(1) is essentially linear up to the highest pressures 
for which observations were heretofore reported. 

Curve 2 represents data obtained using Pyrex 


*Stranathan, Rev. Sci. Inst. 5, 315 (1934). 
® Reported at the Washington meeting of the Am. Phys. 
Soc., 1934. 





insulators. The break from linearity is somewhat 
larger and occurs at a somewhat lower pressure 
than for quartz, though it is even here not 
comparable with previously reported breaks. 
The slope of the line below the break is inde- 
pendent of the kind of insulator used. Fig. 2 
shows similar data taken at 29.6°, curve 2 being 
the author’s data obtained with Pyrex insulators, 
and curve 4 observations reported by Zahn. 
Even using Pyrex insulators, the curve is accu- 
rately linear up to pressures higher than the 
highest for which observations have been re- 
ported previously. While data using quartz 
insulators were not actually taken at this temper- 
ature, they have been taken for various nearby 
temperatures, and the results leave no doubt but 
that such a curve would be linear up to a still 
higher pressure, and the extent of the break less. 

In spite of the fact that several observers who 
have measured the conductivity of the vapor 
condenser have attributed the nonlinearity of the 
dielectric curves to the conductivity effect, it has 
always appeared impossible to the author to 
reconcile the magnitude of the break with the 
measured conductivities. The author has there- 
fore made three tests to show that the con- 
ductivity is not the cause of even the small 
residual breaks found here. First, since all 
previous workers seem to have measured d.c. 
conductivities, and since it is now recognized 
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from the works of Yeager and Morgan’? and of 
Roberds" that the high frequency conductivity is 
larger than the d.c. conductivity, it seemed 
advisable to measure the conductivity of the 
vapor condenser at the same frequency as was 
used in making dielectric constant measurements. 
The author has made such measurements (at 545 
kilocycles), by a method since described and used 
extensively by Roberds." Since the large con- 
denser (840 uuf) used in dielectric measurements 
is not conducive to accurate measurement of high 
leakage resistances at this frequency, no claim is 
made of accuracy other than as to order of 
magnitude. It can be stated definitely that the 
high frequency leakage resistance of the con- 
denser was never less than 5 megohms at any 
vapor pressure of water less than. saturation; 
except at pressures barely below saturation, it 
was very much higher. This value agrees with 
calculations from the specific conductivities 
quoted by Yeager and Morgan and by Roberds. 
Theoretically it appears that the magnitude of 
the break is entirely inconsistent with this high a 
leakage resistance. 


10 Yeager and Morgan, J. Phys. Chem. 35, 2026 (1931). 
1! Roberds, Physics 6, 227 (1935). 


Second, the author has made several experi- 
mental measurements of the effect of conduc- 
tivity, measurements similar to those reported by 
Knowles.*® While it is impossible to introduce an 
artificial leakage resistance in parallel with the 
test condenser without introducing at the same 
time an appreciable capacity, it is possible by 
choice of appropiate resistors to introduce various 
leakage resistances all with essentially the same 
capacity; thence the effect of conductivity can be 
deduced. This has been done, using artificial leaks 
varying from 10° down to 10+ ohms. Groups of 
such artificial leaks of several different types 
have been used, and the results obtained have 
been consistent as regards the magnitude of the 
conductivity effect. Experimentally the effect of 
conductivity is here in the proper direction but of 
entirely too small a magnitude to cause even the 
small residual breaks observed by the author. 
The leakage of the vapor condenser would 
apparently have to be at least 200 times as great 
as the actually measured value in order that the 
conductivity might produce breaks comparable 
with those observed. 

Third, a test was made by placing an artificial 
leak directly in parallel with the vapor condenser 
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and leaving it there throughout the determi- 
nation of the dielectric constant curve. Thus the 
vapor condenser was always of quite high leakage 
regardless of the pressure of the vapor within. 
It can easily be shown that, under such con- 
ditions, any change in leakage accompanying 
increases in vapor pressure would produce an 
entirely different effect upon the measured 
capacity than would be the case were the low 
resistance leak not permanently in parallel. As 
shown in both Figs. 1 and 2, dielectric observa- 
tions, including the residual breaks near satu- 
ration, were identical with and without a 
permanent parallel leak of 25,000 ohms. Thus all 
three pieces of evidence here cited seem to show 
definitely that even the small residual breaks in 
the dielectric curves just below saturation cannot 
be caused by the conductivity effect. 

In the test cited just above, one in effect 
changes the conductivity of the vapor condenser 
without changing the nature or extent of the 
insulator surfaces on which water molecules 
might be adsorbed. A test was next made under 
such conditions as would change markedly the 
extent of the adsorbing surface within the con- 
denser without changing the leakage resistance. 
This was done in the following way: In the 
condenser used (see work of Knowles® for dia- 
gram), the two sets of plates were built as two 
rigid metallic units. The one unit is then sup- 
ported by four cylindrical compression insulators, 
1.2 cm diameter, and 1.6 cm long. Three of these 
support the insulated unit on the bottom of the 
condenser housing, while the fourth is placed on 
top of the unit and has applied to it the com- 
pressional force which clamps the unit into place. 
While the normal plate spacing is approximately 
1.5 mm, the insulators are all mounted between 
plates whose spacing is at least 10 times as 
large, thus placing them in a relatively weak 
field. Into the upper side of the lower of the two 
top plates (those between which the upper 
insulator was located) there were drilled 47 
holes passing only part way through the plate. 
Into each of these was wedged in an upright 
position a piece of 8 mm Pyrex tubing carefully 
cut to such a length that it failed by approxi- 
mately 2 mm to touch the adjacent plate. In this 
manner the effective insulator surface area which 
may adsorb water was increased by a factor 


approximately 6 but, since the new insulator 
material touches only one plate, the leakage 
resistance of the condenser remains unchanged. 
Dielectric data” taken under these conditions 
are shown by curves 3 of Figs. 1 and 2. Again the 
same curve is obtained with or without the use of 
a permanent parallel leak of 25,000 ohms. The 
presence of additional adsorbing insulator surface 
thus emphasizes the break; the break is of much 
greater magnitude, the ratio of magnitudes being 
roughly the ratio of effective adsorbing surfaces, 
and it occurs in measurable magnitude at a 
lower pressure. Dielectric curves obtained under 
these conditions show breaks quite similar to 
those reported by Zahn and Wolf. One can 
scarcely escape the conclusion that all breaks 
obtained here are caused almost wholly by the 
added polarization contributed by the presence 
of a water film adsorbed on the surfaces of the 
insulators. While one cannot say definitely that 
previously reported breaks have not been due to 
conductivity, for it is always possible that an 
oscillator whose frequency is unusually dependent 
upon leakage has been used, still it appears 
unlikely indeed. It is far more probable that all 
such breaks have been caused by the added 
polarization contributed by the water film 
adsorbed on the insulators. 

Two further experiments have been carried out 
to see whether the presence of such an adsorbed 
film of reasonable thickness would produce an 
added capacity of the proper magnitude to 
account for the observed break. An effort was 
first made to calculate by rough approximation 
the capacity contributed by a film of given 
thickness, taking into account the fringing of the 
electric flux caused by the insulator and by the 
long thin adsorbed film of dielectric constant 80. 
Actual measurements of the capacity increases 
caused by similar though much thicker films (a 
few tenths of a mm thick) confined between two 
closely fitting pieces of glass tubing placed in a 
position between plates similar to that of the 
insulator, though touching only one of the plates, 
showed that the crude method of calculating the 
capacity contributed by the thin film gave 
surprisingly the right magnitude. Such experi- 
ments and calculations indicate that an adsorbed 


2 Reported at the Washington meeting of the Am. Phys. 
Soc., 1935. J. D. Stranathan, Phys. Rev. 47, 794A (1935). 
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Fic. 3. For the lower temperature curves both ordinate and abscissa have been multiplied by 
the factor indicated. 


water film a few hundred molecules thick would 
contribute a capacity increase sufficient to cause 
the deviation from linearity of the dielectric 
curve just below saturation. This requirement is 
quite consistent with the thickness of adsorbed 
films reported by various workers, \lcHaffie and 
Lehner® for example, for glass near room 
temperature. 

The second experiment was carried out with a 
“condenser’’ designed to emphasize any effect, 
either conductivity or capacity, of an adsorbed 
film. The actual piece used has been described by 
Roberds,'' and consists of a series of conducting 
rings fitted tightly around a large Pyrex tube, 
alternate rings being connected together. The 
total capacity of this element, including that due 
to solid dielectric, was approximately 60 wuyf. 
Calculations involving the number, length, and 
spacing of these rings show that for a given 
thickness of adsorbed film this element should 
show an increase in capacity 324 times the 
increase which would be caused by the same 
thickness of film on the Pyrex insulators of the 


'S McHaffie and Lehner, J. Chem. Soc. London 127, 
1559 (1925). 


condenser used for dielectric measurements. 
Incidentally, it would also have 324 times the 
leakage. The increase in capacity of this element 
due to the presence of water vapor at 23.3°C was 
measured for a great number of pressures ranging 
up close to saturation. Measurements were made 
both with and without the permanent artificial 
leak of 25,000 previously referred to, and the 
results were entirely independent of its presence, 
indicating that even here the change in capacity 
observed was affected scarcely at all by any effect 
of conductivity. It is interesting that the re- 
sultant curve of capacity increase versus pressure 
of water vapor was quite similar to the curves 
obtained by plotting the number of molecular 
layers of water adsorbed versus the pressure from 
McHaffie and Lehner’s™ data. As a_ specific 
example, the increase in capacity was 11.0 uuf for 
a water vapor pressure of 1.8 cm of Hg. This 
increase may rightfully all be attributed to the 
film, for calculation shows that the vapor, with 
its observed dielectric constant, could account 
for less than 0.02 uuf of the increase. If the 
capacity due to the film here is divided by 324, it 
appears that the same film present on the Pyrex 
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insulators of the condenser used for dielectric 
measurements should cause a capacity increase of 
0.034 wuf. The actual deviation from linearity of 
the dielectric curve at 23.3°C and 1.8 cm of Hg 
pressure, under which conditions this identical 
film exists, represents (from curve 2 of Fig. 1) a 
capacity increase of 0.030 wuf. There can be no 
doubt that the added polarization due to the 
actual film of water adsorbed on the insulator 
surfaces will account quantitatively for the 
breaks observed in the dielectric curves for water 
vapor. The author suggests that probably all 
such breaks observed for various vapors are due 
mainly to this same effect. He has actually found 
that the addition of adsorbing insulator surface 
as mentioned above emphasizes the breaks for 
ethyl alcohol and for toluene, though no quanti- 
tative calculations have been attempted. 


On THE MOMENT OF THE WATER MOLECULE 


Dielectric data have been taken at 14 different 
temperatures from 21.3° to 197.9°C, using quartz 
insulators. For the sake of clarity only half of 
these have been shown in Fig. 3. p’ represents the 
pressure the vapor would exert were it an ideal 
gas; it is obtained from the observed pressure p 
and the van der Waals constants a and 3, thus 


p’=p\{1+(p/RT)(a/RT—5)}. 
The molecular polarization is then given by 


(K —1)M/(K+2)d=(K—1)RT/(K+2)p’ 
=A+B/T, 


the last equality being a result of the Debye 
theory. The constant A should be equal to the 
molecular refraction at infinite wave-length, and 
B is related to the permanent electric moment u 
of the molecule by 


u=1.27210-2%/B. 


The molecular polarization at any temperature 
can be determined from the slope of the curve in 
Fig. 3. The best slopes of these lines were 
obtained by least squares, using all experimental 
points below the breaks. 

Fig. 4 shows the product of molecular polariza- 
tion and absolute temperature plotted against 7°. 
Previous data by Jona,' Zahn? and Sanger" are 
indicated. The present data, including that at low 
temperatures, fall accurately on a straight line as 
demanded by theory; no point deviates by more 
than one percent. Sanger’s data, though confined 
to the upper half of the temperature range here 
investigated, fall quite close to the present curve. 
Zahn’s data, for which he makes no claim of 
accuracy, are not inconsistent with the present 
unless it be at the lower temperatures. A least- 
squares solution for the best straight line through 
the author's data gives 


(K —1)RT/(K+2)p’ = (4.03 +0.39) 
+ (20,710+140)/T. 


4 Sanger, Physik. Zeits. 31, 306 (1930); Helv. Phys. 
Acta. 5, 200 (1932). 
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The least-squares value for the electric moment 
of the water molecule is (1.831+0.006) x10 > 'S 
e.s.u. The polarization 4.0 contributed by other 
than the permanent moment agrees within the 
probable error with the molecular refraction 3.7. 
In view of the range of temperature over which 
measurements were made, and in view of the fact 
that not one of the 14 points deviates by more 
than one percent from the Debye line, the 
moment here found should be quite reliable. 
Previously reported values have varied con- 
siderably. Jona! reported 1.87 X 10~'8, from early 


AND H. W. WEBB 

vapor data. Stuart’ reported 1.7910 "8, Wil- 
liams'* (1.70+0.06) 1078 by dis- 
solving water in benzene. Sanger" obtained 
(1.842 +0.008) x 10-'8 from vapor data. Sanger’s 
value appears the most reliable of the group, 
though it was calculated using only three points 
all within a temperature range of 60°. Sanger’s 
value and the present value, obtained from 14 
points covering a temperature range of 177°, 
agree well within the sum of the probable errors. 


obtained 


% Stuart, Zeits. f. Physik 51, 490 (1928). 
16 Williams, Physik. Zeits. 29, 204 (1928). 
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A study has been made of the behavior of a plasma in 
which no new ions are found, as regards electron tempera- 
ture and concentration, to determine the sources of the 
electron energy. Measurements were made with a movable 
probe in a stream of ionized mercury vapor passing out 
from an arc discharge. The velocity and vapor concentra- 
tion were determined by stroboscopic measurements on 
moving striations produced by an impressed alternating 
voltage. The vapor concentration, expressed in pressure 
reduced to 0°C, varied between 0.05 and 1.5 mm of mer- 


INTRODUCTION 


NVESTIGATIONS of the velocity distri- 
bution of free electrons in the plasma of an 
arc discharge! and in the plasma which diffuses 
out of the direct discharge as the afterglow, show 
that the Maxwellian law is followed, so that the 
distribution may be stated in terms of tempera- 
ture. The conditions in the afterglow have been 
studied by a number of investigators? and the 
rate of recombination has been found to be 
11. Langmuir and H. M. Mott-Smith, Gen. Elect. Rev. 
27, 449, 538, 616, 762, 810 (1924); 1. Langmuir, Phys. Rev. 
26, 585 (1925); H. M. Mott-Smith and I. Langmuir, 
Phys. Rev. 28, 727 (1926); T. I. Killian, Phys. Rev. 35, 
1238 (1930); R. Seeliger and R. Hirchert, Ann. d. Physik 
11, 7, 817 (1931). 
2F. L. Mohler, Phys. Rev. 31, 187 (1928); Mohler and 
Boeckner, Bur. Standards J. Research 2, 489; 3, 303 (1929); 


H. W. Webb and S. C. Wang, Phys. Rev. 33, 329 (1929); 
H. W. Webb and D. Sinclair, Phys. Rev. 37, 182 (1931). 


cury. The rate of cooling of the electrons was very rapid 
at first, but below 2500°K this rate decreased rapidly. 
Finally the electron temperature reached a constant value 
which was several hundred degrees above that of the vapor, 
the difference being proportional to the vapor concentra- 
tion. The most reasonable source of energy necessary to 
keep up the electron temperature seems to be either the 
local recombination process, directly or through the meta- 
stable atoms formed, or the metastable atoms brought 
down the tube by the vapor stream. 


closely connected with the electron temperature. 
The variation of the electron temperature in the 
afterglow indicated that the electron energy was 
markedly affected by the various conditions in a 
plasma in which there was present no new 
excitation or ionization. The purpose of the 
present investigation was to study the important 
factors in more detail by observing under varying 
conditions the time-rate of decrease of the 
electron temperature and concentration in the 
afterglow. 


EXPERIMENTAL 


The discharge tube used is shown in Fig. 1. 
Between the hot cathode* C and the anode Q 


3 For this we are indebted to the General Electric Vapor 
Lamp Co. 
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to the axis of S. 











Fic. 1. Diagram of the experimental tube. 


was an arc discharge of 5 to 10 amperes. Mercury 
evaporating from the well W, heated by the 
furnace F, caused a rapid flow through the arc 
and down the side tube S (about 30 cm long) to 
the condenser D, where a reflux tube brought the 
mercury back to the well. The streaming mercury 
carried with it a large quantity of electrons and 
positive ions, filling the side tube with a moving 
plasma. S was super-heated so as to prevent con- 
densation of mercury except at the condenser D. 

The probe P was of tungsten wire, 0.0025 cm 
in diameter and 4 mm long, mounted at right 
angles to the axis of the tube. It could be slid 
along the axis by means of an external electro- 








accumulation of gas. 


wn 
wn 


magnet, and was mounted eccentrically so that 
by rotating it, it could be moved at right angles 


A nickel mesh G was mounted in S, several 
centimeters from the edge of the arc discharge. 
A magnetically controlled gate in G made it, 
- possible to move the probe into the arc. 

To measure the quantity of mercury passing 
through S in a given time, a movable steel ball B, 
controlled magnetically, was used to cut off the 
reflux tube, so that the condensed mercury was 
collected in the small tube U, of known volume. 

During the measurements the tube was con- 
nected to a mercury diffusion pump to prevent 


The vapor velocity and concentration in S 
were measured, using the method designed by 
Ricker. An alternating potential of 25 to 50 
volts, supplied by a vacuum tube oscillator, was 
impressed between the anode Q and the mesh G, 
which was biased so that it was always slightly 
negative with respect to Q. As G was made more 
negative to Q, the concentration of the electrons 
and positive ions in its vicinity was diminished, 
accompanied by a decrease in the afterglow in 


+N. H. Ricker, Phys. Rev. 17, 195 (1921). 
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Fic. 2. Electron temperature plotted against time, as measured from the arc. 
Vapor concentrations in atoms/cm*: curve (a) 0.7910"; curve (b) 1.9210"; 
curve (c) 4.91 X10"; curve (d) 7.8810". Vapor velocities in cm/sec.: curve (a) 
5200; curve (b) 1300; curve (c) 464; curve (d) 300. 
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546 R. H. RANDALL 
the vapor flowing through the mesh. As a 
consequence, with the alternating potential, the 
afterglow was broken up into moving striations, 
the spacing of which was measured with a 
stroboscopic wheel. From the frequency of the 
alternating voltage the velocity of the vapor was 
determined. From the measurement of mass flow 
already described the concentration of the vapor 
in the tube was then calculated. The striations 
were fairly plane and evenly spaced, showing 
that for any given conditions the vapor concen- 
tration and velocity were uniform throughout 5S. 

By controlling the difference in temperature 
between the well and the condenser, vapor 
velocities between 300 and 7000 cm/sec. were 
obtained. By raising or lowering both well and 
condenser temperatures, vapor concentrations in 
S were obtained corresponding to pressures, 
reduced to 0°C, of 0.05 to 1.5 mm. The electron 
concentration was controlled independently of 
pressure and velocity conditions and of the 
distribution of electron temperature by varying 
the are current. 

The electron temperatures were measured by 
using the probe current-voltage characteristics, 
according to the methods of Langmuir and Mott- 
Smith. The electron concentrations were ob- 
tained from the probe currents at the “break- 
points.” 


RESULTS 


The percentage variation of the electron tem- 
perature over any one cross section of the tube S 
was found to be small. Therefore the results given 
are temperatures taken on the axis. In Fig. 2 are 
shown four typical electron temperature-time 
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Fic, 3. Final electron temperature level 7; plotted against 
vapor concentration in atoms/cm*. 
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curves for four different vapor concentrations. 
The electron temperature 7°, is plotted against 
elapsed time ¢, as measured from the edge of the 
are discharge. The corresponding vapor concen- 
trations and velocities are given below the graph. 
Since the vapor concentration could not be 
increased without decreasing the vapor velocity, 
the total time covered by each curve is different. 
Each curve includes data taken over the full 
length of the straight portion of the tube, S. 

In the early part of the curves 7, falls rapidly. 
In the later portions, 7, reaches a constant value 
7T,;, considerably higher than the temperature of 
the vapor 7,. For all these curves 7,=480°K; 
while in curve }), 7;=550°K, in curve c, 700°K 
and in curve d, 975°K. 

The vapor velocity was not 
enough to show the final flattening out of the 
7,.—t curves. Fig. 3 shows the final electron 
was clearly 


always slow 


temperature level 7;, where it 
reached, plotted against vapor concentration C. 
The relation between C and 7; is approximately 
linear, the scattering of the points being ac- 
counted for by the small range of temperatures 
involved and the difficulty of making precise 
temperature measurements. For zero concen- 
tration 7;=470°K, which was very closely the 
temperature of the vapor 7,. 

The value of 7,, i.e., 480°K, in the region 
where 7°, = 77;, was taken as the wall temperature. 
This was checked in the following manner. 
Because of the low heat capacity of the gas and 
the streaming of the vapor, the temperature 
within the vapor could not be measured directly. 
A thermocouple was therefore placed on the 
outside of the wall of tube S in this region, and 
that portion of the tube lagged. No rise in 
temperature, such as would be expected if the 
vapor were hotter than the wall, was obtained. A 
computation of the rate of cooling of the vapor 
as it flowed down the tube showed that at this 
point of the tube the difference between 7), and 
the wall temperature was less than the experi- 
mental error of measuring 7’. 

The measurements of the electron concen- 
tration were the most uncertain. Over the length 
of S the change of concentration varied between 
100 and 10,000-fold, according to conditions. As 
the vapor moved down the tube, the electron 
concentration on the axis, mo, fell off rapidly at 
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ELECTRON TEMPERATURE 


first; near the arc there was, in many cases, a 
1000-fold change within a few centimeters. The 
decrease in m») continued throughout the whole 
length of S, but at a continuously decreasing 
rate. The electron concentration was nearly 
uniform over any one cross section of S. Half-way 
between the axis and the wall the electron 
concentration had fallen off less than 20 percent. 
This is of importance in the later discussion. 


DISCUSSION 


The most interesting characteristic of the 
results is that 7, drops to a constant value 7; 
which is above the gas temperature. Since the 
factors determining the cooling of the electrons 
which are effective in the falling part of the 
7.—t curves must still be operative in this 
region, there is evidently an important mecha- 
nism in the plasma tending to raise the tempera- 
ture of the electrons. The complete solution of 
the problem presented here is difficult because of 
the variation of the electron temperature along 
the axis and because the mechanism of the 
energy transfer by the electrons and of the 
maintenance of the Maxwellian distribution is 
not understood. However, the measurements give 
the distribution of electron temperature and 
concentration over the greater part of the tube S, 
so that the distribution of space potential, wall 
potential and current to the wall can be approxi- 
mately computed. 

In making these computations use was made 
of the theory worked out for the plasma by 
Tonks and Langmuir.’ Let 7 and z be the radial 
and axial coordinates and let R be the radius of 
the wall. The space potential in each cross- 
sectional element is taken as zero on the axis. 
Let n=el'/k7T.=11,600 V/T., where e is the 
electronic charge, |’ the space potential in volts 
for any value of 7, and k the Boltzmann constant. 
For any value of 7, the concentration of electrons 
or ions is m= nye”, where my is the concentration 
on the axis. 

From kinetic theory the electron current 
density to the wall is 

T,=enye™ (RT. /22m.)', 


where 7,, is the value of » at the wall and m, is the 
electronic mass. If it is assumed that the loss of 


°L. Tonks and 1. Langmuir, Phys. Rev. 34, 876 (1929). 
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IN MERCURY AFTERGLOW : 
number by recombination is negligible, ° 


0 7 0 
Ia(2R)=eo— f nye—"2rrdr = ev—(A ny) rR?, 
Oz / 9 Oz 


where v is the vapor velocity, and Amy is the 
average electron concentration in any one 
section. 

To determine the value of », (and hence the 
wall potential) from the above equations, it is 
necessary to determine A. Since the electron 
concentration was found to be nearly uniform 
across the tube for all values of z, A certainly lies 
between 1.0 and 0.5. If A=1, the maximum 
possible value, n,, is 12.1, 12.5, 15.8, and 11.6 for 
curves a, b, c, and d, respectively, in Fig. 2, while 
for A=0.5, yn» has values smaller by 0.6. This 
variation in , is negligible for the purposes of 
these computations. 

The variation of n,, with s was found to be very 
small. This, together with the fact that e~” varies 
only slowly with 7, until near the wall, justifies 
for a first approximation the assumption that 
€” is independent of z. This gives a condition 
similar to case 4 on page 889 of the paper by 
Tonks and Langmuir. In this the potential 
distribution was determined in a plasma, with 
short ionic mean free path and ion generation 
proportional to electron density. The generation 
of new ions is replaced here by the net flow in the 
z direction of electrons and ions, vd" /dz, into any 
element, since if «~” is a function of 7 only, 


on 0 1 dno 
v— =v—(nNoe—") =v— ——n=constant n. 
2 oO ny ds 


The solution* is 


1 dny : 
era {o(— 7 ) Jara} r| 
No dz 


where g=0.895 1,/(k7T,m,)', l, being the molecu- 
lar mean free path and m, the mass of the 
positive ion. 

As a check on the accuracy of this expression, 
the value of €~ determined above and experi- 
mentally determined values for the variables in 
the argument of J») were substituted, and the 
equations solved for R. The average value found 
6 kK. T. Compton and I. Langmuir, Rev. Mod. Phys. 2, 


197 (1930). 
* Reference 5, cf. equation at bottom of p. 889. 
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was approximately 2 cm, the discrepancy be- 
tween this and the actual value, 1.77 cm, being 
accounted for by the precision in the measure- 
ment of (1/n)dno/dz. Using the expression for 
e”, the value of A was found to be 0.95. 

It must be emphasized that all the above 
results are based on the assumption that a 
Maxwellian distribution is always maintained 
throughout the space in question. The probe 
current-voltage curves taken here give the 
energy distribution of the electrons only up to a 
few tenths of a volt and any departure from the 
assumed distribution for higher energies could 
not be detected. 


e107 





aR 


oz any 
E,= | | r=—| | 
v2 dz d/o 


The first part of E, is the work done by the 
electrons in reaching the wall against the 
retarding potential; the second part, which is 
small in comparison (about 10 percent of the 
whole), is the energy carried out by the electrons 
as kinetic energy. Neglecting » with respect to 
Nw, since 7 is nearly zero except near the wall, this 


11,600 


equation becomes 


e107 3 le 
Be= Arka — meth) | T dno. 
11,600 2 ny 


By substituting experimental values in the equa- 
tions for E. and E,,, it was found that in all cases 
the loss to the wall was much greater than the 
energy brought in by convection. For the data of 
curve a, Fig. 2, E.=435 ergs/sec. and F,.= 1500 
ergs/sec. for a point z near the arc. For a point 
near the condenser, E.=7 ergs/sec. and E,,=38 
ergs/sec. For the data of curve'c, where the vapor 
velocity was only 0.1 as great, near the arc 
E.=0.25 ergs/sec. and £,=1.0 ergs/sec.; near 
the condenser E,=9.1X10~ ergs/sec. and Ey 
=73X10~ ergs/sec. 

To account for this difference in energy AE, 
several sources must be considered. The more 
obvious are the transfer of energy to electrons by 
metastable atoms swept down the tube S or 
produced throughout S by the recombination 
process, the recombination process itself, and the 
transfer of energy in the z direction between the 
electrons. 
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The nature of the mechanism of the cooling of 
the electrons as they pass down the tube can be 
best discussed by considering the total changes in 
energy in any fixed volume. Calculations were 
made of the energy gained and lost in the whole 
of a cylinder formed by that part of the tube S 
lying between a transverse section at z and a 
second section at z,, at the lower end of S (see 
Fig. 1). The net energy brought into such a 
cylinder by convection is E.=(3/2)kT vA[no 
— (no), |r R? ergs/sec., where my and (mo), are the 
electron concentrations on the axis at z and 2,, 
respectively. The energy carried out per second 
by the electrons escaping to the wall is 





ee 3 dno R 
—(e-*)(ne— 9) 2erdr| Nast f [er (f ari) f= 
2 v4 ds ee, 





The magnitude of this last process, which by 
analogy can be referred to as conduction, is 
difficult to estimate, since it involves the mecha- 
nism of the reestablishment of the Maxwellian 
distribution, which is not clearly understood. 
When the electron temperature is falling rapidly, 
as in the early part of the flow, conduction is 
probably important. If it is assumed that all of 
AE is due to this factor and if the small quantity 
of energy leaving the cylinder 22, at the lower end 
is neglected, since any process of this type must 
depend upon the concentration at the section 
across which the energy transfer is taking place, 
it is convenient to compute AF/Anp. Starting at 
the edge of the arc discharge, the value of the 
quotient falls rapidly with the temperature, 
decreasing to approximately one-fourth its initial 
value when the final temperature level 7; is 
attained. The fact that this quotient decreases 
with the electron temperature and with the 
temperature gradient is consistent with the 
assumption that such a process is important in 
the region of falling temperature. In the region in 
which the electron temperature is constant, 
however, there should be no energy transfer of 
this type, since the electrons are all at the same 
temperature. Here AE must be due to other 
sources of electron energy. 

The probability of the transfer of energy by 
metastable atoms to electrons on impact is high. 
No estimate of the concentration of metastable 
atoms or of the energy which they might give to 
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the electrons has been attempted, since no data are 
available for such a computation. Furthermore it 
is not known to what extent these metastable 
atoms are produced locally by recombination and 
to what extent they are swept down by the 
vapor from the region of the arc. Webb and 
Wang? found that both factors are important. 
They further showed that if into a moving 
plasma similar to that studied ‘here, they intro- 
duced sodium vapor, a marked drop in electron 
temperature immediately resulted; in the case 
cited from 2300° to 800°K. The most important 
action of the sodium on the flowing vapor was to 
eliminate the metastable atoms by excitation of 
the sodium, and the most probable explanation 
of the temperature drop is that the elimination of 
the metastable atoms removed a source of elec- 
tron energy, resulting in lowered electron 
temperature. 

Webb and Sinclair? found that the rate of 
recombination depended to a large extent on the 
electron temperature, and came to the conclusion 
that a somewhat complex interaction between 
ions and fast and slow electrons was involved in 
the recombination process. Furthermore, a study 
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of the intensity distribution in the afterglow 
spectrum leads to the conclusion that consider- 
able kinetic energy may be freed by the process. 
It therefore seems probable that the recombi- 
nation process is largely responsible for the 
energy supply which tends to hold up the electron 
temperature, either through some process in- 
volving many-body collisions, or by the meta- 
stable atoms which constitute a large part of the 
atoms formed by recombination. It may be noted 
that the recombination has been found by Webb 
and Sinclair and by Mohler and Boeckner’ to 
increase with the pressure, which corresponds to 
the increase in the final electron temperature 
level with increase in pressure found here. On the 
other hand, for a given rate of production, the 
concentration of metastable atoms also tends to 
increase with increased pressure, as under these 
conditions the life of these atoms is practically 
determined by the diffusion to the wall. 

The authors wish to express their appreciation 
to Mr. T. C. Hardy for his assistance in the 
taking of the measurements. 


7 Mohler and Boeckner, Bur. Standards J. Research 2, 
489 (1929). 
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In order to put Zwaan’s scheme for deriving the Kramers 
connection formulas on a rigorous basis, differential equa- 
tions are set up governing the variation in the coefficients 
used to fit a linear combination of the B. W. K. type ap- 
proximation functions to an exact solution of Schrédinger’s 
equation in one dimension. Approximate solutions of the 
differential equations are worked out which lend them- 
selves to the setting up of the connection formulas and 


1. INTRODUCTION 


Zwaan’'s' scheme for deriving the Kramers 
connection formulas? used in the Brillouin- 
Wentzel-Kramers (B. W. K.) method is simple 
and illuminating, but not entirely rigorous. In 


1A. Zwaan, Intensitaten im Ca-Funkenspektrum, Dis- 
sertation, Utrecht, 1929. 

2H. A. Kramers, Zeits. f. Physik 39, 828 (1926); H. A. 
Kramers and G. P. Ittman, Zeits. f. Physik 58, 217 (1929). 


give definite upper bounds to the errors involved in their 
use. The method is also used to set up the Sommerfeld 
phase integral quantum condition independently of the 
connection formulas. An upper limit to the error in the 
energy is worked out. A similar treatment of the problem 
of the transmission of matter waves through rounded po- 
tential barriers is formulated. 


the attempt to put Zwaan’s argument on a 
rigorous basis the author has developed a new 
way of studying the relation of the Stokes 
phenomenon to the B. W. K. method whose value 
extends beyond the scope of the connection 
formulas themselves. 

Consider a one-dimensional problem in wave 
mechanics with the potential energy )’(z) and the 
classical local momentum 
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p(s, E)= \wlLE-— V(s)]}}. 
The Schrédinger equation takes the form 
v" +X py=0 (1) 


with A set equal to 27h. The basic approximation 
functions of the B. W. Kk. method are 
fi(2, E)=p(2, E)'e"* *); 
; — (2) 
J. (2, E)=p(2, E)~ te" *); 


w(s, &) =f p(t, E)d¢. (3) 


Here « denotes ¥ —1. The point z9 can be chosen 
at pleasure but is ordinarily placed on the axis of 
reals. f; and f, are readily shown to be solutions 
of the differential equation 


f’ +P —-QIf=0, (4) 


where Q denotes the function 


379.2? gp” WSs V’ yy? ~ «2: 

AS) pata tet © 

4\ p 2p 4L4\E-1 E-| 
Wherever |Q|<)*p?, Eq. (4) is very nearly the 
same as (1) and f;, f, may be said to be good 
approximate solutions of (1). 

Let zs’ denote a simple real zero point of the 
function E—V(s) (primes and double primes 
denote differentiation except when applied to the 
symbol z itself). For definiteness, we assume that 
E— V is negative to the left of z’ and positive to 
the right (cf. Fig. 1). At 2’, fi, f, and Q all 
become infinite. Hence f;, f, have no value as 
approximate solutions of (1) in the neighborhood 
of z’. We shall suppose, however, that there exist 
regions I and II to the left and right of 2’, 
respectively, where |Q|<)?p? and where the /'’s 
are good approximate solutions. Then if ¥(z) is an 
exact solution of (1), we may reasonably assume 
that in any not too large portion of I, ¥(z) is 
approximately equal to some suitable linear 
combination of f; and f,. Thus, if a; and a, are 
suitably chosen constants, 


afitaf.=yY(z) in region I. 
Similarly 
in region II. 


Bifit+B.fo=V(2) 


It is well known that the coefficients 6;, 8, are 
not equal to aj, a,, respectively. The shift in the 
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values of the coefficients is called the Stokes 
phenomenon and the problem is to determine the 
magnitude of this shift in terms of the values on 
one side of 3s’, say in I. 

Kramers method of attack is based on the 
assumption that I(s) is sensibly linear in the 
interval between I and II. Under these circum- 
stances ¥(s) can be approximated by a suitable 
Bessel’s function over the interval in question 
and this function constitutes a bridge between 
the B. W. K. approximations in I and II. Langer* 
has introduced a different type of approximation 
function involving a Bessel’s function of the 
argument w(x) defined in Eq. (3). The Langer 
approximation function degenerates into the 
B. W. K. type in I and II, but does not go to 
pieces near 3’. This method is very powerful, but 
has some disadvantages with respect to that 
developed in the present article. 

Zwaan’s derivation of the connection formulas 
follows closely a line of thought previously 
developed by Stokes.‘ It is based on the intro- 
duction of complex values of the independent 
variable and the existence of a path IT joining I 
and II, but passing around 3s’ in the complex 
plane so that the inequality |Q <)°p* is good 
along its entire length. If I(s) is analytic on the 
axis of reals between I and II, and if its analytic 
extension into the upper half plane yields no 
roots of V(s)—F near zs’, we may assume that 
such a path can be found.® 

The essential feature of the situation to which 
Zwaan calls attention is that along a portion of 
any such path I the approximation f; becomes 

3R. E. Langer, Trans. Am. Math. Soc. 33, 23 (1931); 34, 
447 (1932); Bull. Am. Math. Soc. 40, 545 (1934). 

4Sir George Gabriel Stokes, Math. and Phys. Papers, 
Vol. 4, pp. 77-109 and 283-298. 

5 Let us assume that Il’(s)—£ is a polynomial of degree 


n, or can be accurately represented by such a polynomial 
in the interval between I and II. Then, if zs; denotes the 


n 
kth root, E—V=CII (s—2,). It follows that 
bel 


n 
(E-V)' /(E-V)= 5 —_:; 
k= 


12-23% 








> W\" CR Vie E—v"\? § 1 
(E-—V)"/(E-V) (a7) = - 
Hence 


$- enteo (Sata) Feta] 
P 8uCI,(s—2x) 4\ p21 2-2 k (2—2%)? 


Inspection of this expression for Q/p? proves the above 
statement for the case that is small. 
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very small, while f, becomes very large. Hence f, 
and f; are said to be dominant and subdominant, 
respectively, in the region in question. Where this 
is true changes in the subdominant coefficient can 
take place even when |Q|<)*p?. On the other 
hand the coefficient of the dominant approxi- 
mation function cannot change much, since a 
small change in this coefficient makes a very 
large difference in the function. Zwaan assumes 
that the latter coefficient is constant. 


2. THe DIFFERENTIAL EQUATIONS FOR 
THE COEFFICIENTS 


Granted that the changes do take place, we 
propose here to investigate them by setting up 
explicit differential equations for the variation in 
the coefficients. The formulation of these equa- 
tions, which forms the essential contribution of 
the present paper, is reminiscent of a method of 
attack on one-dimensional problems in quantum 
mechanics due to Hill.® 

Let ¥(s) be an exact solution of (1). Let 
u=a;f;+a,f, be a linear combination of f; and f, 
which is fitted to y at the arbitrary point z=). 
The best fit will evidently be obtained if u and 
du ds are made equal to ¥ and dy/dz, respec- 
tively, at the point s;. In that case 


ai fit) +a,f.(s1:) =W(21); (6) 
aif; (s:) ta,f.’ (21) = (21). (7) 


As the Wronskian determinant /,f,’—f,f;’ does 
not vanish (it has the constant value —4m/h 
= — 2X.) it is possible to choose a; and a, so as to 
satisfy these equations at an arbitrary point 2, 
where f; and f, are uniquely defined. The 
coefficients so obtained will then vary continu- 
ously as 2; moves around the path I’, reducing to 
a; and a, where [ meets the axis of reals in I and 
to B;, 8, where T' meets the axis of reals in II. 

Dropping the subscript 1 in Eqs. (6) and (7), 
solving for a;, and differentiating with respect to 
z, we obtain 

da; 


“= (0/29) (fof). (8) 


ads 


On reduction with the aid of Eqs. (1) and (4) 


*E. L. Hill, Phys. Rev. 38, 1258 (1931). 
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da; 1 . Q 
—_—= —Of wv =— —[aite aug ok (9) 
dz 2x 2X p 
Similarly 
da, 


t . Q 
—= —-—Ofw= —— -[a,+e?"a;]. (10) 
2r 2d p 


az 


These differential equations are exact. Their 
integration is equivalent to the integration of Eq. 
(1), although they blow up at the zeros of 
p(s, E). If an exact integration were possible, it 
would be a simple matter to work out the exact 
eigenfunctions and eigenvalues of (1). The pri- 
mary usefulness of Eqs. (9) and (10), however, 
lies in the possibility of deriving useful approxi- 
mations from them. 


3. THE APPROXIMATE INTEGRATION OF EQs. (9) 
AND (10) 
Introducing a fixed point £, let us make the 
transformation 


b(z, &) =a,(z)e~*F(* &). 
(11) 
b,(z, £)=a,(z)et tF(2, &) 


1 77Q 
F(z, p-—f —dz. (11a) 
207 p 


This change of variables leads to the following 
simpler system of differential equations: 


where 


F ) 
db, on tates by; 
dz 2x 
. (12) 
db, 
umn n. Gan canta —e2(F+w)h., 
dz 2d p 


Clearly 6; and b, will be sensibly equal to a; and 
d,, respectively, at any point 3 which can be 
connected with ~ by a path along which |Q/p! is 
small enough. It is convenient to define an index 
of quality uw, for a path A in the complex plane by 
the equation 


1 | Q| iri?) . 
ws=—| \ds= f lw}. (13) 
Mal P| Nos pF 


We shall hereafter refer to a path as good, provided 
that its index of quality is much less than unity. 
Any short path which gives a wide berth to the 


~ 
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zeros of p(s, /) can ordinarily be assumed to be 
ae 2 aed . “* 
good”’ in this sense. Since 


ae 
IX. ge 


it follows that | F(s, £)| <}ua if s and & are any 
two points on A. 6; and b, are thus sensibly equal 
to a; and a, at any point z which can be con- 
nected with & by a good path. 

It will simplify the discussion of the integration 
of Eqs. (12) if we assume at first a definite path of 
integration A, of total length /, leading from an 
initial point Z» to a final point Z;, and having the 
property that the imaginary part of w is a 
monotonically decreasing function of the path 


< MAs 
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assume that |e’| increases monotonically as we 
move from Zp» to Z,; along A. With these re- 
strictions we introduce the real variable of 
integration s and the notation 


10(s)/2dp(s) =X(s), 


w(z)=w/(s), 


(14) 
dz/ds=g(s), F(z, £)=6(s, &). 
Eqs. (12) take the form 
db;/ds = X(s) ¢(s)e—2""*H,, (15) 
db,/ds = —X(s)¢(s)e***tD;. (16) 


Under favorable circumstances 6; can be com- 
puted approximately for points on A by replacing 
b, in (15) by its initial value 6,°=0,(Z»). Let 
b,(s) denote the corresponding approximation for 





length s measured from Z». In other words, we 


bj(z). Then 





bi(s) =b +b," fy X(sie(sve 2s toon ldsy, (17) 
This approximation can be tested by computing an upper bound on |6;(Z;) —),(Z,)|. We have 
b:(Z:) —b(Z1) = fy'[bo(s) —b.° JE(s) p(s)em* ds. (18) 
But from (16) 
by(s) —b,°= — fob i(si)X (si) e(sie* Cov ten Idsy, (19) 
Hence 
b(Z:) —b(Z)) = — f'ds (B(s) o(s)e- 20 +1 B51) EB (81) o(s en OY Hed} 
(20) 


= — SS rbi(s1)=(Si)=(s) e(S1) o(s eB PCD Fo) “2 Ids ids, 


where 7 is a triangle in the s, s; plane defined by the inequalities 0<s;<s; O&s<1. Since |e”| 
increases monotonically along A, |e~*!-#*)]| C1 at every point in 7. ¢(s) is of unit absolute 
value and 


| e~ 2+ (6¢5)—OCa1) J | = le ( ») Say LOC") [68] 908") ds" | < e#A, 


Introducing the symbol ./; for the maximum value of |);! on A, we readily derive the inequality 


|b:s(Z1) —b,(Z) | <Mie*s Sfr| E(s1) | (21) 





D(s) |dsids. 


The function | ¥(s)| |2(s:)| being symmetrical with respect to its arguments, the above integral over 
T is half the corresponding integral over the complete square 0< s<1; 0&5,<1. Thus 


|bs(Z1) —b:(Z1) | <3. Mye#s{ f"|E(s) |ds}2=pa2e*A(M;/8). (22) 


We conclude that in the case of a good path leading uphill with respect to | e‘’| the approximation (17) 
will lead to an error which is small in comparison with the maximum value of |;| on the path. 

Let wo denote the value of w at the point Zo. The maximum absolute value of the factor e~***” of 
the integrand of Eq. (17) is |e~*““*|. Furthermore on the path A|e~??| < e?!@! <e#0t#4 where uo denotes 
the index of quality of the best path leading from £ to Zp. We assume that yo is a small quantity of the 
order of wa. Eq. (17) yields 


|b(Z:) —b| <Suserot#al b,°| |e Zuo | ' 





e 
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Combining this inequality with (22), we obtain 
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|b:(Z1) | <|b,°| + 3uaenot#a| b,°| |e-2"| + wa2eea( M;/8). (23) 


Let Z,, denote the point on A where b;(z) becomes equal to M,, and let ya’ denote the quantity 
1/ASz,7"|Q/p|ds, where the integration is carried out along A. Then ua’ Cua. Clearly an equation 


similar to (23) applies to the point Z,,. Thus 


Mi < |bi°| + (ua’/2)emt# 4" | b,9| |em2"9| + ye44'(M,/8). 


Transposing and neglecting terms of higher power in the small quantities uo, wa, we obtain M,;< |6,°| 


1 | 
2KA | 


b,°| |e-*e|. It follows that in this approximation 


|bi(Z1) —b:| < | b(Zs) —b(Z1) | + |b:(Z1) —B:| <3 ua] b.| [em2"| +3 ua? d,°|. (24) 


Consider next the change in a;. We have 
las(Z1) —a,"| - | erF (Zo. é) 


where 


\b(Z ie 2rSi Q pdz_ h 0| ; 


[1 7 @ | 
| F(Z, €)|<}u0 and \- | oe < ma. 
A 


Hence, neglecting small quantities of higher order, 


|ai(Z1) —a,°| <ere!? fens’? | bi(Z1) —B,°| + | b,°| (ua/2)} <dmal|Bi°| + | b.°| |e-2”9| J. 


Finally, it follows from (11) that to this same approximation | b,°| — |a,°| < |a,°| (uo/2). Consequently 


|ai(Z1) —a°| < dual |ai®| +] a,°| |e?! J. (25) 


We conclude that the variation in a; is small in comparison with the larger of the two quantities 


la;°| and |a,°| |e-29!. 


In applying the above result two special cases 
are to be considered. The first is that of a good 
uphill path restricted by the additional require- 
ment that |e’| 21 along its entire length. (It 
will be convenient to refer to regions in which 
le’! >1 as mountains, and to regions in which 
e’| <1 as valleys. Clearly f; is dominant and f, 
subdominant on the upper levels of the moun- 
tains. These roles are reversed for the lower levels 
of the valleys.) In this case |e-?"°| <1 and 
'a;(Z;) —a,°| is small compared with the average 
of the quantities |a,;°| and |a,°|. Under these 
circumstances we shall say that a; is sensibly 
constant along the path, although if |a,°|<ja,°), 
it may vary by a large percentage of its initial 
absolute value. 

In order to deal with the variation in a; along 
a good path leading from Zo clear over a moun- 
tain to Z;, but not entering the adjacent valleys, 


it is convenient to think of the path as made up 
of two parts leading to the crest from opposite 
sides. Applying the above argument to each part 
of the path we come to the conclusion that a; is 
sensibly constant over the entire path in com- 
parison with the largest of the four quantities 
|ax®|, |av®|, |as(Z1)|, |a(Z1)|. 

A second important case of the application of 
Eq. (25) is that of a good path leading uphill, but 
beginning in a valley, where |e’| <1 and f, is 
dominant. If the valley is deep |e~**’*| is a very 
large number, but if |a,°| |e~®*"*| < |a,°|, we have 
the inequality 


|ai(Z1) —a;°| <yala,°| (26) 


and can say that a,(z) is sensibly constant along 
the good path A. Evidently a;(z) will also be 
sensibly constant along a good path A extending 
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across a valley from one mountain to another, 
aj(s)| at the 


| 2uw(z) 


provided that |a,(s)| |e")! < 
lowest point of the path. 

Turning to a,(z), we readily derive formulas 
parallel to those for aj,(s) provided that we 
perform the integrations in the opposite direction 
on A, i.e., downhill with respect to |e”. Thus, if 
we take a;° and a,° as the values of a; and a, at the 
upper end of the path, we obtain instead of (25) 
the inequality 


|a.(Z1) —a,°| <3ual|ae®| +|a,°| |e] }. (27) 


In general the coefficient of the dominant 
approximation /; or /, is sensibly constant over a 
good path A on which |e” —1 does not change 
sign. The coefficient of the subdominant term is 
also sensibly constant over such a path provided 
that the absolute value of the dominant term at 
the point of maximum dominance is zero, or 
smaller than the product (evaluated at the same 
point) of the absolute value of the coefficient of 
the subdominant term multiplied by the smaller 
of the two quantities e***”|. However, since we 
have defined the phrase ‘sensibly constant” 
somewhat loosely, it is perhaps best to base all 
rigorous arguments directly on the inequalities 
(25) and (27). 


4. DERIVATION OF THE CONNECTION FORMULAS 


We assume that the potential function in the 
neighborhood of the classical turning point 3’ has 
the form indicated in Fig. 1. Introducing polar 
coordinates r, @ with the origin at 3’, we write p? in 
the form 

p? =2p(E— V)=re"o(z), 
where g(z) is assumed to be bounded from zero 
in a wide region of the complex plane about 2’. 
We further assume the existence of a path 
passing around 2’ in the upper half plane (cf. 
Fig. 3) enclosing no complex zeros of p? and good 
in the technical sense that wa<1. Let P; and 
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Fic, 2. 


P, denote the end points of A located in regions I 
and II, respectively. 

To avoid ambiguity due to the multiple- 
valued character of p, p', and w, we restrict the 
angle @ to the range 0 < @& 27 and define p and p' 
by the equations 

p=rie®loi: plarie®/4o}, 


¢? can be treated as single-valued. Let the point 
zo LEq. (3)] be identified with s’. Then in the 
region I we have 
=lol: ploe*/4l pli: 
p=tip\: pr= |P\*: 
w= —irf.*'| pl dé. 


In the region IT to the right of 2’ 


p=|pls p=|pl!: w=dfet |plde. (29) 





Ul 


A cut extending from 3z’ to the right along the 
axis of reals insures that we have a single 
branch of w(z) to deal with. 

In order to visualize the problem it is con- 
venient to map the level lines of |e‘’|. Eqs. (28) 
and (29) show directly that the axis of reals 
passes through a mountain to the left of 2’, but 
that the portion of the axis of reals to the right of 
3’ is a level line with |e”| =1. Writing s=x+ vy; 
p= | pie, the differential equation for the level 
lines is seen to be 

(p—p*) 
dy/dx=t —=-—tan x, 


(p+p*) 


where p* is the complex conjugate of p. Treating 
¢ in first approximation as constant over the 


(30) 








it 
1¢ 
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Fic. 3. Level lines of |e”! for linear potential function. 
e*”!| is unity along cut and dotted lines radiating from Z’. 
Arrows indicate the uphill direction. 


neighborhood of s’ we obtain the equation dy/dx 
=-—tan3é@ for the neighborhood in question. 
Fig. 3 shows the approximate form of the level 
lines derived from the above differential equa- 
tion. Three lines radiate from the branch point 
s’ and from every other simple zero point of 
E—V(s). The exact forms of the contour lines 
will, of course, depend on ¢g(z), but the distor- 
tion produced by ¢ is sure to be small if ¢ is 
smooth in the neighborhood of 2’. 

We identify the point € of Eqs. (11) with P,, 
thus insuring that the a’s and b’s are sensibly 
equal on A. Let D denote the point where A 
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un 
wn 
uw 


upper half plane. Then a, will be sensibly 
constant between D and P» because f, is domi- 
nant over this portion of A. It will also be 
constant over P,D if a; vanishes, or is very small, 
at P,. It follows that if a;=a;(P;)=0, a, =). 
This fact, together with the reality condition 
gives the first connection formula. 

Let us cast the argument in rigorous form. 
Denoting the point on A where |e’! has its 
minimum value by £ and the corresponding 
values of a; and a, by y; and y,, respectively, we 
apply Eq. (27) first to the arc P,E and then to 
PE. Let uw; and pe be the values of u for these two 
sections of the path A. Then 


| 1 f | | 29ww(P))/| 1 
[ye are| < Bur { [ore] + Jovi] [errr |}, 


lye —Bv| <3ue2{ |Be| +] Bi! }. 


(31) 


Consider the application of these equations to a 
solution of Eq. (1) which is real on the axis of 
reals and which is of the f, type at P;. From the 
reality condition it follows that 8;=8,* and 
|8:|= |8.|. If the function is to be of f, type at P; 
we may assume that |a;| |e?) | < |a,|. Hence 


la»—By| <|yr—a@, + |y»—Bo| 
<pilar| tme!Br| <pafla.|+|B,!}}. 


We can accordingly identify a, and 8, in first 
approximation. 

It follows from (28) that if a, is given the value 
e‘*/4 the product a,f, is real in I and describes a 
real y function. Then 8,=8,;*=e'"/*. Thus the 





crosses the line |e“’|=1 which extends into the connection formula 
az T 
ettl4 f= | p|—le Se" Pi B42 -! cos E pdé—- | (32) 
e/ 2’ 4 
eens . “sd 
Region I Region II 


is established. 

In this formula, following Langer, we draw the 
arrow from left to right to indicate that the 
approximate validity of the left-hand member 
implies that of the right, but that the converse 
statement is not true. In order to make clear the 
justification of this one-way street sign and to 
prepare the way for the derivation of the second 
connection formula, we observe that in virtue of 
the homogeneous linear character of Eqs. (9) and 
(10), the functions a;(z), a,(s) must be homo- 


geneous linear in the constants of integration aj, 
a,. Hence 8; and 8, must be homogeneous linear 
functions of a; and a,. Thus a complete exact 
solution of the whole connection problem would 
involve a complete and exact determination of 
the matrix |g!) of the equations 


Bu = LriGit Lov. (33) 


Our approximate solution of Eqs. (9) and (10) 
for the special case that a;=0 shows that 


Zev=1+6; giv= —e(1+5*), (34) 


B; = Lisi t Livdy, 
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where 6 is a small quantity of the order of magni- 
tude of ua which we have neglected in (32). 

Additional information regarding the matrix 
|g! isobtainable from the fact that the Wronskian 
of any two exact solutions of the Schrédinger 
equation, say y and y¥, is constant along the axis 
of reals.’ Thus if 


Y=aifita.fi; 


in J, it follows that in J 


_ do _dy 
Wy, ¥ ]=¥— —¥— = 2.\(a,a;—a;a,)= constant. 
Oz Oz 


The Wronskian of the same pair of solutions in JJ 
takes the form 


WLy, y j= 2iA(a,ai— QAy) (Lor¥ ii — ZivBvi). 


Equating these two expressions for the Wron- 
skian, we see that the determinant of ||g|| is unity. 
Hence the inverse equations to (33) are 


ay= — £viBit ZiiBo. (35) 


It follows from the reality condition that if 
6:=8,*, the quantities a,e'*/* and a,e‘*/4 are real. 
Hence 


a= LovBi cae LirBo, 


Zoi= — Zii*; Liv= — tLe". (36) 


Combining the last relation with the requirement 
that |g| be equal to unity, we obtain the 
following relation between g;; and 6: 


gii(1+6)+¢:,*(1+6*)=1. (37) 


If we neglect small quantities of the order of a, in 
comparison with unity, we may conclude that the 
real part of gi; is 3. The imaginary part of gi; 
cannot be determined even approximately with- 
out an explicit evaluation of an integral of the 
form given in Eq. (17) over a path on which the 
exponential factor is very large. However, if we 
apply the inequalities (31) to the evaluation of 
|8.| when a,=0 we can derive the upper bound 


az 


lgis|=|gvi] <male*P | =y4 exp| 20 olde 
Pp 


i 


7 The writer is indebted to Dr. Eugene Feenberg for the 
observation that the second connection formula can be 
derived from the first by means of the properties of the 
Wronskian. 
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It does not seem necessary to include the details 
of this simple proof here. 

Let us now consider the validity of the inverse 
of Eq. (32). If the relation in question were 
exact, we could, of course, invert it; but it is 
actually only approximate. We know from the 
derivation that if the left-hand member fits the 
function ¥(s) exactly at P;, the right-hand mem- 
ber fits the same function approximately at P2 
and along the neighboring portion of the axis of 
reals. The inverse relation would be the statement 
that if the right-hand member fits ¥(z) exactly at 
P2, the left-hand meu.ber is approximately 
correct at P;. With the aid of our partial determi- 
nation of the matrix |'¢!| we can test the validity 
of this statement. We accordingly assume as in 
the right-hand member of (32) that 8;= .8,=e'"’*. 
It follows that the corresponding values of a, and 
a; are not e‘*/4 and 0, but e‘*/4[ 1 —6*+-¢;;(6* —5) ] 
and (6—6*)e~‘*/*. Although the correction to the 
value of a; is small, its product by the dominant 
approximation function f/f; is not necessarily 
small. Hence the left-hand member of (32) may 
be entirely incorrect if the right-hand member is 
exact. 

The second connection formula is a direct 
consequence of (35) and our information regard- 
ing |'g'!. Consider an exact solution of (1) which 
has the form 


¥(x) =2p-! cos (w+y)=erfitef, 
in the neighborhood of P:. By (35) the corre- 


sponding values of a; and a, are, if we neglect 
small quantities of the order of ya, 


TT 
a;=2e'*!/4 cos (1-*) : 
+ 
T us 
amen 2x cos (v-*) —sin (+-=) | 
4 4 


Here we have introduced the symbol x, for the 
unknown, but finite, imaginary part of g;;. If y 
is not too close to the critical value — 7/4 (cf. the 


first connection formula) so that | cos (y—7/4) | 


is not much smaller than |sin (y—7/4)|, the 
product a,f, is small compared with a;f;. Hence 
we have the conneetion formula 
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T ; z 
cos (»-*) p| et S21 Pl p-! cos Bi pdr+y. (38) 
4 ev 2’ 
i Region I Region Il : 


This is usually specialized by writing y= 7/4, but 
the possibility of a generalization has been 
indicated by Langer (last reference, Note 3). 

In addition to the restriction on y, the only 
condition for the validity of (32) and (38) is that 
there shall exist a path A in the upper half plane 
connecting P; with P2, enclosing no complex 
zeros of p*? and having the property that wsa<1. 

Both connection formulas in the form given 
above presuppose that the slope of the real curve 
y= (sz) is negative at 2’. In order to obtain the 
connection formulas for the case in which V(z) 
has a positive slope at the classical turning point 
s’ it is only necessary to make the transformation 
——s in (32) and (38) and to interchange the 
roles of the regions I and II. 


u 


5. LimiTinG CASE WHERE ONE OF THE CON- 
NECTED Points Is AT INFINITY 


So far we have assumed that the points P; and 
Py» are located at finite distances from z’ and can 
be connected by a good path A of finite length. In 
applying the B. W. K. method to energy level 
problems, however, we frequently wish to assume 
that ¥(s) vanishes at s= +x. In this case it 
becomes necessary to discuss the changes in the 
coefficients a;, a, which occur in a path which 
extends to infinity along the axis of reals. 
Consider, for example, a case in which p(s) is 
imaginary and bounded from zero for all real 
values of z to the left of the turning point 2’ of 
Fig. 1. Let it be required to find the values of 8;, 
8, appropriate to P, when a; is known to vanish 
at s= — ©. The required transformation of the 
a's can be thought of as made in two steps, of 
which the first gives the values at P; in terms of 
those at — © and the second is the one already 
studied giving the values at P2 in terms of those 
at P. 

Consider first the value of a; at P;, which as 
before we shall call a;. It follows from Eq. (9) 
that 
Pi 


Of ds. (39) 


t 
a;=— 
9 


= _ 





As both f, and y approach zero exponentially for 
large negative values of z, a;| is clearly very 
small. We compute an upper limit. Let 2, 
denote the value of z at P;. Eq. (39) is equivalent 
to 


Hence 





w( 21) 
jai <3) my lett |p| dw]. 


"(—oo) | 2p! | 


Let m denote the maximum value of |Q/\*p! in 
the range of integration. This we may assume to 
be small. Since | ¥(s)| decreases monotonically as 
we pass from s; out to —~* along the axis of 
reals, 


m1 
| av; | <— view f 
2 w'( 24) 


m 
<—!¥(s)) eveGn! 
? 


e«ld|\ w!| 


ml 
<————_[ J ai| + |p em2 (2017, 39a) 
2 p(z1)?) 


If m,/ p(s:)}, is much less than unity, a, 
</la,|e~ vGn!, 

By analogy with the corresponding problem 
for a finite path one expects that a, is sensibly 
equal to the limiting value of a,(s) as s moves out 
to —x. It is not necessary, however, to in- 
vestigate this point in order to deal with the 
eigenvalue-eigenfunction problem. 

Consider the application of Eqs. (33) to the 
determination of 8; and 8, when a;(z) vanishes at 
the negative end of the axis of reals. In view of 
the upper bounds on a;! and | g,;, 

m 


! | 
rid; < MA| ArZor 
i 


2 P(21)?| 


Thus, by Eqs. (33), if P: is connected with P: 
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by a good path, and if |Q(s)/dMp(s)'p(s,)!! is 
small between P; and — ~, 8, has sensibly the 
same value as if a; were zero, i.e., 8, is sensibly 
equal to a,. We conclude that the connection 
formula (32) remains valid when the point P; is 
moved out to infinity along the negative axis of 


reals, provided that m/| p(s,)?| <1. 


6. EVALUATION OF THE LoW ENERGY LEVELS OF 
A ONE-DIMENSIONAL OSCILLATOR 


The usual wave mechanics derivation of the 
Sommerfeld “phase integral’’ quantum condition 
(with half-integral quantum numbers) is based 
on the use of the connection formulas and loses 
its validity when applied to the lower energy 
levels. Thus in the case of the Planck ideal linear 
oscillator, where (sz) = }k2°, we have 


Q=k(E+6V)/(E-V)?. (40) 


The minimum value of |Q/\?p"! consistent with 
the reality of p occurs at the origin and is 
k/(8us*). In the case of the lowest energy level 
the value is } and the corresponding minimum 
value of A" f°|Q/p!ds for any path joining the 
region of classical vibration with any point 
outside this region is 4. As this is of the order of 
unity, appreciable departures from the con- 
nection formulas are to be expected for this 
lowest energy level, and it is surprising to note 
that the value of & given by the phase integral 
formula is in exact agreement with the rigorous 
value obtained by other methods. Thus one is 
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led to inquire whether there is not some way of 
validating the Sommerfeld formula’ without 
assuming the existence of good paths joining the 
region of classical vibration with the regions of 
imaginary momentum on the axis of reals. 

Birkhoff® has recently sketched a proof of the 
phase integral formula applicable to the lower 
energy levels. Partly, however, due to the fact 
that Birkhoff's technique is unfamiliar to 
physicists, and partly because his proof rests on a 
plausible, but as yet unverified assumption, it 
seems worth while to give an independent 
derivation here. 

Let the problem under discussion be that of 
finding the discrete energy levels of a one- 
dimensional oscillator with the fundamental 
interval —*x <s<+~« and a single potential 
valley in which the classical local momentum 
p(s, £) takes on real values. We shall assume at 
first that I(s) has the special parabolic form 
|= }ks and consider later on the modification of 
the argument for an anharmonic oscillator. 

Let the complex s plane have a cut extending 
from the left-hand classical turning point 3’ 
through the right-hand turning point 3” to + 
along the axis of reals. The functions p, p' and w 
are uniquely defined over the cut plane by the 
requirement that they conform to Eqs. (28) for 
all negative real values of z—s’. They will then 
conform to Eqs. (29) on the upper lip of the cut 
between 3’ and 3”. For positive real values of 
s—s’’ we have 


Upper lip: p=—clp!; pi=er™4| p}]; w=rf lplav—af |plar: 


(41) 


Lower lip: p=—ilp{; pi=—e-'*4| p}]; w=—rf iplas—af \p|dé. 


We designate by III that portion of the region of 
imaginary momentum on the axis of reals to 
the right of s’’ where |Q/\*p?| <1 (cf. Fig. 4). It 
follows from Eqs. (41) that |e’| is very large in 
both of the regions I and III. The level lines of 

e’’| are readily sketched in and show that any 
path connecting I with III and avoiding the 
classical portion of the axis of reals must cross a 
single valley between the ridges on which it 





terminates. This is true whether the path lies in 
the upper or lower half plane. 

Let s; and z; be any two points on the axis of 
reals in I and III, respectively. It follows from 
(40) that in the case of the parabolic potential 
function under consideration it is always possible 
to join 3; and s; by a path A for which uas<1. In 
more general cases such a path will exist provided 


8G. D. Birkhoff, Bull. Am. Math. Soc. 39, 696 (1933). 
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Fic. 4. Level lines of |e*”| for harmonic oscillator. 
there are no complex zeros of V —E near to 2’ and 


Let us now assume that we have to do with an 
eigenvalue of the energy H=£, and a corre- 
sponding eigenfunction y,(s). Then y, vanishes 
at both ends of the axis of reals. Since f;(z) 
becomes infinite at both ends of this axis, a; 
must vanish there. It follows from our previous 
analysis that a; is very small at both of the points 
s; and s;. Let E denote the point where the path 
A crosses the bottom of the valley, i.e., the point 
on A where |e“! takes on its minimum value. 
Integration from zs, to E shows that a, is sensibly 
constant along the corresponding portion of A. 
Similarly integration from 3; to E shows that it is 
constant along the rest of A. If we reflect the path 
A in the axis of reals we get a second path 
connecting 3; with 23; along which a, is constant 
neglecting small quantities of the order of usa,(z;). 
The existence of these two paths shows that a, 
takes on the same values on the upper and lower 
edges of the cut at 3). y, is, of course, single-valued 
along the entire axis of reals. Hence f,(z) must 
have the same value on the upper and lower 
edges of the cut at z;. Equating the two values of 
f, in question and using the corresponding values 
of p' and w given in Eqs. (41), we obtain 


exp [ASL ?"| p| dg} = —exp [—ASy*”"| plds}, 
or the familiar equation 
£ pdz=(n+43)h. (42) 


The accuracy of the result is limited only by the 
quality of the path A. But in the case of the ideal 
linear oscillator it is possible to reduce the value 
of ua below any assignable quantity by choosing 
the path A as a semicircle of sufficiently large 
radius with the origin as its center. Hence the 


above equation is exact for the case under. 
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Turning now to the problem of the anharmonic 
oscillator, we observe that in order to secure 
single-valued functions f;, f, it is necessary to 
introduce cuts extending from the complex roots 
of E—V to infinity. The argument then goes 
through exactly as in the case of the harmonic 
oscillator except that the path A must lie inside 
all the roots above mentioned in order to avoid 
crossing the cuts. Consequently it is not possible 
to reduce the value of ua below a certain mini- 
mum determined by the location of the roots of 
E—V. The formula (42) is essentially inexact for 
anharmonic oscillators. 

The uncertainty 6£ in the energy due to the 
fact that ua and m do not actually vanish is 
readily calculated if we retain terms of the first 
order in these two infinitesimals. Let @,(z;) and 
a,(z3) denote the values of a;(z) and a,(z) on the 
lower edge of the cut at 23. Equating the approxi- 
mate expressions for ¥ on the two edges of the 
cut at s; in terms of fi, f, we obtain 


. a x * a Se 
a,(23) +e? & VPl dS (g3) 


AD) ae eee a a 


= ’ e 2 , 
By(Z3) +e? S11 Aha ,(z5) 


Let «, and é, denote the small quantities a,(z;) 
—ad,(%;), &(s3) —a,(z,), respectively. To a first 
approximation 


umr+rAJ) = 


€y —€, + [a (23) —ay(23) JerJeririat 
loge | 1+— ae... 


a y(2,) 





Introducing appropriate upper bounds for e,, é,, 
a;(z3)|, |G@(z3)' from the inequalities (27) and 
(39a), we obtain 


|a@+ASJ—2rk| <uatm/ | p(zs)'|, k=O, 1, 2, ---, 


where m is the maximum value of |Q/d*p!! 
between 3; and + «. Let AE denote the quantity 
hdE/dJ which measures the approximate spacing 
of the energy levels given by the Sommerfeld 
formula. Then 


bE 1 m 
-—(mt+—"_), 
AE 2r | p(z3)}| 


The writer has carried through a calculation of 
ua and m/|p(s3)'| for the Morse curve potential 
function corresponding to the normal state of the 
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II, molecule. For simplicity of computation the 
energy / was so chosen as to bring the classical 
turning points 2’ and z’’ together. A was taken 
to be a semicircle with radius equal to half the 
distance from zs’ to the nearest complex root of 
P(E, s). Under these circumstances yy is ap- 
proximately 0.015 and m/ p(s3)!, is 7.5X10~°. 
The corresponding value of 6£ AF is 0.0025, 
showing that the B. W. Kk. approximation should 
give the energy correct to at least two or three 
tenths of a percent of the spacing of adjacent 
levels. Of course the actual error involved should 
be appreciably less than this upper limit. 


7. HIGHER ORDER APPROXIMATIONS 


It is well known that our functions f;, f, can be 
obtained from (1) by making the transformation 
y=e%S v=, expanding y formally in powers of h, 
or 1/X, and taking only the first two terms of the 
series. By taking » terms of the series, one obtains 
functions f;, f,°” having properties similar to 
those of f; and f,. The series are semiconvergent 
and hence the f'”’s, regarded as approximate 
solutions of (1), at first improve in quality as n 
increases, and then grow worse. They obey 
equations of the form 

d? f( 

I a Uf mG, 

dz" 

where Q(z, E) is a polynomial in 1/\ beginning 
with a term of degree n—2. If Q=Q™ is small at 
any point of the complex plane, Q° will be 
smaller than |Q) provided that n >2, but not too 
large. Thus if |Q/p*! is small along a path A 
joining 2; and z3, |Q%/p*! will be still smaller. 
We infer that the argument used to prove that, 
in the case of an anharmonic oscillator and an 
eigenvalue of /, the function f, has sensibly the 
same values on the upper and lower edges of the 
cut in the region JJJ, is applicable with a higher 
degree of precision to f,“. It follows at once that 
the energy level condition® 


Sy Pdz=mh, m=0,1, 2, +: (43) 


is even more accurate than (42) for small values 
of m greater than 2. 


Cf. J. L. Dunham, Phys. Rev. 41, 713 (1932). Dun- 
ham’s derivation of (43) is fundamentally sound, though 
no more rigorous than Zwaan’s proof of the connection 
formulas. 
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8. MopIFICATION OF THE B. W. Kk. METHOD FOR 
THE RADIAL EQUATION IN THE 
Two-Bopy PROBLEM 


In applying the B. W. Kk. method to the radial 
equation of the two-particle problem one meets 
with an apparent difficulty in that the B. W. k. 
approximations do not have the right character 
to fit the exact of the differential 
equation at the left-hand boundary point r=0. 
As Kramers has pointed out, however, it is 
possible to fit the boundary conditions at both 
ends of the fundamental region if we modify the 
approximation formulas by the addition of the 
term h?/327°ur? to the potential energy. The 
added term is negligible except in the immediate 
neighborhood of the origin, and it is ordinarily 
possible to find a good path A leading from the 
origin around the classical turning points (there 
is only one when the angular momentum is zero) 
in the complex plane to a point on the axis of 
reals outside the region of classical motion where 
the approximation function is good. In this 
manner it is possible to apply to the two-particle 
problem essentially the same reasoning as we 
have used above for the anharmonic oscillator in 
the interval — x <s<+ 2. One obtains the 
energy level formula (42) with p(r, £) modified 
by the substitution of (/+3)* for /(/+1) in the 
term which gives the contribution of centrifugal 
force to the effective potential energy of the 
radial vibration. Further details regarding this 
application of our technique will be published 
elsewhere."® 


solutions 


9. THE TRANSMISSION OF PROGRESSIVE MATTER 
WAVES THROUGH A POTENTIAL BARRIER 


The transmission of matter waves through 
potential barriers has been the subject of much 
discussion,'! but so far as the writer has been able 
to discover there has been no application of the 
B. W. K. method to this problem that is properly 
applicable to cases in which the energy of the 
incident particles is not very different from the 
maximum potential energy of the barrier. 

Let us assume a rounded potential hill of 
approximately parabolic form. Consider first the 

10 In a book on the general principles of quantum me- 
chanics now in preparation. 


1 Cf.,eg., N. H. Frank and L. A. Young, Phys. Rev. 
38, 80 (1931). 
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case in which the top of the hill projects above 
the energy level of the incident particles. We 
divide the axis of reals of the complex z plane into 
three parts F, G, I/, separated by the classical 
turning points 2;, 22 as indicated in Fig. 5. 
Single-valued functions f;, f, are obtained by 
introducing a cut along the axis of reals from 2; 
to +e and by choosing those branches of p 
and p' which are negative real and negative 
imaginary, respectively, in F. p, p',and w are to be 
analytic over the cut plane. p is then real and 
positive in /7 and positive imaginary in G. On 
both sides of the hill the function 


2re . 
fe 2riEtih =p j exp | f pdé —-Ft | 
1 2’ 


describes an outgoing wave, while f,e~?™*?!/" 
describes an incoming wave. The level lines of 

e‘“ fora parabolic potential function are shown 
in the figure. f, is dominant in most of the first 
and third quadrants, while f; is dominant in most 
of the second and fourth quadrants. 

Let us assume that progressive waves are 
incident on the potential hill from the left side 
only. There will then be reflected as well as 
incident waves in F, but only transmitted waves 
in 77, We accordingly assume that a, vanishes in 
the first quadrant where f, is dominant. Let the 
regions F and J/ be connected by a good path A. 
Then a; will be sensibly constant along this path. 
Thus we have a connection formula of the type 


fitho—fi, (44) 


where c is a constant whose value has yet to be 
determined. As a matter of fact we can evaluate 
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only the modulus of c. This is done by using the 
constancy of the current density 


h dy  dy* 
Dtv}=—(v* ene 4 ) 


11 dz dz 





on the axis of reals. Choosing 2; as the origin of 
the integral w, we give e‘” the absolute value 1 in 
the region F and the absolute value e~* in J/ 
where K denotes the integral (27/h)J/-,*| p'd¢. 
The right-hand member of (40) now represents a 
current of density e°* in the direction of the 
positive s axis, while the left-hand member 
represents a current with the density cc*—1 in 
the same direction. Equating these two currents, 
we find that cc*=1+e~°4. The corresponding 
transmission coefficient is 


T =e~*K /¢c* =1 (1+e?*), (45) 


This formula is free from all restrictions depend- 
ing on the height of the hill provided that 
E— Vmax <0. Like Eq. (42) it is easily seen to 
hold for any analytic potential hill which yields 
no zeros of E— V' near 2, and 22 and hence permits 
a good path A joining the regions F and /1/. 

If the maximum potential energy is less than 
the energy of the incident particles, a parabolic 
potential function yields two imaginary roots of 
E—V(s),2= +¢|2(E—Vimax)/k}', where k denotes 
the negative of 1”. Denoting these roots by 2; 
and zz and the integral (27/h) S.,*:| p| |d¢| by K’, 
it is not difficult to prove that the transmission 
coefficient is 


T=1/(1+e-24’). (46) 


The reader will observe that the formulas (45) 
and (46) join continuously in the limiting case 
where 2,;=22=0, giving the common value 3 
for 7. 

The writer is greatly indebted to Dr. Eugene 
Feenberg, Professor G. D. Birkhoff and Professor 
J. H. Van Vleck for helpful suggestions. 
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Surface distortion and rupture permits field emission 
from liquid surfaces at field strengths less than those 
effective for equally smooth solid surfaces. An approximate 
mathematical theory of the rupture of a plane liquid sur- 
face in a uniform electric field has been developed. The 
relation between initial distortion, rupture time, and field 
strength has been calculated for fields large compared to 
that which just renders the surface unstable. This critical 
field is EL, =22'/*(pgT)!/4=53 kv cm™ for mercury, where 
p is density and T is surface tension. Typically, the theory 
shows that a hump initially 4X 1075 cm high and of diam- 
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eter 9X10~4 cm will lead to rupture in 5X10~° sec. in a 
field of 10° v cm™. Relative to initial humps in the surface 
whose linear dimensions vary inversely as the square of 
the field, the time to rupture varies inversely as the cube 
of the field strength. This calculation shows that a lowered 
sparking potential to liquid mercury can be ascribed to 
surface rupture and shows that it is possible that surface 
rupture plays a part in Beams’ low field emission from 
liquid mercury. Possible application of the theory to the 
high field condition at the cathode spot of the Hg arc is 
not clear. 











N the course of a study I have been making on 
the breakdown between a metal electrode and 
a liquid (Hg) surface close to it, I have observed 
that sparking occurs at considerably lower volt- 
ages just above the freezing point than below it. 
This suggested that the plane surface of the 
liquid might be deformed by the mechanical 
force accompanying the electric field, thereby 
concentrating the field, augmenting the disrup- 
tive force, finally building up a point sharp 
enough to give field emission and initiating the 
discharge. The same theory might apply to ex- 
periments performed by J. W. Beams! in which 
he applied 1 microsecond pulses of increasing 
voltage across an iron-mercury vacuum gap. 
Whereas theory requires 3107 v/cm for ap- 
preciable field emission, he found that an arc 
started at 1.8X10° v/cm with clean mercury. 
The duration of the voltage pulse was deemed 
too short for the liquid to respond mechanically, 
but no analysis supporting this view was offered. 
According to the following approximate treat- 
ment, a field of 1.8 10° v/cm™ would suffice to 

1J. W. Beams, Phys. Rev. 44, 803 (1933). In a very 
recent letter J. W. Beams describes tests performed by 
Mr. Quarles which indicate that in the time duration and 
field strength range where he was working surface rupture 
of the kind discussed here was not occurring. It is too early 


to say just what the relation is between those experiments 


and the present theory. 
Note added in proof: Since submission of this paper, an 
article by L. R. Quarles, Phys. Rev. 48, 260 (1935), 


describes experiments showing that at the high fields 
emission did not depend upon a time response of the sur- 
face. In a private communication, however, Professor 
Beams states that the sparking voltage did begin to de- 
crease markedly for application times of several micro- 
seconds, so that it would be there that the theory developed 
in the present paper should find application. 
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give field emission in 10~ sec., so that it is quite 
possible that point formation on the liquid was 
responsible for initiating the arc. 

Furthermore, as Dr. Langmuir has suggested, 
peak formation may be occurring in the intense 
field at the cathode spot of the Hg are. 

The bursting condition, -?>1677 R, for a 
liquid droplet of radius R and surface tension 
T subjected to a radial surface field FE is well 
known, and the distortion and disruption of 
bubbles has been frequently investigated.? The 
present problem differs essentially from the 
older one in that the stability of a surface of 
zero curvature is to be investigated and we are 
interested in the development of distortion as a 
function of time. 

The problem can be considered in two parts; 
first, what is the relation between applied field, 
surface distortion and resulting net pressure on 
the surface; and second, what will the motion of 
the liquid be under the influence of this pressure ? 

The rigorous solution is impossibly com- 
plicated, so that an approximate treatment will 
be followed from the very start. Consider an 
axially symmetrical hump in the liquid surface 
as indicated in Fig. 1. Let us use the three quanti- 
ties, R, y and r to characterize it, R being the 
radius of curvature at the vertex of the hump, y 
the height of the vertex above the undisturbed 

2C. T. R. Wilson and G. I. Taylor, Proc. Camb. Phil. 
Soc. 22, 728 (1925); J. J. Nolan, Proc. Roy. Irish Acad. 37, 
28 (1926); W. A. Macky, Proc. Camb. Phil. Soc. 26, 421 
(1930); Proc. Roy. Soc. A133, 565 (1931); J. J. Nolan and 
J. G. O’Keeffe, Proc. Roy. Irish Acad. 40, 86 (1932); G. D. 


West, Proc. Phys. Soc. 44, 336 (1932); J. Zeleny, Phys 
Rev. 16, 102 (1920); 47, 638 (1935). 

















LIQUID SURFACE 


Fic. 1. Shape of hump assumed in liquid surface. 


level, and 2r the effective diameter of the hump, 
taken here as the diameter of the small circle 
cut from the y=0 plane by the sphere R. By 
geometry 


yR= 3(? +"). (1) 


We now analyze the component pressures on 
the vertex of the hump. The impressed uniform 
field will be denoted by Eo and the enhanced 
field at the vertex by E. The pressure trans- 
mitted to the vertex from the undisturbed por- 
tion of the liquid surface is made up of three 
components, (1) the pressure A of any atmos- 
phere above the liquid, (2) the hydrostatic 
pressure —pgy, p being the density of the liquid, 
and (3) the electric pressure — E,?/87, with a 
total value for the three of A —pgy—E,?/8r. 
This is the upward pressure exerted against the 
hump vertex by the liquid beneath it. The down- 
ward pressure there also has three components; 
(1) that of the atmosphere, (2) that arising from 
surface tension, 27/R, where T is the surface 
tension of the liquid, and (3) the electric pressure, 
—E?/8x, totaling A+27/R—E?*/8xr. Subtract- 
ing the latter from the former gives the net 
upward pressure P at the vertex: 


P=(E?—E,?)/84r—pgy—2T/R. (2) 


If P is positive for any assumed conditions, then 
the surface tends to greater distortion; if nega- 
tive, it tends to flatten out. 

To relate E to y and r we note that for y=r, 
that is, for a hemispherical boss on a plane, the 
field is three times the normal. Accordingly, we 
assume the linear relation 


E=E,(1+2y/r) (3) 


which fits the two values y/r=0 and y/r=1 and 
lies sufficiently close to the theoretical values for 
our treatment. 
Substitution into Eq. (2) of Eq. (3) to elimi- 
nate EF, and of Eq. (1) to eliminate R gives 
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Fic. 2. Stability of surface relative to impressed field Fy 
and hump radius r. Exponent of (pg7°) should be 4. 


‘ytr 4nT/E? xpg 
i 

r- y 4+ r? E?? 
This equation shows that we are dealing with a 
disturbance that has no equilibrium condition. 
When the bracket is positive, y tends to increase, 
which in turn increases the value of the bracket. 
When the bracket is negative and y decreases as 
a result, the bracket becomes still more negative. 
We must suppose that the disturbance arises 
as a continuous transition, so that in investigat- 
ing the incipient distortion, y<r, and we may 
neglect y in the brackets. Thus we find that the 
growth or subsidence of a small distortion ensues 

according as 


1/r>or< (4947 /E,*)/r?+rpg/E? 
or, in terms of the field, as 
E?>or <4xT/r+xper. (5) 


The stable and unstable regions in the Eor plane 
are mapped in Fig. 2. The minimum value E,, of 
Eo which will cause instability and the corre- 
sponding value r,, of r are seen to be 


Em=20*(pgT)', tm=2(T/pg)!. (6) 


Taking 7=450 dynes cm and p= 13.6 g cm" 
for mercury and g= 980 cm sec.~*, the numerical 
values are 52,700 v/cm™ (176 e.s.u.) for E£,, and 
0.367 cm for r,,. The interpretation to be made of 
the plot is that on a large free surface subjected 
to a given Eo which is greater than E,, only 
those distortions will grow those radii lie within 
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the limits indicated. This suggests an experi- 
mental test in which ¢ is limited to values less 
than that given by Eq. (6). 

The answer to the first portion of the problem 
is now seen to be: Below a certain field strength, 
F,, no distortion of the surface arises, but above 
that field strength the distortion increases pro- 
gressively under the net unbalanced pressure P. 

In particular, it can easily be shown that for 
r=r,, the bracket in Eq. (4) becomes 


(1—(E,, Eo)*) 2A, (7) 
where A=8r7'/F,?. As the (-,, /o)* term repre- 
sents the gravitational force, it is evident that 
when /,'>E,4 gravitation can be neglected, a 
condition which is sufficiently well satisfied 
when Eo >3En. 

As the hump develops, liquid flows into it from 
its immediate neighborhood. On the average, this 
tlow will extend over a distance comparable with 
r, so that associated with a surface area a at the 
vertex of the hump, a volume of liquid ar is in 
motion with a velocity y. The corresponding 
momentum is thus pary. Its rate of increase is 
equal to the force acting, whence 


aP=pad(ry) dt. 


By substituting from Eq. (4), we find for the 
equation of motion 


rt+y A/2_ «pg 
d(ry) /dt= By(~ —— = *) » (8) 
Pr Pty Ee 


where B?=E,?/rp. If we confine attention to 
fields greater than 3£,,(=0.158 MV/cm for 
mercury) the gravitational term can be neg- 
lected in accordance with Ea. (7). If, further, new 


variables 





7=(B/A)t, w=y/A,o=r,A (9) 
be substituted, Eq. (8) becomes 
o+w 1,2 
d(ois)/dr=w( -— ) (10) 
Pe o+w 


where, as later, the dot represents the derivative 
with respect to the time-like variable. 

To solve this equation, the relation between 
o (r) and w (y) must be known. It seems reason- 
able to assume that r varies in such manner as to 
give the right member of Eq. (8), and hence of 
Eq. (10), the maximum value. In other words, 
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Fic. 3. Method of determining relation between hump 
height (~w’) and radius (~o). 
the development of the disturbance “follows the 


path of least resistance.” 
The routine analysis for this condition leads 
to the parametric equations 


(1+A7)?(2+A), o=Aw. 


»— }4 
w= 


(11) 


It is certainly to be expected that as the deforma- 
tion progresses, the volume of the hump shall not 
decrease. But plotting ow, as a measure of the 
volume, against w, these equations lead to curve 
A of Fig. 3. We must, accordingly, modify our 
original assumption so as to have it apply only 
to the initial stage of deformation. Subsequent to 
this we postulate that the volumetric change ap- 
proaches zero as y (w) approaches r (¢). This 
anticipates the second stage in the surface 
deformation which begins when y (w) exceeds 
r (o), and throughout which it will be assumed 
that the hump volume is constant. The extra- 
polation of the beginning of curve A, Fig. 3 
in accordance with this change in assumption was 
carried out by plotting first differences along 
A up to w=0.08, and from there fairing in a 
first difference curve to cross the axis where 
o=w, which occurs at about w=0.45. Curve B 
is the continuation of curve A based on this 
method. The values of o which have been 
calculated from curve B are shown plotted 
against w in Fig. 4. This made it possible to 
plot f(w), the right member of Eq. (10), as a 
function of w in Fig. 5. 

We can now proceed to solve Eq. (10). In so 
doing we shall confine ourselves to that solution 
whose initial condition is w=0 and w=0 at 
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Fic. 4. Hump radius (~¢) as a function of height (~z). 





























Fic. 5. The function f(w). 


7=-— x. the magnitude of the time interval for 
disruption will not depend markedly enough on 
the initial relation of velocity to displacement 
to warrant carrying out a number of numerical 
integrations for different initial conditions. 

For small values of w the solution 


W= wy exp (7/23) (12) 


is exact to the degree that 2w» is negligible com- 
pared to unity. This equation therefore serves to 
establish the initial conditions, w= wy, w= wy ‘2! 
at r=0, for the numerical integration, and to 
permit extrapolation to find the times involved 
when the initial disturbance is less than w. 

For the numerical integration the value 


Wy = 0.04, so that wo= 0.0283 (13) 


was chosen, and the calculation was carried out 
by well-known methods.’ Fig. 6 shows the height 
of the hump as a function of time. Denoting the 
values of w, o and r when w=c by the subscript 


3R. E. Root, Mathematics of Engineering (1927), § 74. 
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Fic. 6. Hump height as a function of time. 


1, it was found that w,;=o0,= 0.450, 7;= 2.41, and 
w= 0.431. These are the initial conditions for the 
second phase of the rupture in which no new 
liquid enters the hump but the concentration of 
force at the vertex creates a rapid peaking. 
Thus, at the beginning of the second stage we 
have for y, R, r, and y: 


yi= Ri =r, =W,A = 0.450A 
t}=7,A B=2.41A ‘B, (14) 
dy, dt=w,B=0.431B. 


In this second stage we shall suppose that 
the shape of the hump near the vertex is that of 
a cylinder terminated by a sphere, as shown in 
Fig. 7, and that the growth of the peak is by the 
upward flow of liquid from the cylindrical zone 
of altitude, R 2, just below the base of the 
hemisphere to form an additional segment of 
narrower cylinder capped by a smaller hemi- 
sphere. Again this is a rough approximation, for 
the assumption itself leads to a cone frustum 
rather than a cylinder below the hemisphere. 











Fic. 7. Hump vertex assumption in second stage. 
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Fic. 8. Evolution of hump. 


Of more importance will be the neglect of 
differences in velocity at different points of the 
hump through treating the motion as that of a 
rigid body of variable mass when the equation of 
motion is written. It appears probable that taking 
this factor into account would explain the 
thread and ultimately the drop formation which 
has been observed in the bubble and droplet 
disruption experiments. The present hypothesis, 
however, leads to an infinitely sharp peak when 
y reaches 3y;, as shown by Eq. (24), and it is 
this condition which shall be deemed to be 
surface rupture. 

Incidentally, this relation permits the delinea- 
tion of the sequence of hump shapes from the 
initial small distortion of stage one through to 
rupture itself. This is illustrated in Fig. 8 where 


no attempt has been made to show regular time. 


intervals. 

Although the postulated shape of the hump 
peak during the second stage is that of a hemi- 
sphere upon a cylinder, the configuration left 
by the termination of the first stage is a 
hemisphere upon a plane. Yet it does not seem 
worth while to devise a special treatment for 
this transitional phase, particularly as the whole 
duration of stage two is but a fraction of that of 
stage one. 

Equating of volumes above the PP plane of 


Fig. 7 before and after deformation gives 
dy= —2dR, (15) 


dR being intrinsically negative. The volume is 
(5/6)7R° and its center of gravity moves upward 
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1G. 9. Method of calculating force on the vertex. 


—(3/2)dR, so that its momentum is 


—(3/2)dR(5 6)rR*(y dy)p=(5/8)rpR°y. (16) 


To calculate the force increasing this mo- 
mentum we shall treat the hump as an ellipsoid 
of semi-major axis y and radius of curvature at 
the end of this axis R. The force will then be 
taken as the resultant acting on the portion of 
the equivalent ellipsoid within a circle of radius 
R on its surface, as indicated in Fig. 9. 

The semi-minor axis, 5, of the ellipse and the 
eccentricity e are given, respectively, by 

b?= yR, e& =1-—R/y. (17) 
The charge density C at the vertex is 
Eo e° 
C=— —— — — -—— ——, (18) 
4m (1—e?) {(1/2) In [(1+e)/(1—e) ]—e} 


Elsewhere it is C sin y. The vertical component of 
force acting on the surface element ds is 


dF=2nC? sin® yds=2nC? sin? ydo. (19) 


Using r to denote distance from the hump axis, 
we have 
(1—e*)[y°(1 —e*) —1°*] 


(20) 
[y(1 —e?)?+ e*r? |] 





sin? y= 


Integration of dF from r=0 to r=R after the 
substitution of Eq. (20) into (19) gives 


E,*y* In (1+?) —e?(1 —e?) 
F= e . (21) 
8 {(1/2) In [(1+e)/(1—e)]—e}? 
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By equating this to the rate of change of the 
momentum given by Eq. (16), we obtain the 
equation of motion: 

d( R°y) /dt= (1/5) B*y*g(e), (22) 


In (1+?) —e?(1 —e?) 
g(e) = ee —$_$_$_—$—$—______—_—_——_.. (23) 
{(1/2) In [(1+e)/(1—e) ]—e}? 
The relation between y and R is given by inte- 
grating Eq. (15) and using the condition im- 
posed by Eq. (14). It is 


2R=3y-—y. (24) 


where 


Transforming variables by substituting 
n=9/y1, T= (t—t)(8/5)'B/y1 (25) 


and using Eq. (24) to eliminate R from Eq. (21) 
we have 
d(3—n)*n/dT=n?*g(e) (26) 


with »=1 at 7=0. From Eqs. (17) and (24) we 
find 

n= 3/(3—2e?’). (27) 
Thus, by calculating 7 and g(e) for a sequence of 
values of e from 0 to 1, values of the right mem- 
ber of Eq. (26) can be tabulated as a function of 
n, and the equation can be solved by numerical 
integration. 

Such a tabulation appears in Table I and a 
plot is shown in Fig. 10. It is seen that at 
n= 2.51 a tenfold increase of field is attained, and 
that up to this value of n, n? g(e) is fairly constant 
at the value 5. By using this value in Eq. (26), we 
can integrate directly to obtain an approximate 


solution 
(3—n)*7=57T+<¢, (28) 


where c is a constant of integration. Now at 7=0, 
we have n= 1, and 
"1 = dni. ldT= (dy, ‘y,) /dt(dT /dt) 

= (dy,/dt)(dt/yidT) =w,(5/8)* (29) 














TABLE I. 

e y/R E/E n n’g(e) T 
0 1.00 3.00 1.00 4.50 0 
0.4 1.19 3.24 1.12 5.10 0.257 
0.6 1.56 3.63 1.32 5.44 .504 
0.8 2.78 4.78 1.75 5.81 .746 
0.9 5.26 6.70 2.17 §.12 .814 
0.95 10.3 10.02 2.51 4.03 .829 
0.98 25.2 18.1 2.78 2.81 _- 
1.00 x“ x 3.00 0 .833 
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Fic. 10. The force function 7%g(e) of Eq. (26) and elapsed 
time (~7T) in second stage. 

















from Eqs. (14) and (25). Using these conditions 
gives for the first order equation 


(3 — )*9=5T+2w,(10)!. 


_ Another integration, adjusted to initial condi- 


tions, gives 
107°+8(10)!w,7 —16+(3—n)*=0. (30) 


Values of T for different values of 7 are given 
in Table I and shown in Fig. 10. Whether we 
seek a sixfold or greater increase in field makes 
but little difference in the value of T which 
rapidly approaches its final value 7;=0.83. 

In converting the mathematical quantities 
w, o, tr and T into the physical quantities, y, 7 
and ¢, we note from Eqs. (9) and (25) that the 
quantities A and B are of primary importance. 
Expressing field strength in megavolts cm™, 
denoted by Exyy, we have for mercury 


A=8rT/E?=1.02 X10“Emy~™, (31) 


B= E,/(xp)'=509Emy. (32) 
We then obtain directly: 
Yo= WA = 4.08 X10 Ey, (33) 
rp=0.913A = 0.931 X10 Ey, (34) 
_ (y1=n=0.450A = 0.459 X10 Fwy, — (35) 
ie dy; /dt=0.431B=229Ewy, (36) 
(23) 


t= 71A/B=4.8X10-8Eyy, (37) 
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Eq. | te—t:=(5 8)? 0.450 X0.834,B 
(24) = 0.5910 7 aN .. (38) 
and it follows that 

to= 5.4 x 10 $ Euy. (39) 


Thus an initial disturbance of height 410% 
cm and diameter 210-* cm in a field of 1 
megavolt cm~! will give rise to surface rupture in 
5.4u sec. The increase in time attendant upon 
smaller initial unevenness y, is found from Eq. 
(12) to be 


to= 2.83 x 10 SFausy 3 In (Vo ‘Va), (40) 


so that at 1 megavolt cm™ again, the rupture 
time is increased 2.8u sec. for every power of « 
that ya is less than yp. 

The striking feature of the result of this 
calculation is the extreme shortness of the time 
in which a marked mechanical response is 
achieved. To be sure, the peak which is raised is 
very small, being only 3y,=1.4X10“%Eyy-? cm 
high at rupture. 

In the experiments on sparking between mer- 


cury and a solid electrode the mercury surface . 


was so constrained by a meniscus that no short 
circuiting occurred when the sparking voltage 
was applied. Under these conditions there was 
a concentration of field at the mercury meniscus, 
and there was also an initial distortion of the 
surface. The field was produced by an Mo strip, 
edge down, 0.025 cm above the mercury. A 
60-cycle a.c. voltage of somewhat over 11,000 v 
peak was required for sparking at room tempera- 
ture. This rose about 20 percent as the mercury 
was cooled to its freezing point and there almost 
doubled when solidification occurred. At room 
temperature the sparking voltage rose 60 to 70 
percent when a second Mo strip was substituted 
for the liquid mercury surface, so that the spark 
occurred between two solid electrodes. Thus the 
effect of lowering the sparking voltage by sub- 
stituting a liquid for a solid electrode is definite. 
But the geometry was too indefinite for a 
quantitative test of the theory. Suffice it to 
point out that the field of the order of 0.4 MV 
cm! gives a value of 80u sec. to f2, which is very 
short compared to the duration of 1/2 cycle. 


TONKS 


In Beams’ experiment a field of 1.8 MV cm7! 
applied to a clean mercury surface caused elec- 
trical breakdown in 1 sec. Eqs. (33) and (39) 
state that for this field an initial amplitude of 
1.3X10°-° cm will give rupture in 0.93u sec. 
It is impossible to say whether or not dis- 
turbances of this magnitude were present. 
The corresponding diameter, 5.7X10-* cm 
by Eq. (34) is the wave length \ of capillary 
waves of frequency (277)/p)'A-?=10° sec.-', a 
frequency usually considered as outside the range 
of mechanical vibrations. Certainly, experiments 
to determine the magnitude of the waves on the 
mercury surface and the dependence of electrical 
breakdown on pulse duration can be performed. 

The possibility suggests itself that thermal 
roughness of the surface may account for the 
incipient rupture. R. Gans has calculated* the 
mean square amplitude by means of a Fourier 
analysis. The order of magnitude of the displace- 
ment for a hump is (k@/7)* x independent of 
diameter, @ being the absolute temperature. At 
room temperature this is 3X10-° cm. By Eq. 
(40), to~30u sec. at 1 MV cm”. Further con- 
sideration shows that the thermal velocity 
rather than displacement is the pertinent 
quantity, and it seems possible at this writing 
that this may furnish ‘‘strong’’ enough nuclei. 
Nuclei of optimum diameter, which have been 
the only ones considered here, are still too 
“weak,” but smaller nuclei may gain in strength 
(because of the operation of statistical laws) 
more than enough to compensate for what they 
lose in growth acceleration by the field. 

The application of the results of the present 
calculation to the cathode spot is less clear. The 
cathode drop of 10 v in a distance of s cm, the 
cathode sheath thickness, would give a_ field 
strength Eyy~10~ s. The height of peak at 
rupture would then be 3y,~1.410's? by Eq. 
(35), and for the hump to lie entirely within the 
sheath at this stage 1 >3y,/s~1.410"s, whence 
s<7X10-° cm, a value far below that at which 
field emission would occur without peaking and 
also far below the 10-° cm which has been calcu- 
lated in other ways. 


*R. Gans, Ann. d. Physik 74, 231 (1934). 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important 
discoveries in physics may be secured by addressing 
them to this department. Closing dates for this 
department are, for the first issue of the month, the 


The Photographic Absorption Spectrum of Gaseous 
Ammonia 


The absorption spectrum of the ammonia molecule, 
which is of the symmetrical top type, has been investigated 
by many authors. The four bands in the photographic 
region, 10,230A, 7920A, 6470A and 5490A,! have generally, 
been considered to be harmonics of the fundamental band 
at 3u(v:). Because of the complexity of the structure, 
only partial analyses of these bands have been made. 
A higher resolution has been attained in the present work, 
but a satisfactory analysis has not yet been completed. 
Nevertheless a brief preliminary report will be given in 
the following. 

The 10,230A band appears in the present work to 
consist of two separate bands which can be distinguished 
by their intensity in absorption and which slightly overlap 
at 10,150A. The one of longer wave-length is the stronger 
and its center is approximately at 10,230A, while the 
center of the weaker band appears to be at 9900A. In the 
former, there are several groups of strong lines on the high 
frequency side with intervals of about 10 cm™', half the 
spacing of the parallel tvpe band. This suggests similarity 
to the perpendicular type band such as 1.974, studied by 
Stinchcomb and Barker.* The structure of the weaker 
band is relatively simpler. 

In the 6470A band, a group of strong lines which are 
nearly equally spaced in the center of the band suggests 
the Q branch of a perpendicular type band, while there are 
also other groups of lines resembling P and R branches. 
Its appearance is similar to that expected for a perpen- 
dicular type band derived from rotation energy levels 


E,/he = BJ(J+1)+8BR°?, 


where B=9.92, BB=1.4 cm™ (1.4=half the average 
spacing of the lines in the Q group). There are, however, 
several discrepancies. For instance, the spacing of the 
lines due to changes of A with fixed J-values are smaller 
in the RX branch and are larger in the P branch than those 
in the Q branch. 

Recently Adei® proposed that the 6470A band is the 
fourth harmonic of the parallel type band, », with a 
doubling similar to the 10u band. But the structure of 
the band does not agree with this type, even with different 
values of B and @ in the normal and excited states. As 
already mentioned, the structure looks much more like 
that of a perpendicular type band, although it has not 
vet been possible to make it fit this type entirely. 

In order to account for perpendicular type bands in the 
photographic region, the following possibilities may be 
considered. The first possibility is that, if the 1.974 band 
is considered, as was done by Barker,? to be one of the 


twentieth of the preceding month; for the second 
issue, the fifth of the month. The Board of Editors 
does not hold itself responsible for the opinions ex- 
pressed by the correspondents. 


fundamentals of the perpendicular type, v2, then the 
10,230A band or else the 9900A band may be the first 
and the 6470A band the second harmonic of this funda- 
mental. The second possibility is that, if v2 is in the 
neighborhood of 3300 cm™ as suggested by Amaldi and 
Placzek,‘ then the four bands will be harmonics of this 
perpendicular type fundamental band. Although this 
3300 frequency has been found only in the Raman spectrum 
of liquid ammonia, it should be noted that Howard® 
recently showed by calculations based on the valence force 
potential function that one of the perpendicular type 
fundamentals is probably approximately at 3450 cm”, 
which is very close to »;. 

It has been shown’ that the harmonics of perpendicular 
type bands in general include a set of several parallel and 
perpendicular type bands, so that their structure should 
be rather complex. 

I wish to express my appreciation to Professor R. S. 
Mulliken for the proposal of the problem and for constant 
inspiration and valuable advice throughout the investi- 
gation. 

Sic-HunG CuHao 

Ryerson Physical Laboratory, 

University of Chicago, 
August 23, 1935. 
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Fine Structure in the K X-Ray Edge of Gallium 


Previous investigations have shown no extended fine 
structure in the K x-ray absorption edge of Ga. In view 
of Hanawalt’s work! on the fine structure in the K edge 
of Fe as a function of temperature, it was thought that a 
structure might exist in the Ga edge at temperatures suf- 
ficiently below the melting point, 30°C. Accordingly, 
absorption spectra in the vicinity of the Ga K edge were 
investigated at low temperatures, and a preliminary report 
of the results is given herewith. 

The Ga absorber was a rolled foil 0.001 inch thick, sup- 
ported on stiff paper. The source of x-rays was a tungsten 
target, water-cooled, metal tube run at 38 kv and 24 ma. 
A Société Génevoise double crystal spectrometer was used 
as a single crystal instrument, observations being made by 
the ionization method. Blank runs were taken to ensure 
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that the results were not spurious. The foil was cooled by 
a stream of vapor from boiling liquid air. 

At the lowest temperature, —150°C, five peaks were 
found at energies greater than the drop at the main Ga K 
edge by 77, 127, 154, 201, and 248 volts, respectively. The 
peak at 127 volts appeared as an unresolved shoulder on 
the 154-volt peak and the 248-volt peak is weak and 
unverified. The other three peaks are quite strong and 
definite. 

A single run through the structure at —67°C showed 
only peaks at 63 and 123 volts from the main edge, both 
of which were broader and less intense than the corre- 
sponding peaks at the lower temperature. Runs at room 
temperature agree with previous observations in showing 
no structure. 

The crystal structure of Ga is unique, and an attempt 
is to be made to compare the experimental results with 
Kronig's theory for the Ga structure. A full report of the 
work will be published soon. 

A value of the wave-length of the Ga K edge of 1192.5 
X.U. is found, agreeing with the value given by Kievit 
and Lindsay? (1192.9 X.U.) rather than the accepted value 
of 1190.2 X.U. 

W. W. Mvutcu 

Sloane Physics Laboratory, 

Yale University, 
August 26, 1935. 


J. D. Hanawalt, Zeits. f. Physik 70, 293 (1931). 
*B. Kievit and G. A. Lindsay, Phys. Rev. 36, 648 (1930). 





On the Recombination of Neutrons and Protons 


The process of the photoelectric disintegration of the 
deuteron, and the reverse process (recombination of a 
neutron and a proton with formation of a deuteron and 
emission of a y quantum) have been treated theoretically 
by Bethe and Peierls.' These authors consider the emission 
or absorption of y-rays by the ordinary mechanism of the 
electric dipole radiation. For the recombination cross 
section they obtain an expression which vanishes in the 
limit of low relative velocities of the neutron and the 
proton. This last result seems to be in contradiction with 
general experimental evidence, that slow neutrons are 
captured in hydrogenated materials such as paraffin or 
water. Recently the lifetime of a slow neutron in paraffin 
has been evaluated by a direct experiment and found to 
be of the order of magnitude of 10~* second.* 

I wish to point out that if we take into account also the 
radiation processes due to oscillations of the magnetic 
dipole of the neutron-proton system, one finds a proba- 
bility of capture which is in satisfactory agreement with 
the above value. The theory is based on the customary 
assumption that the fundamental state of the deuteron 
is a °S (neutron and proton with parallel spins). 

The ordinary theory of impact shows that in the case 
of slow neutrons only the capture from S states in the con- 
tinuum is of importance. The usual radiative processes of 
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electric dipole and quadripole are forbidden for these 
S—S transitions. Instead transitions from a 'S state in 
the continuum to the fundamental 4S state can occur by 
the mechanism of magnetic dipole radiation. 

The ordinary electromagnetic theory gives for the 
inverse mean life 1/7 of a neutron with velocity v moving 
in a medium containing protons per unit volume the 
following formula: 

1 64 rth 


=n-——— uo Bp — 8n)?| S-fer*dr|?. 
7 omy 





m is the neutron mass: v is the emitted frequency; jo is 
the nuclear magneton; gp and g, are the nuclear g-factors 
for proton and neutron. The last integral contains the 
normalized eigenfunction f(r) of the fundamental state 
and the continuum eigenfunction g(r). This last is normal- 
ized in such a way that for large r it goes over into (1/r) 
sin (2rpr/h-+const.). 

This last integral can be evaluated numerically by 
methods similar to those used by Bethe and Peierls. The 
parameter needed for the description of the 1S eigenfunc- 
tion can be obtained on the assumption that the high 
elastic scattering cross section of slow neutrons in hydrogen 
is mainly due to the 'S eigenfunction (real or virtual 'S 
level close to energy zero). We find for slow neutrons 


1 _ 16731? po? (gp — Bn)? 
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where / is the mean free path for the elastic scattering of 
slow neutrons (for paraffin about 0.5 cm)’, and W'=2.1 
X 10° e.v. is the binding energy of the deuteron. Assuming 
3 and 1, nuclear magnetons, as magnetic moments of the 
proton and the deuteron,‘ we get g,=6 and g, = —4. This 
gives 1/r=5.2 108 in agreement with the experimental 
result. 

The magnetic dipole radiation affects also the process of 
the photoelectric disintegration of the deuteron. The cross 
section for the photoelectric absorption of a y-ray by a 
deuteron results on the same assumptions 
(y—1)! 


47° a7 a7 u07( gp — gn)? 
Eee 
. She 








y{1+a%a?(y—1)}' 





' — , +2°mW hy 
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The cross section must be added to that already cal- 
culated by Bethe and Peierls in order to give the total 


cross section. 
Details of this theory will be published elsewhere. 
Enrico FERMI 


University of Michigan, 
August 27, 1935. 


1H. Bethe and R. Peierls, Proc. Roy. Soc. A148, 146 (1935). 
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40. Frisch and O. Stern, Zeits. f. Physik 85, 4 (1933); Estermann 
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Alpha Radioactivity of Argon Formed by Radiochlorine 


During a study of the beta-radiation from a sample of 
silver radiochloride by means of the screen-wall counters 
previously used in this laboratory! it was found possible to 
follow the decay of the radiochlorine when the sample was 
at a distance from the counter so that none of the beta- 
radiation could have reached the counter itself. However, 
the sample and counter were both enclosed in a chamber 
filled with argon at 2 or 3 cm of Hg pressure so that a gas 
evolved by the sample might have reached the counter. 
Following careful establishment of these observations 
direct evidence for the properties of the gas was obtained 
by passing gas which had stood for several minutes in 
contact with the hot silver radiochloride sample at a rate 
of a fraction of a liter per minute through a copper coil 
three or four feet in length immersed in liquid air. This 
gas was used to fill a solid wall tube counter which was 
immediately placed in operation. A fall in rate was observed 
which corresponded to an activity of half-life 1.9+0.3 
minutes. The initial effect was generally about fifty or 
seventy-five per minute superimposed on a background of 
fifty per minute. Changing the size and material con- 
struction of the counter had no appreciable effect on the 
behavior of the gas. From these results the conclusions 
were drawn that, first, the active gas must be formed by 
radiochlorine when it disintegrates and, secondly, that the 
gas must be argon because no other gas which radiochlorine 
might form would pass the liquid air trap used. 

The character of the radiation from the gas was partially 
established by covering the screen-wall counter used in 
the first experiment with an aluminum foil of approxi- 
mately 3 cm air equivalent and repeating the experiment 
in which the counts due to the diffused radioactive gas 
were used to follow the decay of the radiochlorine. The 
result was almost complete disappearance of the effect, 
indicating that the radiation from the gas had a range less 
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than 3 cm of air and was therefore probably alpha in 
character. 

To investigate further the gaseous radiation, a low 
pressure cloud chamber operating at 1/8 atmosphere of 
air and water vapor was used. The gas was passed from 
the warm silver radiochloride sample into the chamber and 
the chamber expanded at the rate of 20 expansions per 
minute. A polonium source was frequently placed near a 
mica window to furnish comparison alpha-particle tracks. 
In the course of the investigation of approximately two 
dozen silver radiochloride samples we obtained between 
50 and 75 tracks with both ends located in the gas and 
approximately 1.5 to 2.5 cm in length. All the tracks were 
observed within 8 minutes after transference of the gas. 
The densities of the tracks were comparable with those of 
the polonium alpha-tracks and definitely greater than 
those of beta-tracks formed by recoil electrons from a 
strong gamma-radiation source. 

Fig. 1 shows a picture of the direct and two mirror 
images of the chamber and one of the tracks. The pointer 
shown was being used to mark the position of a track in 
the preceding picture when the track shown occurred in 
almost exactly the same position. Reprojection of the 
three images showed the track to have had its upper end 
6 mm from the top of the chamber and its lower end 12 mm 
from the bottom, to have been inclined at a vertical angle 
of 49.5° and to have been 22 mm long. Reduction to 
standard air gives 3.0 mm as the equivalent length of this 
track. 

Data on four tracks gave standard air ranges of 2.9, 3.0, 
3.1 and 3.2 mm, indicating a probable value of 3.1+0.2 mm 
in rather close agreement with the prediction of 1 to 3 mm 
of the Gamow, Gurney and Condon theory for an alpha- 
particle from argon of 1.9 minutes half-life. 

The silver radiochloride was prepared by precipitation 
from the water phase resulting from the extraction of 





Fic. 1. Sample picture of track. Arrows indicate positions of the direct and two mirror images. 
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ethylene dichloride in which small amounts of FeCl. and 
FeCl, had been dissolved prior to irradiation by neutrons 
from a radon-beryllium source. 

Because of the relative scarcity of A®* we suggest that 
it is the active argon isotope formed, the principal reactions 


in the series being 


C8 =e—-+A*, A%*=a+S%, 
Cl** being formed by Cl7+n =Cb»*, 

We wish to thank Dr. L. R. Taussig and the University 
of California Hospital for the donation of the several 
hundred millicuries of radon used in the course of this 
research. The members of the Radiation Laboratory, in 
particular Dr. F. N. D. Kurie and Dr. M. G. White, also 
have our thanks for the advice and assistance they have 
given us. 

W. F. Lippy 
M. D. PETERSON 
\W. M. LATIMER 
Department of Chemistry, 
University of California, 
Berkeley, California, 
August 26, 1935. 


' Libby, Phys. Rev. 46, 196 (1934). 


The Arc Spectrum of Tungsten 


The structure of the arc spectrum of tungsten has been 
studied in this laboratory in connection with the re- 
measurement of the spectrum on plates made last vear by 
the Bureau of Standards. The number of known levels is 


as follows: 
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j-Value Even Odd 
0 3 2 

1 5 29 

2 7 53 

3 9 62 

4 & $1 

5 7 41 

6 2 13 
Total 41 251 


including those levels recently published by F. Poggio! 
which we could verify. Our list includes besides the levels 
published by Laporte in -lfomic Energy States,* those tound 
by Dr. J. E. Mack and kindly communicated by him. 

Approximately 3000 lines are accounted for, but about 
2000 are lines remain unclassified, including about a dozen 
which show self-reversal or furnace of underwater spark 
absorption. Some lines of the latter class arise from rather 
high levels with j-values greater than the lowest 7) term; 
for instance, from the 5-level at 15,069.94 cm™!, which 
was discovered by Laporte and Mack but remains un- 
published. Other important even levels found more recent! 
by the author are: 


4+— 16431.30 5 — 19535.00 5—49354.58 
6— 17008.46 6-— 19648.48 4—51123.08 
+-- 17107.01 4-—47975.45 5—-55333.03 


The last four of these account for some of the most 
intense lines in the spectrum and an odd level has been 
found at 26,676.38 with j-value equal to 5 which combines 
with the first of these to account for an intense line which 
King® lists as possibly showing absorption in the electric 
furnace. 
D. D. Lats 
Marquette University, 
August 24, 1935. 
1F,. Poggio, Anales Sociedad Espafola de Fisica y Quimica 33, 171 


(1935). 3 ee : - 
2 Bacher and Goudsmit, Alfomic Energy States (1932), p. 590. 
3A. S. King, Astrophys. J. 75, 379 (1932). 
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